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SEQUENCE COMPLEMENTARITY OF T2-DNA AND T2-SPECIFIC 
RNA* 


By BENJAMIN D. HALL AND S. SPIEGELMAN 
DEPARTMENTS OF CHEMISTRY AND MICROBIOLOGY, UNIVERSITY OF ILLINOIS 
Communicated by H. E. Carter, December 8, 1960 


Investigations of the functional interrelations among DNA, RNA, and protein 
are most conveniently performed under conditions which limit the synthesis of each 
macromolecular class to a few chemical species. A situation of this type obtains in 
E. coli cells infected with bacteriophage T2. Volkin and Astrachan! examined 
the nature of the RNA synthesized in the T2-coli complex by means of P*?-labeling. 
Estimation of the relative P*? content of the 2’,3’-nucleotides isolated from an 
alkaline hydrolysate led Volkin and Astrachan to deduce that the RNA formed in 
the infected cell possessed an apparent base ratio analogous to that of T2-DNA. 
Subsequently, Volkin? obtained data suggesting that the synthesis of a specific 
RNA is a prerequisite for the intracellular production of bacteriophage. 

Nomura, Hall, and Spiegelman*® confirmed the chservations on the apparent 
base ratios. In addition, they offered independent evidence for the existence of a 
“T2-Specific RNA” by demonstrating that RNA molecules synthesized after in- 
fection differed from the bulk of the £. coli RNA in electrophoretic mobility and 
average sedimentation coefficient. Because the procedures employed (zone electro- 
phoresis and sedimentation) led to a selective separation of T2-specific RNA from 
the normal RNA of E. coli, they open up possibilities of further experiments relevant 
to an understanding of the nature of T2-RNA. 

The fact that ‘““T2-RNA” possesses a base ratio analogous to that of T2-DNA 
is of interest because it suggests that the similarity may go further and extend to a 
detailed correspondence of base sequence. The central issue of the significance 
and meaning of ‘“T2-RNA”’ is whether or not this is in fact the case. <A direct 
attack on this problem by complete sequence determination is technically not 
feasible at the moment. However, some recent findings of Marmur, Doty, et al.4* 
suggest the possibility for an illuminating experiment. These authors demonstrated 
the specific reformation of active double-stranded DNA when heat-denatured 
DNA is subjected to a slow-cooling process. Such reconstitution of the double- 
stranded structure occurs only between DNA strands which originate from the 
same or closely related organisms. Presumably, the specificity requirement for a 
successful union of two strands reflects the need for a perfect, or nc-ur-perfect, 
complementarity of their nucleotide sequences. We have here then a possible 
method for detecting complementary nucleotide sequences. The formation of a 
double-stranded hybrid during a slow cooling of a mixture of two types of poly- 
nucleotide strands can be accepted as evidence for complementarity of the input 
strands. 
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We have used this procedure to examine for complementarity of sequence between 
“T2-RNA” and T2-DNA._ Purified T2-RNA was used in order to provide an 
optimal opportunity for the T2-RNA to combine with its DNA complement, un- 
hindered by non-specific interactions involving irrelevant RNA. Since the hybrid 
would have a lower density than uncombined RNA, a separation of the two might 
be attainable by equilibrium centrifugation in CsCl gradients.6 To insure a 
sensitive and unambiguous detection of the hybrid, should it occur, double labeling 
was used. The T2-RNA was marked with P* and the T2-DNA with H*. Two 
isotopes emitting 8-particles differing in their energies are conveniently assayed in 
each other’s presence in a scintillation spectrometer.’ This device, coupled with 
the use of the swinging-bucket rotor for the equilibrium centrifugation, permits the 
actual isolation of the pertinent fractions along with a ready and certain identifica- 
tion of any hybrids formed. 

The primary purpose of the present paper is to present the results of such experi- 
ments. The data obtained demonstrate that specific complexes are indeed formed 
between “T2-RNA” and its homologous DNA. Their occurrence offers strong pre- 
sumptive evidence for a detailed complementarity of the nucleotide sequences in 


these two macromolecules. 


1. Preparation and denaturation of DNA: Tritiated phage were prepared by the addition of 
H*thymidine to a T2-infected culture of #. coli B. The cells were treated with 5-fluorouracil 
deoxyriboside (0.5 zg/ml) prior to infecton. The phage were purified by treatment of the lysate 
with DN Aase and RN Aase followed by three cycles of high- and low-speed centrifugation. DNA 
was extracted from the purified phage by treatment with sodium dodecyl sulfate followed by 
chloroform-iso-amyl alcohol deproteinization and ethanol precipitation of the DNA.’ This 
preparation will be designated by H*-DNA(T2). DNA from other sources was similarly purified. 

Tritiated FE. coli DNA was prepared from cells of a thymineless mutant (15T~) grown in a 
synthetic medium supplemented with tritiated thymidine. This preparation will be designated by 
H?-DNA (E. coli). 

DNA used for complex formation with RNA was first denatured by heating for 15 minutes at 
95°C in 0.15 M NaCl + 0.01 M sodium citrate (pH 7.8), after which the tube containing the 
DNA was quickly placed in an ice bath. In all cases, the denaturation was carried out at a DNA 
concentration of 130 yg/ce. 

2. Preparation of T2-specific RNA labeled with P*?:  P**-labeled ribosome RNA was obtained 
from E. Coli B grown in synthetic medium, infected with T2 at a multiplicity of 3.8, and labeled 
with 10 millicuries of P%? between three and eight minutes after infection. The infection and 
radioisotope incorporation were done at 37°C in medium C (Roberts e¢ al.) modified to include 
5 gm NaCl, 0.37 gm KCl, and 1 gm casamino acids per liter. The phosphate concentration was 
lowered to 10~-% M and 0.1 M tris (hydroxymethyl) aminomethane (tris), pH 7.3, was used for 
buffering. The number of infective centers and uninfected survivors (2.5°7 ) agreed with the multi- 
plicity of infection. The procedures used for stopping incorporation, washing and disrupting 
cells, and preparing ribosomal RNA were those described previously (Nomura, Hall, and Spiegel- 
man®), 

3. Purification of T2-specific RNA: Enrichment of the ribosomal RNA preparation in its 
content of T2-specific RNA (as judged by an eightfold increase in specific activity of P32) was 
obtained by zone centrifugation through a sucrose gradient. One ml of a 1.5% sucrose solution 
(w/w) + one ml of P*? ribosome RNA soiution (1 mg/ml, 106,000 cpm/ml) were layered, with an 
inverted gradient of RNA on 20 ml of a 2 to 15 per cent sucrose gradient. All solutions were 
0.05 M in KCl and 10-2 M in tris buffer at a pH of 7.3. Following centrifugation for eight hours 
at 25,000 rpm in the SW-25 rotor of the Spinco preparative ultracentrifuge, the contents of the 
tube were removed by dripping through a hole punctured in the bottom of the tube. Fractions of 

.2 ml were collected by drop counting. The ultraviolet absorption at 260 mu and P* content of 
the fractions are shown in Figure 1. The two fractions at the peak of P* activity (corresponding 
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to 18 and 19.8 ml) were used for hybrid formation. These will be referred to in the text as P*?- 
RNA(T2). In the experiments described below, the two fractions exhibited identical properties. 

4. Slow cooling of DNA and RNA: Slow cooling of RNA and DNA was done in solutions 
0.03 M in sodium citrate and 0.3 M in NaCl at a pH of 7.8. An insulated water bath having a 
capacity of 40 liters was used to provide slow cooling as follows: 


Time (hr. Temperature (°C) 
65 ‘ 
5 52 
5 t4 
13.5 36 
24.5 28.5 
30.5 26 


When the bath temperature reached 26°, the tubes containing RNA and DNA were removed and 
brought to a volume of 5.1 ce and a density of 1.74 gm/ce by addition of suitable amounts of 
water and saturated CsCl solution. Twenty-five yg of unlabeled, undenatured DNA were 
added to the solution as a reference density marker. 

5. Separation of RNA from DNA by density-gradient cenirifugation: The solutions of RNA and 
DNA containing CsCl were centrifuged at 
33,000 rpm in the SW-39 rotor at a temper 
ature of 25°C. At the end of each run, 


fractions corresponding to various density 40 3s RNA OF To - INF. E. COLI 
levels in the tube were obtained by piercing * 
the bottom of the tube and collecting drops, 
30 for each fraction. These were diluted 
to a volume of 1.2 ce for measurement of 
ultraviolet absorption and radioisotope 
concentration. 
6. Counting of H®-DNA and P**-DNA 

To an aliquot from each swinging-bucket 


= 

a 

oO 
1000 


fraction 250 yg herring sperm DNA was 
added as carrier. The nucleic acid was 


° 
° 


then precipitated with trichloracetic acid 
(final concentration 10°%) in the cold, Lu L 


collected, and washed on a millipore filter 10 15 
ae eer 


| 

“a 

| 
—A © 

25 


(course, 50 mm dia.). The filter was air- 
dried for one hour and placed in a eylindri- Fig. 1.—Separation of P*-RNA(T2) from E£. Coli 
RNA by sucrose-gradient sedimentation. One ml 
of ribosome RNA solution containing 1 mg RNA 
and 106,000 cpm P*? + 1 ml 1.5% sucrose solution 
(5-phenyloxazolyl)benzene (POPOP) and were layered on 20 ml of a 2-15% sucrose gradient. 
60 mg of 2,5-diphenyloxazole (PPO). P32? Centrifugation: 8 p hours at 25,000 rpm. Cpm 
Packard Tri- shown refer to 0.05 ml fractions of the swinging- 
bucket fractions. 


cal glass v‘al filled with 15 ml redistilled 
toluene containing 1.5 mg of 1,4-bis-2- 


and H® were counted in a 
Carb liquid scintillation counter. 


Separation of T2-Specific RNA in a CsCl Gradient.—It was first necessary to 
establish the conditions required for an adequate separation of T2-specific RNA 
from T2-DNA. Whereas EF. coli ribosome RNA formed a narrow band within two 
days, T-2 RNA, because of its smaller size, required five days to form a band near 
the bottom of the tube. 

Figure 2 shows the result of a five-day run carried out under the conditions speci- 
fied above. The mixture being separated consisted of 6.5 ug of heat-denatured 
H*-DNA(T2), 25 wg of unlabeled and undenatured T2-DNA, and 14 ug of the 
purified P**-RNA(T2). Here, the three nucleic acids were not exposed to a slow- 
cooling operation but were mixed at 25°C, immediately put in the CsCl solution, 
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and centrifuged. It will be noted that there is no appreciable interaction between 
the RNA and DNA as evidenced by the absence of any appreciable overlapping of 
the P**- and H*-containing regions. The small “tail’’ of P*? which extends to the 
top of the tube is presumably a consequence of the low molecular weight of the 
T2-specific RNA. 

Hybrid Formation between Denatured T2-DNA and T2-Specific RNA.—The re- 
sults described in Figure 2 show that CsCl density gradient centrifugation permits 
a clear separation of H*-DNA(T2) from P*?-RNA(T2) and provides, therefore, a 
test for interactions leading to the formation of RNA-DNA hybrids. Any distor- 
tion of the distribution of H*-DNA or P*?-RNA from that observed in Figure 2 which 
leads to regions of overlap between H* and P*? would be indicative of such interac- 
tions. 


-—— [P°*-RNA(T, )+H°-ONA(T, )] 
UNINC CONTROL 1 


| 











ie) 





30 35 
FRACTION NUMBER 
Fic. 2.—Separation of P*-RNA(T2) from H*-DNA(T2) by CsCl- 
gradient centrifugation. A mixture of 6.5 wg heat-denatured H*- 
DNA, 14 ug P?-RNA and 25 ug undenatured, unlabeled T2 DNA 
was made at 25°C, immediately diluted with CsCl, and then centri- 
fuged for five days at 33,000 rpm. 


1. The effect of temperature during slow cooling on hybrid formation: The influence 
of the starting temperature of the slow-cooling process was examined in a number 
of runs. In all cases, the nucleic acid mixture incubated consisted of 6.5 yg of 
heat-denatured H*-DNA(T2) and 14 ug of P#?-RNA(T2). The rate and conditions 
of the cooling were as described earlier. 

Three tubes containing this RNA-DNA mixture were placed in the slow-cooling 
bath at starting temperatures of 65°, 52°, and 40°C respectively. Slow cooling was 
followed by CsCl gradient centrifugation. 

Figure 3 shows the optical density profiles and distributions of H* and P* 
obtained at the three temperatures. Comparison of the profiles of H* and P* 
with those of the control (Fig. 2) shows clearly that in all three cases, slow cooling 
of the DNA and RNA has produced a new peak of P*? approximately centered on 
the band of H* (denatured DNA). This new P*?-containing band must contain an 
RNA-DNA hybrid having approximately the same density as denatured T2-DNA. 
The amount of complex formed on cooling from the three temperatures was the 
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a : +780 
[P??- RNA (Ty) + H9-ONA(T, )] 65°C 1 
| 


— 400 
| 


P?*_ PNA (T,) + H7-DNA(T,) | 52°C 
2 


fe Mion 
35 
+600 | 


FRACTION NUMBER 


Fig. 3.—Formation of DNA-RNA hybrid at various temperatures. 
CsCl-gradient centrifugation analysis. P*-RNA(T2) (14 ug) and 
H*-DNA(T2) (6.5 wg) were mixed in 0.6 ml 0.3 M NaCl and 0.03 M 
Na citrate, (pH 7.8); then the solution was immediately placed in the 
slow-cooling bath. Three identical solutions were made; (a) was 
placed in the bath at 65°, (b) at 52°, and (ec) at 40°C. When the 
bath temperature reached 26°, CsCl and 25 wg T2 IYNA were added to 
each solution; then they were centrifuged for five days at 33,000 rpm. 
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same within experimental error. The three differ slightly in the density of the 
complex relative to DNA, the complex formed at low temperature being apparently 
more dense. This may be explained by the occurrence of partial renaturation of 
the H*-DNA at the higher temperatures. 

2. Requirement for presence of single-stranded DNA during cooling: In order 
to successfully complex with T2-RNA, the molecules of T2-DNA must be present 
in the single-stranded state. This was shown by an experiment in which a mixture 
of native H*-DNA(T2) (13 ug) and P*?-RNA(T2) (15ug) was subjected to slow cool- 
ing, starting from 40°C. No evidence of hybrid formation is observed (Fig. 4). 
In a companion run (a repetition of the experiment of Fig. 3c) with denatured H*- 
DNA(T2), approximately 10 per cent of the P**-RNA was included in the hybrid 


region. 


| [P??-RNA(T,)+H?-ONA (Tz) Native] 
r 40°C 


250-— 
| 











QO Rxncncl nae cloeemee ees 


ie) 5 10 15 20 
FRACTION NUMBER 


Fic. 4.—CsCl-gradient centrifugation of a slowly cooled mixture of 
native H* DNA (T2) with P#?-RNA(T2). 13 wg DNA and 15 wg RNA 
were mixed in 1.2 ml 0.3 M NaCl. 0.03 M Na citrate, slowly cooled 
from 40° to 26°, diluted with saturated CsCl solution, and centrifuged. 


3. Stoichiometry of hybrid formation: Assuming the specific activity of P** to 
be equal in all RNA molecules, one can estimate the amount of RNA which formed 
hybrid. From the data of Figure 3, this would be 1.4 ug RNA. (This figure is a 
maximum value, for some of the RNA which failed to form hybrid may be pre- 
existing 2. coli RNA and, therefore, devoid of P*2.) The amount of DNA in the 
hybrid cannot exceed 6.5 wg, the total amount of denatured T2-DNA present. Be- 
cause the hybrid and denatured T2-DNA have the same density, no more precise 
estimate can be made. From these considerations, it appears probable that the 
ratio of DNA to RNA in the hybrid does not exceed 5. That the complex does in 
fact contain considerably more DNA than RNA is suggested by its density, which 
is very nearly that of T2-DNA. A further indication that the entire 6.5 yg of 
T2-DNA has participated in complex formation is the proportionality observed 
between the amount of hybrid formed and the amount of DNA present when the 
ratio of DNA to RNA is varied. Experiments completely comparable to those 
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described by Figure 3 were carried out with the same concentration of P*?-RNA(T2) 
but with '/; the quantity of denatured T2-DNA. In these cases, the amount of 
P8-RNA found associated with the denatured DNA band was approximately two 
per cent of the input RNA which compares with the average of ten per cent ob- 
served when five times as much T2-DNA is included in the cooling mixture. 

4. On the specificity of the interaction leading to hybrid formation: Having es- 
tablished the existence of the phenomenon and the conditions required for its oc- 
currence, it became of obvious interest to examine the specificity requirements of 
hybrid formation. This was tested by carrying out the cooling process with mix- 
tures of P*-RNA(T2) and denatured DNA from heterologous sources. These in- 
cluded DNA of Ps. aeruginosa, EF. coli, and bacteriophage T5. The DNA of the 
E. coli was labeled with H*-thymidine whereas the others were unlabeled. The 
mixtures of P*?-RNA(T2) and denatured DNA preparations were subjected to a 
52°C slow-cooling incubation under conditions identical to those described for the 
experiment of Figure 3b. Upon completion of the incubation, unlabeled native 
DNA was added to each tube as a density marker. 

Figure 5 gives the optical density profiles and distributions of radioactivities. 
There is a suggestion of a very slight peak of P*? in the DNA region of the mixture 
containing £. coli DNA. It coresponds to !/39 the amount of hybrid produced with 
T2-DNA in a similar experiment. This may reflect the presence of the small 
amount of non-infected cells present during the P**-labeling of the material from 
which the T2-RNA was obtained. However, it is too small to be considered seri- 
ously without further investigation. None of the other heterologous mixtures 
tested yielded detectable amounts of hybrid. It is of interest to note that although 


has the same over-all base ratio as T2, no evidence of interaction between 
-] 


{NA and T5-DNA was observed. 


‘5 
a 
Interpretation of the Results.—The data presented here show that RNA molecules 
synthesized in bacteriophage-infected cells have the ability to form a well-defined 
complex with denatured DNA of the virus. That this interaction is unique to the 
homologous pair is shown by the virtual absence of such complexes when T2-specific 
RNA is slowly cooled with heterologous DNA. The fact that T2-RNA and DNA 
do satisfy the specificity requirement must reflect a correspondence in structure 
between the two. Structural specificity of this order in single polynucleotide 
strands can only reside in definite sequences of nucleotides. We conclude that the 
most likely interrelationship of the nucleotide sequences of T2-DNA and RNA is one 
which is complementary in terms of the scheme of hydrogen bonding proposed by 
Watson and Crick.'® 
Extension to Other Systems.—The bulk of the RNA in E. coli corresponds to the 
18S and 25S'! components of the ribosomes. These are metabolically stable,” 
remain firmly attached to ribosome protein at 10-4 J Mgt*,'’ and have a base 
composition!‘ not related in any obvious way to the DNA of the cells. In addition 
to the lack of correspondence in base ratio, two other reasons can be advanced for 
doubting the suitability of the large ribosomal RNA molecules for directing the 
synthesis of proteins as specified by the genetic material. First, the experiments 
of Riley et. al. suggest that the intermediary between the genome and the protein- 
synthesizing mechanism is metabolically unstable. Second, the formation of the 
larger RNA components is virtually absent in T2-infected cells (cf. Fig. 1) despite 
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Fig. 5.—CsCl-gradient centrifugation of slowly cooled heterologous 
DNA-RNA mixtures) (a) P%-RNA(T2) + H*DNA (E. colt), 
(b) P?-RNA(T2) + DNA (Ps aeruginosa), (c) P*#?-RNA(T2) + DNA 
(T5). 

In each case, 14 ng P*-RNA and 6.5 ug heat-denatured heterol- 
ogous DNA were mixed in 0.6 ml 0.3 M NaCl, 0.03 M sodium citrate 
and slowly cooled from 52° to 26°C. Before CsCl-gradient centri- 
fugation, 25 ug T2 DNA was added to (a) and (ec) and 25 wg Pseudo- 


monas DNA to (b). 
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the fact that they are actively synthesizing a variety of new protein species. 
It seems more likely that the RNA molecules directly concerned with specifying 
protein synthesis in normal cells would have a base ratio corresponding to DNA and 


would possess other properties analogous to those found for T2-specific RNA. Its 
principle characteristics may be summarized as follows:':* (1) a weak linkage 
with the ribosome fractions since it can be broken by dialysis against 10-4. Mg**, 
(2) an active metabolic turnover, (3) an average sedimentation coefficient of about 
8S, (4) a base composition which is closely analogous to its homologous DNA 
(considering thymidine equivalent to uridine and similarly for cytidine and hy- 
droxymethyl cytidine), and (5) a sequence complementary to its homologous 
DNA. 

The detection of the complementary RNA in T2-infected cells was greatly facili- 
tated by the fact that the larger ribosomal components are not synthesized. In- 
deed, it would appear as if RNA synthesis in the T2-coli complex is largely confined 
to the class which is complementary to DNA. This advantage is not present in 
uninfected cells. Consequently the search for normal complementary RNA will be 
technically more difficult. That it is nevertheless feasible is suggested by the 
experiments of Yéas and Vincent’ with yeast. These authors used P*? in a manner 
comparable to the procedures of Volkin and Astrachan and, despite surprisingly 
long pulses, they were able to detect the formation of a fraction with a high meta- 
bolic turnover and possessing a base composition analogous to yeast DNA. 

Ultimately, attempts at establishing the presence in normal cells of RNA com- 
plementary to the genetic material will require that it be separated from the other 
RNA components. If all complementary RNA molecules possess physical chemical 
characteristics analogous to those of T2-specific RNA, the same methods which 
effected a successful isolation in this case may well serve in others. Once isolated, 
sequential complementarity to relevant DNA can be examined by the methods 
described above. 

Some Implications of Complementary RN A.—An increasing amount of attention 
is currently being focused on the possibility of forming hybrid helical complexes 
composed of DNA and RNA strands. Interest in this originates quite naturally 
from its obvious implication for translating the genetic information coded in DNA 
to a functional RNA complement. Previous experiments": '* had already demon- 
strated that paired helices were generated in mixtures of the synthetic polyribo- 
this has been ex- 


20 


nucleotides of uridylate and adenylate. More recently,'® 
tended to combinations involving synthetic polydeoxyribonucleotides and _ poly- 
ribonucleotides. It is of some interest that the experiments reported in the present 
paper lend support to the concepts underlying such model experiments by exhibiting 
hybrid formation between natural polynucleotides which are complementary and 
biologically related. 

The demonstration of sequence complementarity between homologous DNA and 
RNA is happily consistent with an attractively simple mechanism of informational 
RNA synthesis in which a single strand of DNA acts as a template for the poly- 
merization of a complementary RNA strand. 

Summary.—Experiments are described showing specific complex formation 
between single-stranded T2-DNA and the RNA synthesized subsequent to infec- 
tion of F. coli with bacteriophage T2. No such hybrid formation is observed with 
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heterologous DNA even if it has the same over-all base composition as T2-DNA. 
It is concluded that T2-DNA and T2-specific RNA form hybrids because they 
possess complementary nucleotide sequences. The generality of the existence of 
complementary RNA and its possible role as a carrier of information from the gen- 
etic material to the site of protein synthesis is briefly discussed. 


We are indebted to Drs. Noboru Sueoka and John Drake for providing the T5 
DNA and Pseudomonas DNA used in this work and to Miss Cherith Watson for 
her assistance in performing the experiments. 


* This investigation was aided by grants in aid from the U.S. Public Health Service, National 
Science Foundation, and the Office of Naval Research. 

1 Volkin, E., and L. Astrachan, Virology, 2, 149 (1956). 

2 Volkin, E., these PROCEEDINGS, 46, 1336 (1960). 

’ Nomura, M., B. D. Hall, and 8. Spiegelman, J. Mol. Biol., 2, 306 (1960). 

* Marmur, J., and D. Lane, these PROCEEDINGs, 46, 453 (1960). 

5 Doty, P., J. Marmur, J. Eigner, and C. Schildkraut, ibid., 46, 461 (1960). 

6 Meselson, M., F. W. Stahl, and J. Vinograd, tbid., 43, 581 (1957). 

7 Okita, G. T., J. J. Kabara, F. Richardson, and G. V. Le Roy, Nucleonics, 15, 111 (1957). 

8’ Marmur, J., personal communication. 

® Roberts, R. B., P. H. Abelson, D. B. Cowie, E. T. Bolton, and R. J. Britten, “Studies of 
Biosynthesis in Escherichia coli,’ Carnegie Institution of Washington, Publication #607 (1957), 
p. 5. 

1 Watson, J. D., and F. H. C. Crick, Nature, 171, 964 (1953). 

1! Kurland, C. G., J. Mol. Biol., 2, 83 (1960). 

12 Davern, C. I., and M. Meselson, thid., 2, 153 (1960). 

13 Tissiéres, A., J. D. Watson, D. Schlessinger, and B. R. Hollingworth, zhid., 1, 221 (1959). 

14 Spahr, P. F., and A. Tissiéres, ibid., 1, 237 (1959). 

% Riley, M., A. B. Pardee, F. Jacob, and J. Monod, ibid., 2, 216 (1960). 

1 Yéas M., and W. 8. Vincent, these PROCEEDINGs, 46, 804 (1960). 

17 Warner, R. C., Fed, Proc., 15, 379 (1956). 

18 Rich, A., and D. R. Davies, J. Am. Chem. Soc., 78, 3548 (1956). 

19 Rich, A., these PROCEEDINGS, 46, 1044 (1960). 

2” Schildkraut, C. L., J. Marmur, J. R. Fresco, and P. Doty, J. Biol. Chem. (in press). 


FORMATION OF HETEROZYGOTES BY ANNEALING A MIXTURE 
OF TRANSFORMING DNAS* 
By Rocer M. Herriorr 
DEPARTMENT OF BIOCHEMISTRY, JOHNS HOPKINS SCHOOL OF HYGIENE AND PUBLIC HEALTH 
Communicated by Paul Doty, December 27, 1960 
An important discovery was made when Marmur and Lane! observed that the 
genetic properties of pneumococeal transforming DNA which had been lost through 
heating to 100°C were restored by slow cooling. The loss of function was con- 
sidered to be due to denaturation, a rupture of the inter-strand hydrogen bonds 
followed by a separation of the strands,? and Marmur and Lane suggested that 


restoration occurred by a highly specific pairing of complementary strands—a 


process favored by slow cooling. Since recovery of genetic function of a small 
quantity of transforming DNA was enhanced by the addition of unmarked but 
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homologous DNA from recipient pneumococci and not with heterologous DNA, 
the authors suggested that a marked strand might pair with the complementary 
strand of an unmarked but otherwise equal molecule and form a mixed molecule 
or “heterozygote.” Physical chemical evidence in support of this suggestion was 
presented by Doty, Marmur, Figner, and Schildkraut.* : 

To test their suggested mechanism, we decided to attempt the formation of 
doubly marked heterozygotes by heating 





and annealing a mixture of differently marked 
transforming DNAs. Markers known to link 
were chosen, for they would represent allelic 
or similar molecules and their strands would 
be expected to complement each other to a 
high degree. Although segregation of mark- 


IONS 


ers might well have obscured detection of the 
heterozygotes, the present results show that 
a physical unit carrying both markers was 
formed by annealing a heated mixture of 
these two DNAs. 


Materials and Methods.—In general the methods 


TRANSFORMAT 


sc 


used were described previously.':* A detailed ac 


OF 


count of all the methods will appear elsewhere. 
Genetically marked DNAs: Four differently 
marked DNAs were employed in this work. They 


ON 


were isolated from cell stocks of Hemophilus influ 


na 
rs) 
a 
x 
u 


enaze resistant to 250 ug/ml of streptomycin (8S), 2.5 
ug/ml cathomycin (C) (novobiocin), 6-10 wg/ml of 
erythromycin (EF), and finally a stock resistant to 
both streptomycin and cathomycin (S-C). The 








S-C stock organism was obtained by transforming 





competent streptomycin-resistant cells with C DNA. 
3.4 = «(C5 


It could also have been prepared by the reciprocal : 
DILUTION OF DNA 


transformation. Each DNA preparation, after first 
lysing the organisms at pH 11, was purified by the Fig. l.- Effect of concentration of 
DNA before and after heating on the 
number of double (S-C) transformations. 
cy Pama cae goad The materials, concentrations, and treat- 
cipitation from 1.5 M NaCl-35 per cent ethanol ment for this study have been described 
solution. The E260/E230 ratio of the preparations in the section on Methods. After an ini- 
tial '/i9) dilution of the DNA, samples 
were diluted further as noted in the ab- 
J di ? ; . i 2 scissa of Figure 1, after which equal ali- 
of purified DNA dissolved in 0.5 M NaCl-0.01 quots of these dilutions were mixed with 
M Na; citrate was placed in a 15 ml screw-capped competent cells in the usual manner. 


tube and heated in four liters of water at the 


chloroform-octanol method of protein removal fol- 
lowed by digestion with RNAase and repeated pre- 


ranged from 2.3 to 2.4. 
Heating and annealing: Ingeneral, a 1-ml sample 


indicated temperature, after which the tube was transferred rapidly to four liters of water at 90°C, 
which was then allowed to cool without stirring. For those experiments at temperatures above 
100°C, the 1-ml samples were sealed into glass tubes and then heated in a glycol bath. When the 
bath temperature was 100°C, the temperature in the sample rose to 95°C or higher in one minute. 
In annealing, the temperature dropped from 90° to 60°C in 100 minutes and to 37°C in another 
100 minutes. Cooling from 100°C instead of 90°C gave similar but somewhat lower recoveries 
of activity. 

Transformation procedure: Unless otherwise noted, the transformation reaction was carried 
out by mixing 0.1-ml competent receptor cells (1.5 & 10%/ml) with 2.8 ml saline and 0.1 ml of 0.2 
ug/ml DNA, a 1/100 dilution of the annealed DNA solution. After thirty minutes’ incubation, 
an aliquot was diluted appropriately and plated with nutrient-agar after which it was incubated 





148 BIOCHEMISTRY: R. M. HERRIOTT Proc. N. A. S 


two hours for development of antibiotic resistance. The plates were then covered with an equal 
volume of nutrient-agar containing the antibiotic at double strength. The plates were incubated 
twenty-four to forty-eight hours and the resulting colonies counted. Duplicate plates agreed 


within +20 per cent but duplicate experiments were more variable. 


Experimental Results.—Doubly marked transformants from a mixture of singly 
marked units: The results of an experiment in which a mixture of S and C DNAs 
was heated to 100°C and then annealed are shown in (c) of Table 1. As controls, 
four additional tubes were included. The first (a) was an unheated sample; (b) 


TABLE 1 
HETEROZYGOTES FORMED BY ANNEALING MIXTURES OF TRANSFORMING DNAs 


Materials = A as prepared from separate stocks of cells resistant to different antibiotics. S = resistance to 
streptomycin, C resistance to cathomycin (novobiocin). Purified DNA was diluted to contain 20 ug/ml in 
0.5 M NaCl 0.01 VM Nas citrate. 

1 ml DNA in a capped tube was placed in boiling water for 3 minutes, after which it was either 
Samples were then diluted '/500 (5-fold greater than described in the section on 
Aliquots were then screened 


Procedure 
chilled in ice water or annealed. 
Vethods), after which 0.1 ml was mixed with competent cells in the usual manner 


for transformations to 8, C, and S-C. 
> 
S-C cells 
expected 
Titer in (product of 
Temperature of Cooling reaction frequencies 
Sample heating procedure Marker Dilution Count mixture of singles) 


a 25°C S Lx 16" 26 24 Xe 
2K i1¢ 233 ix 1¢ 
Undiluted ( i) 
Slow Ss 1 x 10? 
5 X 10? 
Indiluted 
100°C Slow 5 5 X 10! 
L x AG 
Undiluted 


100°C Rapid 5 2 X 10! 


5 X 10! 
Undiluted 
100°C Slow S 2 xX 10! 


Each DNA heated and cooled 5 X 10! 
separately, then mixed. 
Undiluted 


was the mixture heated to 90°C and then cooled slowly; in (d), the mixture was 
heated and chilled rapidly instead of annealing; and finally in (e), the two DNAs, 
S and C, were separately heated to 100°C and cooled before they were mixed. The 
number of double transformants expected from the random uptake of the indi- 
vidual markers was calculated in Table 1 from the product of the individual fre- 
quencies. In the above system, where the cell concentration was 5 X 107/ml, the 
expected number of doubles equals the product of the S and C transformations 


divided by the number of cells. 





Vou. 47, 1961 BIOCHEMISTRY: R. M. HERRIOTT 


An examination of the results in Table | shows the following: 

1. Heating the mixture of S and C DNAs to 100°C followed by a slow cooling 
(c) produced many more S-C transformations, i.e., to both markers, than expected 
from a random distribution of 8 and C units. 

2. The number of S-C transformations from the annealed mixture was greater 
than observed in the unheated mixture (compare ¢ and a); i.e., there was a net in- 
crease. 

3. Formation of the S-C unit required a temperature above 90°C, for at this 
temperature no unexpected double transformations were observed (0). 

4. Slow cooling is important to formation of this S-C unit for there were few if 
any unexpected doubles when the heated sample was chilled suddenly (d). 

5. S-C units do not form if the DNAs are mixed after heating and annealing. 
They must anneal together (e). 

The Nature of the New Agent Producing Double Transformations.—The effect of 
concentration of annealed DN A on the number of double transformations: The results 
in Table 1 show that the increase in doubles after heating and annealing could not 
have been produced by the random uptake of residually active individual units, 
but conceivably some type of reactivation or nonrandom uptake of these markers 
following the heat treatment might produce the observed effect, although the 
former seems unlikely in the light of the net increase in doubles over those present 
in the unheated sample. Except for an extreme case of nonrandom uptake in which 


the uptake of one marker produced a strong preference for the other, it would be 


expected that doubles produced by two physical units could be distinguished from 
those produced by a single unit carrying both markers by the manner in which the 
double transformations varied with DNA concentration. Doubles produced by 
two units would vary with the square of the DNA concentration, whereas those 
produced by a single unit would vary directly. The results of such a test are 
shown in Figure 1, where the double transformations obtained before and after 
annealing are plotted against the dilution of DNA. The solid lines are the curves 
expected if before annealing, the doubles arise from two independent units and 
after annealing, from a single physical unit. lor comparison and to bring out the 
difference in the two systems, the fraction of the double transformations is plotted 
on a logarithmic scale against the simple dilution of the DNA. 

The results in Figure 1 show clearly that annealing changed the agents respon- 
sible for double transformations. After annealing, dilution of the DNA produced 
the effect expected if the doubles were produced by a single unit. Titers of the 
doubles obtained after heating were adjusted when the contribution from random 
doubles was greater than 10 per cent. The results are strong evidence that a single 
physical unit carrying both genetic markers (a heterozygote) has been produced by 
heating and annealing a mixture of these two DNAs. 

Effect of ultraviolet radiation on the heterozygote: Ultraviolet radiation was used to 
detect a difference between the hetero- and homozygous doubly marked units. 
The results are shown in Figure 2. The difference is clear but difficult to interpret 
other than to state the obvious fact that the heterozygous unit was less sensitive to 
ultraviolet radiation than was the homozygous S-C. 

The mechanism of heterozygote formation: Most of the conditions favoring 
heterozygote formation in the highest concentration were those which Marmur and 
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Lane! reported for maximal restoration of transforming activity. Thus, about 20 
ug DNA/ml in 0.5 MW NaCl heated to 100°C for a few minutes and then cooled 
slowly from 90°C appeared to be optimal, but some heterozygotes were formed in 
0.05 M NaCl. In studying the mechanism of formation, increasing quantities of 
one marker were mixed with decreasing quantities of the other marker, so the total 
DNA was constant, and after heating and annealing, the mixtures were scored for 
the number of single and double transformations. Tligure 3 shows a plot of the 
fraction of the total of each marker present as a double transformant in the various 

mixtures of DNA. If the annealing 





process is highly efficient and the 
mechanism of reformation is the 
specific strand pairing suggested by 
Marmur and Lane, then, when a 
small quantity of one marker (3) 
was heated with a relatively large 
quantity of the second (C), it would 
be expected that a high proportion 
of the S would be found in the het- 
erozygotes, for there would be more 
C strands to unite with S strands 
than there were complementary 8 


SC MARKER 


SURVIVING 


strands. As seen in Figure 3, the 


OF 


proportion of the total S as S-¢ 
(SC/S) did increase as the fraction 
of S in the initial mixture of DNA 
decreased, but even in the extreme 
case of 0.5/19.5 or 0.1/19.9, the 
proportion of S appearing as heter- 


oO 
= 
S) 
a 
« 
u 


ozygote did not exceed 0.05. In a 








comparable way, the proportion of 





; 30 40 5 : ‘ . ‘ 
C appearing as heterozygote did 

U.V. DOSE (sec.) 9 mo . i 
; oo . not exceed 0.02. The experiment 
Fig. 2.—A comparison of the sensitivity of : : 
heterozygote and natural (S-C) DNA to ultraviolet has been repeated with many varla- 
radiation.—Doubly marked S-C DNA and a ti yreedi we fie 
ee , and ions without exceeding these fig- 
heated and annealed mixture of S and C DNAs ; < ts 
each at 1 wg/mlin0.15 MV saline were exposed to UV ures, which suggests that perhaps 
light from at 15-watt G.E. germicidal lamp at a dis- thev are upper limits. Neither the 
tance of 388 em. This delivers a dose of 29 ergs ; 

mm2/sece to the solutions. 


linear effect in Figure 3 nor the low 

recovery of heterozygotes is in 
keeping with the simple bimolecular reaction such as specific strand pairing. In 
view of the relatively high (often 50 per cent or more) recovery of the individual 
markers it appears that an unrecognized feature of this system must be control- 
ling heterozygote formation. Several possible explanations will now be considered, 
some of which permit experimental testing. 

A. The number of heterozygotes in the annealed solution may be greater than 
observed with the biological assay responsible for the reduced number. Thus, 
segregation (separation of the genetic units into different daughter cells) or re- 
combinational events acting on such an unusual unit as the heterozygote may be 
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responsible for the low recovery. A clonal analysis should throw light on the 
possibility of segregation. 

B., Perhaps due to their complete complementarity, pure singles reform during 
cooling at a slightly higher temperature than do the heterozygotes, so that when 
the temperature falls to that permitting heterozygotes to form, the number of 
separated strands is lower than expected. On this basis, the heterozygotes should 
melt at a lower temperature than the individual markers on reheating. Experi- 
ments thus far on this point are not decisive. This explanation is weak, however, 
in view of the results in Figure 3. 

C. Under the conditions of heating, the strands of most of the molecules may 
not completely separate despite the nearly complete loss of function and a rise in 
chromicity on chilling. However, extending the time of heating or, as shown in 
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Fig. 3.-—-Proportion of marker incorporated into the heterozygote as the 
proportion of markers in the mixture prior to heating is varied. Mix- 
tures of Sand C DNAs at a constant concentration of 20 ug/ml but different 
proportions were heated and annealed in 0.5 M NaCl-0.01 Nay citrate after 
which they were assayed for 8, C, and S-C markers. 


ligure 4, raising the temperature above 100°C to separate the strands of more mole- 
cules failed to increase heterozygote formation, so that this explanation lacks direct 
support, 

D. If the strand complementarity requirements for reformation are less specific 
than originally proposed and strands from nonallelic molecules have sufficient com- 
plementarity (points of hydrogen bonding) to pair with the marked strands and 
thereby prevent reunion of the more complementary strands and also to prevent 
collapse of the marked strands on cooling, this would explain the low recovery of 
heterozygotes in the face of a high recovery of individual markers. //emophilus 
influenzae contains about 50 molecules of DNA per cell, only one of which appears 
to carry a particular marker, so that if all strands are capable of some degree of 
complimentation, a lower limit of a few per cent in recovery of heterozygotes would 
be expected. If a relatively random pairing of strands is the case, then a recovery 
of twice the original number of individual markers would be expected after anneal- 





152 BIOCHEMISTRY: R. M. HERRIOTT Proc. N. A.S 


ing, which is not observed. However, this assumes that the probability of incor- 
poration of a marker is the same for molecules having one or both strands carrying 
the marker. 

If the pairing is not highly specific, other doubly marked heterozygotes would be 
expected. A few experiments with mixtures of streptomycin- and erythromycin- 
resistance markers failed to reveal any doubly marked units over those expected 
from random uptake of singles, but at least one explanation for this failure may be 
suggested. The S and C markers under normal circumstances link or recombine 
together soon after entering the cell,> whereas S and E do not,® and in the absence 
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Fic. 4.—-Effect of high temperature on the recovery of 
heterozygotes. DNA mixtures at the usual concentration 
sealed in glass tubing were heated in a glycol bath for 5 
minutes at the indicated temperature, after which they 
were transferred to a bath at 90°C and annealed. 


of linkage, segregation may obscure the presence of the heterozygote. A clonal 
analysis of an annealed mixture of S + E DNAs would help clarify the possibility 
that nonallelic strands can pair up. 

EF. If, as Doty et al.* suggest, a fraction of the strands develop breaks at the 
high temperatures, then, on annealing, a single unbroken strand could rescue, by 
pairing up, segments from both S and C strands. In such a model, the low recovery 
of heterozygotes in the face of high recoveries of individual markers is accounted 
for, and more important, it provides a possible mechanism for linking the two 


markers. 
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F. Finally, the possibility that the heterozygote consists of two singly marked 
molecules fused together through partial pairing making a unit composed of four 
original strands cannot as yet be excluded. 

Annealed samples were shaken with 80 per cent phenol, treated with chymotryp- 
sin and shaken with chloroform-octanol to remove any protein which might have 
bound them together, but none of these treatments altered the ratio of doubles to 
singles. 

With so many possible explanations for the low recovery of heterozygotes, some 
decisive experiments must precede further consideration of the mechanism. One 
additional result may be noted, however, in closing. Heated mixtures of markers 
which had been chilled rapidly (such as din Table 1), and which consequently showed 
a low level of both markers, returned to 25-50 per cent of the initial transforming 
activity on reheating and annealing just as Marmur and Lane had found. More 
important for this discussion, heterozygotes which were absent from the chilled 
mixture formed during the annealing of the reheated DNAs. This shows that the 
rapid chilling which causes collapse of the separated strands,* does not, however, 
inflict irreparable damage on either marker or their capacity to form heterozygotes. 

Summary. Heat denaturation of a mixture of genetically different transforming 
DNAs from different stocks of Hemophilus influenzae followed by the annealing 
treatment recommended by Marmur and Lane led to the formation of heterozy- 
gotes, physical units carrying both genetic markers. A number of possible mech- 
anisms for their formation have been considered. 

The author is grateful to Drs. Sol H. Goodgal, Claud 8. Rupert, and C. A. Thomas, Jr., for their 
criticisms and suggestions. He also wishes to acknowledge the excellent technical assistance of 
Miss Stella Mayorga-Nestler. 
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THE KINETICS OF CARRIER-MEDIATED ACTIVE TRANSPORT 
OF AMINO ACIDS* 
By Joun A. JACQUEZ 


DIVISION OF EXPERIMENTAL CHEMOTHERAPY, SLOAN-KETTERING INSTITUTE FOR CANCER RESEARCH, 
AND THE SLOAN-KETTERING DIVISION, CORNELL UNIVERSITY GRADUATE SCHOOL OF MEDICAL SCIENCES 


Communicated by Joseph O. Hirschfelder, December 5, 1960 


Amino acids are concentrated in many types of cells by an active transport 
process. However, it is only with the cells of Ehrlich ascites carcinoma that 


sufficient data have accumulated for us to consider setting up models of this process. 


We proceed by first setting up some general models of active transport systems 
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and use the stationary state approximation to derive initial fluxes. These can 
then be compared with the experimental data on amino acid active transport. 

Experimental Characteristics of the Amino—Active Transport System in Ehrlich 
Ascites Cells.—Let. us first briefly review the major experimental features of amino 
acid transport in Ehrlich ascites cells. 

Initial uptake velocity: The initial flux in a 1- to 2-min incubation increases 
rapidly as the extracellular concentration of amino acid rises and asymptotically 
approaches a linear function of the average gradient across the cell membrane. 
Thus, the initial flux may be written F = k,-(c, — c;) + Fr, where k, is the per- 
meability constant for the cell membrane, (c, — c¢,) is the average gradient across 
the cell membrane during the period of incubation, and Fr is the net active trans- 
port. The latter shows saturation behavior as the extracellular concentration 
increases. 

The steady state: After a two-hour incubation, the intracellular concentration 
of amino acid is fairly constant. A graph of the intracellular versus the extra- 
cellular concentration after two hours of incubation shows a high value for the 
ratio c,/c, at lowe,. This falls, approaching 1 as c, increases. However, the data 
from steady-state studies at high extracellular concentrations must be used with 
caution, because there is evidence of some cell damage. 

Competition among amino acids: Although many of the amino acids compete 
with one another for active transport, in some cases the initial flux of one amino 
acid is increased when another amino acid is present in roughly equimolar amounts. 
L-tryptophan and O-diazoacetyl-L-serine (azaserine) represent such a combina- 
tion. When each is present at about 1 mM/L in the extracellular phase, the 
initial flux of L-tryptophan is increased and that of azaserine is decreased as 
compared to the fluxes when each is present alone at the same concentration. 
However, if the extracellular concentration of azaserine is increased to about 
5 mM/L, the initial flux of L-tryptophan is decreased. Similar results are found 
if L-tryptophan is incubated with the Ehrlich ascites cells in the presence of L- 
histidine, L-2,4-diaminobutyric acid, and L-leucine. 

Exchange diffusion: Heinz! reported that the initial uptake flux of C-labeled 
glycine was markedly increased when the cells had first been loaded with unlabeled 
glycine. Subsequently, Heinz and Walsh? reported a hetero-exchange diffusion 
in which the initial flux of alanine was increased by preloading the cells with 
glycine. We have shown that the initial flux of L-tryptophan is slightly decreased 


after preloading the cells with glycine. However, if the ascites cells are first loaded 
with L-tryptophan or azaserine, the initial flux of the other is markedly increased. 

Effect of inhibitors: Many of the metabolic inhibitors markedly impair active 
transport of the amino acids.*~* Dinitrophenol and cyanide markedly inhibit 
the initial flux of L-tryptophan active transport but only slightly inhibit the much 
larger initial flux of L-tryptophan when the cells are preloaded with azaserine.® 


Thus, although active transport is directly dependent on cell metabolism, exchange 
diffusion is not. 

Models of Active Transport.—Two major types of models have been suggested 
to explain transport across a cell membrane. These may be labeled adsorption- 
type and carrier-type models. The former is characterized by the presence of 
fixed sites which can combine with a substrate in the extracellular or intracellular 
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phase and can release the substrate into one or both phases. This type of model 
may be generalized to include a series of adsorption sites, a substrate being passed 
from one adsorption site to the next. Carrier models, on the other hand, are 
characterized by the presence of a mobile carrier which can combine with free 
substrate at one surface of the cell membrane to form a carrier-substrate complex 
capable of traversing the cell membrane by diffusion or by an active process and 
of releasing the substrate at the other surface of the cell membrane. In this case, 
the assumption is usually made that within the cell membrane the equilibrium 
constant for the reaction between substrate and carrier is such as to make dissocia- 
tion of carrier substrate complex negligible. 

Adsorption transport models: It is well to briefly consider the adsorption trans- 
port models first if for no other reason than to justify our dismissal of them as 
inadequate models. Detailed mathematical arguments will not be presented 
because the adsorption models are simpler than the carrier transport models and 
essentially the same type of treatment (but simpler) can be used with them as will 
be used with the carrier models. 

Irreversible adsorption transport models: This model is illustrated in Figure 1. 
The area between the two vertical lines represents the cell 
membrane; e and 7 are the extracellular and intracellular 


phases respectively. All of the adsorption sites C are avail- S$ +| 


able to extracellular amino acid for formation of the CS com- 
plex. This is of interest for two reasons: The differential 
equations describing this model can be solved explicitly under 
the assumption of constant extracellular concentration and 
constant cell volume, conditions which can be approximated 
quite closely in amino acid uptake experiments up to extra- 
cellular concentrations of 10-15 mM _/L. Secondly, it is sur- 
prising how closely one can fit the kinetic data on uptake 
of a single amino acid with this model. It is obviously un- 
realistic, for it cannot give any exchange diffusion because 
of its irreversibility. 

Reversible adsorption transport: The previous model may 
be extended somewhat by assuming that all steps are rever- 
sible. This is essentially the model used byLeFevre’ in Fic. 1.—Schematic 


discussing sugar transport in erythrocytes. This model can pec 
give exchange diffusion in the sense that the adsorption site 

can form a complex with an amino acid from one side of the membrane and release 
it on the other side in exchange for another amino acid. However, if two differ- 
ent amino acids are present, one on each side of the membrane, one can show that in 
this model they will compete for the adsorption site. Thus, it cannot give an 
increased net uptake flux of an amino acid after preloading the cell with another 
amino acid. This is obvious from the assumption that the adsorption sites are 
equally accessible to substrates in either phase. 

The ‘bucket brigade” model: Danielli,s among others, has considered models 
in which the substrate is passed along a chain of non-mobile adsorption sites in 
the membrane or in a pore passing through the membrane. This too can give 
exchange diffusion in the same sense that the previous model can. But again 
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it cannot give an increased uptake flux of one amino acid after preloading the 
cell with the same or another amino acid. If we consider the 7th adsorption site 
at the start of such a process, then amino acid (1) from site (¢ — 1) must compete 
for site 7 with amino acid (2) from site (i + 1). 

Carrier transport models: The carrier transport models are generally more 
complicated than the adsorption transport models. Although one can consider 
extremely complicated carrier models involving various enzymatic reactions,° 
this seems needlessly complex at this stage of our knowledge. It is important to 
consider the implications of the various types of carrier models at the simplest 
possible level. 

The carrier active transport models may be classified on the basis of the nature 
of the linkage between the transport system and cellular metabolism. We differen- 
tiate four main types: 

1. Active transport of free carrier. The linkage to cell metabolism is through 
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Fic. 2.—Schematie of generalized model of carrier active transport. 
e = extracellular phase; i = intracellular phase; S; = substrate 1; 
C = carrier; CS, carrier-substrate complex. 








one or more reactions which speed up the movement of free carrier from the inner 
towards the outer surface of the cell membrane. 

2. Active transport of carrier-substrate complex. ‘Transfer of carrier-substrate 
complex from the outer to the inner surface of the cell membrane is speeded up 
by linkage to cellular metabolism. 

3. Association reaction model. The rate of formation of substrate-carrier 
complex at the outer surface of the cell membrane is increased by linkage to cell 
metabolism. 

1. Dissociation reaction model. The rate of dissociation of substrate-carrier 
complex at the inner surface of the cell membrane is increased by metabolically 
linked reactions. 

One cannot of course rule out the possibility of combinations of these mechanisms. 

Widdas” has used simplified versions of a carrier model in discussing sugar 
transport. Some of the carrier models proposed by Rosenberg and Wilbrandt® 
are similar to those above. However, they introduce the simplifying assumptions 
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(1) that the total concentration of carrier is equal on the two sides of the membrane, 
and (2) that the diffusion constants in the membrane are the same for free carrier 
and carrier-substrate complex. Neither of these is a physically reasonable as- 
sumption. Finally, it should be pointed out that the first two adsorption transport 
models can be viewed as special cases of the carrier models in which the diffusion 
rates in the membrane are much greater than the rates of reaction of substrate 


with carrier. 
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Fic. 3.—Schematie of model of carrier active transport with two 
substrates. 
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A simplified general model: We proceed by setting up a simplified generalize 
model which includes the above types as special cases. We neglect the amount of 
carrier or carrier-substrate complex which is in transit between the two surface 
phases of the cell membrane. Figure 2 is a schematic of the model when only 
one substrate is present. Figure 3 is the same model when two substrates are 
present. 

Notation: 


C = total amount of carrier in mM_/sq cm surface area. 
concentration of free carrier at outer surface of cell membrane in units 
mM/sq em. 
concentration of free carrier at inner surface of cell membrane. 
concentration of CS; at outer surface. 
concentration of it inner surface. 
concentration of CS». at outer surface. 
concentration at inner surface. 
concentration of substrates 1 and 2 respectively in extracellular 
phase (mM /kg water). 
Ca, Cr concentration of substrates | and 2 respectively in intracellular phase. 
kn, h 


A surface area of cells. 


‘pe permeability constant for free substrate 1 and 2 respectively. 


V volume of extracellular phase. 
V intracellular water, assumed constant. 
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The rate constants for the various reactions and transfers are indicated in 
Figures 2 and 3. 

In order to keep the models as simple as possible, we assume symmetry in rate 
and diffusion constants between the two sides of the cell membrane except for the 
portion of the transport system which is linked to cell metabolism. Although the 
linkage to cell metabolism may be quite complicated, we introduce it simply as an 
increase in the pertinent rate constant. Thus, the various types of models are 
characterized by the following conditions on the rate constants. 

1. Active transport of free carrier: 


ky > k 1 


Active transport of carrier-substrate complex: 


Qi = do, i; = Bo, ko > k 
b» 


= VS; 4 = ks Sk 


ky -- k 
Association-reaction model: 


ky = ae 


k 
71 
Bo > Bi, 6» > 6 . 


It should be emphasized that the above conditions do not violate the principle 
of detailed balancing because the rate processes shown in the model are not a closed 
system. It is assumed that the rate processes which are linked to cellular metab- 
olism are linked via reactions which are not shown in the model but which provide 
the free energy to increase the rate of the process involved. The effect of this 
linkage is included at the simplest possible level as an increase in the rate constant 
of the process which is assumed linked to cellular metabolism. The effect of 
inhibitors of cellular metabolism would be to decrease this rate constant; and in 
complete inhibition of cellular metabolism, the inequality signs in the above 
cases would be replaced by equality signs, as would be required by the principle 
of detailed balancing. 

The mathematical model: We assume that V and V, remain constant. The 
equations describing the model when two amino acids are present are then: 


Co Let eit Yo t yi t Ze + 2. (1) 


—(kK_y + aca + Vile) kya; + BiYe + 512. 
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dy, 
dt 


+ Bi)y. +h 


dy 


OC ty + hoy, — th 
dt 


dz, 
dt 


dz 
dt 


dt 


(10) 


The equations applicable when only substrate ‘1’ is present can be obtained by 


setting z. = 2, C2 Co = 0. 

Such a system of ‘ordin: ary differential equations presents great analytical diffi- 
culties because of the non-linear terms. Furthermore, the initial conditions on 
x, and x; are unknown except for the constraint 1, + 2; = Co at ft 0. It is 
difficult to see how sensitive the solutions are to changes in the initial values of 
x, and x; However, it is more important to determine whether the carrier models 
could give the same type of phenomena as are found experimentally in active 
transport and exchange diffusion experiments. The detailed kinetics are of sec- 
ondary importance, particularly since these models are quite simplified. We 
consider only the carrier-mediated flux, since the flux due to simple permeability 
is usually quite small in comparison to carrier-mediated flux. Let F’r(1) represent 
the active transport flux when only substrate 1 is added to the cell suspension and 
let Fy(1) represent the carrier-mediated thux of substrate 1 when the ‘cells have 
been preloaded with substrate 2. We wish to know whether the carrier models 
can give F'y(1) > F7(1) and what the effects of inhibitors are on Fy(1) — Fr(1) 
and Fy(1)/Fr(1) when Fy(1) > Fr(1) 

Although the equations cannot be solved as they are, they can be solved for the 
stationary-state fluxes. In analogy with the derivation of the Michaelis-Menten 
equations of enzyme kinetics, we obtain the stationary-state fluxes for conditions 
corresponding to those at the start of the uptake experiments. This is equivalent 
to performing the hypothetical experiments of keeping ¢ and cj» constant and 
C2 = Ca = 0 and solving for the stationary state. Furthermore, it should be 
stressed that except for possibly some special initial conditions, these fluxes repre- 
sent upper bounds for the initial fluxes in the experiments. Under these conditions, 
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the stationary-state equations corresponding to equations (1) to (10) reduce to 
equations (11) to (16). 


Zotityet vit ee te = Ce (11) 
—(k_y + ayea)e, + hiv: + Bry. + &2 = 0. (12) 
ANCL — (ke + Bilye + key; = 0. (13) 
koye — (k-2 + Bay; = 0. (14) 

—(ks + 6:)z. + k_3z, = 0. (15) 

yoC ot; + ksze — (kK-3 + 62)2; = 0. (16) 


Aside from the tedious algebra these are readily solved to give equations (17) to 
(20) when only substrate 1 is present. 


BokikoaycaC 
A; 


Ai = A + Bayey1. (18) 


Fr(1) = By; = (17) 


(ky + hy) [hoB2 + h281 + 6182). (19) 
B= [heo(ky + Bo) + ki(k 2 + Be) |. (20) 


It is immediately apparent that if the stationary state is rapidly attained this 
should be a good approximation for the initial transport flux; F’7 gives the expected 
saturation behavior in ¢,;. 

The solution obtained when substrates 1 and 2 are present under the conditions 
C1 and cj constant, ¢.2 = Cy = 0, is given by: 


BokvaneaCy) kyks60 a hy6,(k 3 + bo) +6:k syo2e 2 
As 


rah) = 


Ar = D+ Earta + yeelF + Goren]. 
D = [kyd2 + 6,(k_3 + 5) JA. 
E = [hsb. + 8(k_3 + 6) |B. 
= (hy + k_g + 6) (hk koBo + kk 28) + k1B182) + 
k 36; [ko82 + Bi(k_» + Be) ] 
G = (ks + k_3 + 6))hoBo + k 36:(ko + ko + Be). 
The difference between the two fluxes is given by the equation 


Pe(1) — Fp(1) = Peheencaracalo 1 5 ki(ks + ks + 8)] X 
y(1) — Fe = OK_3 — hh bs 
X 7 AoA, 1 1 1 


[k2Boareer +k 1 [AoBe + Bi(k, + Be) |]. (27) 


Since all the rate constants and concentrations in equation (27) are positive, F’y(1) 
is greater than F'7(1) when 


5k ” aes hy (ky + k_3+ 6) > 0. (28) 
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Finally, the ratio of the two fluxes is given by the equation 


Fy(1) — [kykgs + kydi(k 3 + 52) + Sk 


: - (29) 
I r(1) ky, As 


Thus, all the carrier models can give F'y(1) > Fr(1) under the conditions imposed 
by equation (28). To distinguish between the different carrier models, we assume 


Fy > Fr and seek what effect metabolic inhibitors would have on Fy — Fr and 
F’,./F 7 when the two amino acids in the exchange diffusion experiment are different 
and when they are the same except for an isotopic label. It will be assumed that 
the rate constants for a labeled amino acid are the same as for the unlabeled. The 
effect of metabolic inhibitors will be to decrease the rate constant which is linked 
to cell metabolism in the various models. Thus, in the model of active transport 
of free carrier, the linkage to cell metabolism increases k, to give ky > k1. The 
effect of a metabolic inhibitor is to make /; decrease and approach /_, in value. 


Carrying out the above program we obtain the results given in Table 1. 


TABLE 1 
RELATIVE Errect OF METABOLIC INHIBITORS ON Fy AND Fr WHEN Fx > F7 


_ Direction of Change —- 
Meta 
bolically Fx Fx Fr FPx/Fr 
Conditions for linked Hetero Iso Hetero- Iso Hetero Iso- 
Model Fx > Fr constants Fy7 exchange exchange exchange exchange exchange exchange 
Active transport 6; > 2k s * os . 


of free carrier 


Active transport 
of carrier-sub 


strate complex 


Association reac 


tion model 


Dissociation re 


action model 
2k 


6 
v ¥ ¥ ¥ ¥ 


* Tf all diffusion parameters are much less than the association and dissociation rate constants so that equilibrium 


constants can be used 
t Decrease if Bik» ki (2k: + Bi) > O, otherwise increase 
t Not determinable unless relative rates of change of 82 and 4: are known 


First, it should be noted that equation (28) reduces to the same conditions for 
models 1, 3, and 4 in order for Fy(1) > Fr(1). Secondly, only the model of active 
transport. of free carrier predicts an increase in F'y/F7 and Fy — ir in all cases 
when a metabolic inhibitor acts. However, the directional changes in the second 
model could not be determined in hetero-exchange diffusion. Furthermore, the 
signs of d[Fy(1)/Fr(1)|/dky and d[Fx(1) — Fr(1)|/dkz, could not be determined 
in the isoexchange case for the model of active transport of carrier-substrate com- 
plex. However, these can be determined if one assumes that the reaction rates 
are much faster than the diffusion rates so that ‘equilibrium exists between carrier 
and substrate at each of the surfaces of the cell membrane. For this case, metabolic 
inhibitors decrease Fy(1) — F7(1) and Fx(1)/F7(1). By continuity, we see that 
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this must also be true for some range of parameter values where the equilibrium 
assumption is less tenable. For the association reaction model. metabolic inhibitors 
give a decrease in Fy(1)/Fr(1) and Fy(1) — Fr(1) in the isoexchange case. In 
hetero-exchange, F'y(1)/Fr(1) and Fy(1) — F7(1) decrease if Bik. — ky (2k2 + B,) > 0 
and increase. if Bik, — hy(2h. + 61) < 0. The former condition implies Fy(2) > 
F,(2) if the roles of the two amino acids are switched, whereas the latter gives 
Fy(2) < Fr(2). Again for the dissociation reaction model, metabolic inhibitors 
decrease F'y(1) — Fr(1) and Fx(1)/F7(1) in the isoexchange situation, but the 
corresponding changes in the hetero-exchange situation were not obvious from the 
expressions obtained. 

It should be emphasized that the results obtained are strictly applicable only to the 
stationary state which would be obtained if we could maintain the initial conditions 
which apply to measurements of initial uptake fluxes. Nonetheless, they demon- 
strate the theoretical possibility of obtaining Fy > Fr with carrier models. This 
so-called pseudostationary state approximation has been useful in chemical kinetics 
and has been basic to the derivation of the equations of enzyme kinetics. Since 
the amino acid active transport system is quite similar kinetically to enzyme sys- 
tems and since, in analogy with the latter, the amount of carrier present is probably 
far less than the amount of substrate, the stationary state fluxes may be fairly 
good approximations to the measured ‘‘initial’” uptake fluxes. Even if this is not 
true, we can expect that at least the directional changes due to inhibitors will be 
the same for initial uptake fluxes as for the corresponding stationary-state fluxes. 
For this reason, the effects of inhibitors have been considered primarily in terms 
of the directional changes in the fluxes. The induction period in amino acid 
active transport is certainly less than one minute. Thus, the pseudostationary 
state solutions should at the least be a good approximation to the kinetics after 
the process has gone on for a few minutes.'! However, with two amino acids 


present, even the pseudostationary-state equations are difficult to solve for the 


kinetics in the general case where we do not impose the restrictions that ¢,1, ¢,s. 
Ca, and cj» remain constant. 

The meager data available® on the effects of inhibitors on /'y and Fp fit in best 
with the model of active transport of free carrier. Obviously, more measurements 
on the effects of inhibitors on Fy/Fr7 and Fy — Fr, particularly in isoexchange 
experiments, are needed. 

The difference between the model of active transport of free carrier and the other 
models may be seen on a physically more intuitive basis by examination of Figures 
2 and 3. Let us call the closed cycle in Figure 2 the active transport loop and 
the outer closed cycle in Figure 3, involving both substrates, the exchange diffusion 
loop. Then, it is only in the model of active transport of free carrier that the 
metabolically linked step does not occur in the exchange diffusion loop. The flux 
F7(1) is obtained by operation of the active transport loop alone, whereas F’x(1) 
is obtained by the simultaneous operation of the active transport loop and the 
exchange diffusion loop. Since only that part of F’y(1) contributed by the active 
transport loop can be decreased by an inhibitor and even that can be decreased 
no more than can /’7(1), it is obvious why F'y(1) — Fr(1) and Fy(1)/F7(1) cannot 
decrease. In each of the other three models, two steps in the exchange diffusion 
loop are linked to metabolism. One of these steps is in the direction of the flux 





Vou. 47, 1961 BIOCHEMISTRY: J. A. JACQUEZ 163 


being determined, whereas the other opposes this flux. Thus, the effect of an 


inhibitor will depend on the relative changes in the two parameters involved. It 


would appear that the predominant effect is to decrease Py — Fr and Fx/F7, 
although this must be qualified somewhat for the model of active transport of 
carrier-substrate complex, because the findings are applicable only when carrier 
and substrate are near equilibrium. 

In conclusion, the use of the stationary-state approximation with carrier models 
gives results which check adequately with experimental data on the initial uptake 
fluxes in active transport and exchange diffusion. The data on the directional 
effects of inhibitors on these fluxes can also be fitted by the carrier models and 
the inhibitor studies may provide the data required to differentiate among the 
various types of carrier models. The mathematical solutions corresponding to 
the experimental steady-state have not been considered; the algebraic equations 
involved are so cumbersome it is difficult to see the forest for the trees. The 
results on competition also have not been considered. It is obvious from the 
structure of the models that the carrier models can give competition for uptake; 
it is not obvious that they can give an increase in initial uptake flux of one amino 
acid when another is present, as has been found with Ehrlich ascites cells. Un- 
fortunately, even with the pseudostationary-state approximation, the equations 
for this case remain difficult; it would probably be simpler to try a computer 
solution for this case. We conjecture that the amino acids which at low concen- 
trations give an increased initial uptake flux of another amino acid in a competition 
experiment will be those which give an increased initial uptake of that amino acid 
in an exchange diffusion experiment. 

* This work was supported by Grant T36 from the American Cancer Society. 
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THE ROLE OF PROTEIN SYNTHESIS IN EARLY ESTROGEN 
ACTION* 
By Greratp C. Muetier, Joun Gorskt,t anv Yosuro Arzawat 
MCARDLE MEMORIAL LABORATORY, UNIVERSITY OF WISCONSIN MEDICAL SCHOOL 
Communicated by Henry Lardy, December 7, 1960 

Previous studies from this laboratory on the mechanism of estrogen action have 
demonstrated that a single dose of estradiol causes a rapid acceleration of synthetic 
reactions which lead to the accumulation of phospholipids, ribonucleic acids, and 
proteins in the uterus of the ovariectomized rat.'~* In those metabolic pathways 
which have been investigated in cell-free systems (i.e., one carbon metabolism and 


protein synthesis), the early acceleration by estrogen was correlated with a simul- 


taneous increase in the activities of certain participating enzymes.?* These data 
suggested the possibility that the metabolic alterations in general may result from 
hormonally induced increases in enzyme levels. 

Since attempts to activate hypothetical proenzymes were uniformly negative, 
experiments were initiated to assess the role of protein synthesis in the estrogen-in- 
duced responses. The recent report of Yarmolinsky and de la Haba’ demonstrating 
that the antibiotic, puromycin, inhibited protein synthesis in rat liver enzyme prep- 
arations prompted its use in the present studies. This paper will demonstrate that 
with the proper in vivo administration of puromycin, protein synthesis can be blocked 
almost entirely in the rat uterus without inhibition of phospholipid or ribonucleic 
acid synthesis in this tissue. Under conditions in which protein synthesis was 
blocked by puromycin, the early acceleration of phospholipid and ribonucleic acid 
synthesis as well as the characteristic imbibition of water resulting from the in vivo 
action of estradiol failed to occur. In accordance with these findings, the require- 
ment for protein synthesis in the early estrogenic response is discussed. 


Methods.—Female rats of the same age and weighing approximately 180 gm were ovariecto- 
mized through the dorsal approach and maintained on a diet of commercial dog chow (Purina) for 
at least three weeks prior to experimentation. Only rats ovariectomized on the same day were 
used in any one experiment. 

Puromycin, 15 mg suspended in 0.5 ml of buffered saline (0.154 M NaCl in 0.04 M NaH2PO,- 
Na2HPO, buffer at pH 7.4) was injected intraperitoneally at 0, 1, 2, and 3 hours.§ Controls 
received only the buffer. Estradiol-178, 10 wg in 1.0 ml of buffered saline containing one per 
cent ethanol, or the buffered saline-ethanol control solution was injected via the tail vein at zero 
hour as indicated. Glycine-2-C™, 0.25 or 0.50 mg per dose, (specific activity 1.19 wcuries/umole ) 
dissolved in buffered saline was injected intraperitoneally as indicated at 0, 1, 2, and 3 hours in 
two trials and at 0 and 2 hours in one trial. In other experiments, inorganic orthophosphate-P®, 
600 or 750 ueuries dissolved in 0.5 ml of 0.154 M NaCl (5.9 x 10% weuries/umole), was injected 
intraperitoneally at zero time. 

The animals were killed at 4 hours; the uteri were removed, stripped of fat and connective tissue, 
and weighed on a Roller-Smith Torsion balance. The uteri were then placed in ice-cold five per 
cent trichloroacetic acid (TCA) and either pulverized from the frozen state in a stainless steel 
mortar cooled with liquid air or homogenized in an all glass homogenizer in cold five per cent TCA. 
Tissue residues were washed four times in ice cold five per cent TCA to remove the acid-soluble 
fraction. The lipid fraction was extracted at room temperature by successive washes of 80 per 
cent ethanol, 100 per cent ethanol, chloroform-ethanol (2:1), ether and ether.} 

In experiments in which glycine-2-C' was the radioactive precursor, the nucleic acid—protein 
residues were plated and counted in an internal gas flow counter. The results were expressed as 
counts per minute (C.P.M.) per mg of residue after correction for self absorption. Previous 
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studies have demonstrated that over 90 per cent of the radioactivity was due to the label in the 
protein of this residue.® 

The isolation of the nucleic acid purines was carried out according to the procedure of LePage.” 
For this purpose, the residues from two uteri were combined and extracted with four per cent per- 
chloric acid at 90° for 30 minutes. The mixed purine bases which were separated from the extract 
by adsorption and elution on ashort column of Dowex 50 (H* form) were purified to constant 
specifie activity by chromatography on paper first in a solvent system consisting of isopropanol, 
isobutyric acid, H.O, and 30 per cent NH,OH (60:25:14:1) and then in a system of N-butanol 
and H,O (86:14). Horizontal suspension of the paper during the chromatography gave the 
optimal separation of the bases. After elution of the bases from the final chromatograms, aliquots 
were plated and counted in an internal gas flow proportional counter and the quantity of adenine 
determined by the absorption at 260 mu. 

In experiments in which inorganic orthophosphate-P*? was the radioactive precursor, the lipid 
fraction remaining after evaporation of the combined organic solvent extracts was chromatographed 
on silicic acid-impregnated paper according to the method of Marinetti by using diisobutyl ketone, 
acetic acid, and HO (40:25:5) as the solvent system.!! The phosphotides were located by radio- 
autography. The ethanolamine phospholipid spot was eluted, aliquots were taken for counting 
in an end window counter, and the phosphate was determined by the method of Chen.'? Results 
are expressed as C.P.M. per ug of P. 

The labeling of the ribonucleic acid fraction was determined on the ribonucleotides isolated 
from the alkaline degradation of the RNA. For this purpose, the pooled protein-nucleic acid 
residues from two uteri were hydrolyzed in 0.2 N NaOH at 37° for 14 hours. After precipitation 
of the acid-insoluble residue with cold ten per cent perchloric acid, the nucleotides were adsorbed 
on charcoal (nuchar C) and then eluted with 2.5 per cent pyridine in 50 per cent ethanol as described 
by Tsuboi.!* The chromatographic isolation of the AMP and CMP was achieved by chromatog- 
raphy on paper at 5°C using isobutyric acid and 0.4 M NH,OH (75:25) as the solvent system. 
The nucleotide areas were eluted and aliquots were counted in an end window counter and 
assayed for AMP and CMP by ultraviolet absorption. Results were expressed as C.P.M. per 
umole of nucleotide. 


Results.—The effect of puromycin on the incorporation of glycine-2-C'* into uterine 
protein: In order to establish the in vivo effectiveness of puromycin as a blocking 
agent for protein synthesis, control and puromycin-treated rats were injected 
with glycine-2-C', and the amount of isotope incorporated into uterine protein 
during a four-hour period was determined. In preliminary experiments with 
varying doses and schedules, it was found that the intraperitoneal injection of 15 
mg of puromycin each hour gave approximately a 90 per cent inhibition of the in- 
corporation of glycine-2-C' into uterine protein of the control rats over the test 
period (Fig. 1). This dose of puromycin was similarly effective in blocking protein 
synthesis in the estrogen-treated rats. It was also observed in the course of these 
experiments that the synthesis of uterine proteins had fairly well resumed by two 
hours after a single 15-mg dose, thus necessitating the hourly injections of puro- 
mycin. The experiments were confined to four-hour periods so as to provide 
adequate time for highly significant estrogen responses and yet not so long as to 
incur secondary toxicities resulting from the blockage of protein synthesis with 
puromycin. 

In another report,'* it will be shown that this dose of puromycin is equally 
effective in inhibiting the incorporation of glycine-2-C' into proteins of the liver 
but was less effective in the heart, thymus, and kidney (64, 51, and 60 per cent 
inhibition, respectively). 

The effect of puromycin on the incorporation of glycine-2-C™ into nucleic acid 


purines: In contrast to its action on protein synthesis, puromycin had little or no 
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effect on the incorporation of glycine-2-C'™ into the adenine of the mixed nucleic 
acids of the control uteri (Fig. 2). On the other hand, the characteristic accelera- 
tion of nucleic acid purine labeling during the first four hours of estradiol action was 
completely prevented by the simultaneous administration of puromycin; it is of 
interest that the incorporation rate actually remained at the control level throughout 
the four-hour period. 
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Fic. 1.—Glycine-2-C" incorpora- Fic. 2.—Glycine-2-C'* incorpora- 

tion into protein—nucleic acid tion into nucleic acid adenine. 
residue of rat uteri. Each bar is Same experiment as in Figure 1. 
the average of 3 (—puromycin) or Protein-nucleic acid residues were hy- 
4 (+puromycin) rat uteri. 15 mg drolyzed in 4 per cent PCA at 90°C. 
of puromycin (+) or a control Adenine was isolated by a combina- 
solution (—) were injected intra- tion of column and paper chroma- 
peritoneally at 0, 1, 2, and 3 hr. tography as described under 
0.5 mg of glycine-2-C! (6.2 uC) methods. Specific activity was 
were injected intraperitoneally at determined by counting C" activity 
0,1,2,and3hr. 10 ug of estradiol- and measuring adenine by _ its 
178 (+) or a control solution (—) ultraviolet absorption. Results in 
were injected via the tail vein at this figure are expressed as per cent 
0 hr. Rats were killed at 4 hr and of the control value which was 16 
the uteri removed. The protein C.P.M./umole adenine. 
nucleic acid residue was obtained by 
procedures described under methods. 
The residue was plated and counted 
and the results calculated as C.P.M. 
per mg protein residue after cor- 
rection for self absorption. The 
data in this figure are expressed_as 
per cent of the control value which 
was 428 C.P.M./mg of protein 
residue. Brackets indicate range. 


The influence of puromycin and estradiol on the incorporation of radioactive phos- 
phorus into ribonucleic acids of the rat uterus: In experiments in which inorganic 
orthophosphate-P*? was used as the radioactive precursor, it was demonstrated that 
puromycin treatment resulted in an actual increase of the labeling of RNA (Fig. 3) 
of the control tissues. However, as in the case of glycine-2-C'" the acceleration of 
RNA synthesis by estradiol was largely prevented by the action of puromycin, and 
the labeling was essentially the same as seen in the puremycin-treated control rats. 
The reason for the increased labeling in the uterine RNA of the control rats which 
were treated with puromycin alone is not apparent at this time. 
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Alterations of phospholipid synthesis by estradiol and puromycin treatment in vo: 


In previous experiments with surviving uterine segments from rats pretreated with 


estradiol, it has been demonstrated that the synthesis of phospholipids is an un- 
usually sensitive indicator of early estrogenic action.! This observation has been 
reconfirmed and extended in the present in vivo study. In Figure 4, the results of 
an experiment are shown in which the incorporation of inorganic orthophosphate-P* 
into the ethanolamine phosphatide fraction has been measured in vivo over the 
first four hours following the administration of estradiol. It is evident that the 
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Fic. 3.—The incorporation of Fic. 4.—The incorporation of 
orthophosphate-P*? into uterine orthophosphate-P*? into phospho- 
RNA. Each bar is the average of ethanolamine. Same experiment as 
four rat uteri. 15 mg of puromycin Figure 3. Lipid fraction was 
or a control solution was injected chromatographed on. silicic acid 
intraperitoneally at 0, 1, 2, and 3 impregnated paper as_ described 
hr. 622 ywC-inorganic orthophos- under Methods. fesults are ex- 
phate-P** was injected intraperi- pressed as per cent of the control 
toneally at O hr. 10 yg estradiol- value which was 151 C.P.M./ug of 
178 or a control solution was ethanolamine phosphorus 
injected via the tail vein at O hr 
tats were killed and uteri removed 
at 4 hr. RNA was degraded with 
alkali, and the 2’, 3’ AMP isolated 
by paper chromatography as de- 
scribed under methods. tesults 
are expressed as per cent of the 
control value which was 2360 


C.P.M./umole of AMP 


single dose of estradiol increases the incorporation of radioactive phosphorus into 
this phospholipid fraction over 800 per cent. 

As in the case of the labeling of the RNA, the administration of puromycin 
to the control rats also stimulated the labeling of the ethanolamine phosphatide 
fraction significantly; as yet, this stimulation is not understood. However, the 
administration of puromycin to rats treated simultaneously with estradiol pre- 
vented almost entirely the estrogen acceleration of the phospholipid labeling. In- 
stead, the amount of P*? incorporated into the ethanolamine phosphatide fractions 
was of the same order as obtained in the control rats treated with puromycin alone. 

Effect of estradiol and puromycin on the uterine wet weight: A characteristic 
action of estrogens widely used for their assay is the hormonally induced imbibition 
of water in the uterus of the overiectomized rat. In Figure 5, the interference with 
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this phenomenon by puromycin treatment is shown. While puromycin treatment 
had little or no effect on the wet weight of the control uteri, the estradiol-induced 
inhibition of water was almost completely blocked. This experiment is repre- 
sentative of seven separate experiments. 

Discussion.—The data reported above demonstrate that puromycin at the proper 
dosage is a highly effective agent for blocking protein synthesis zn vivo in the uterus 
of both control and estrogen-treated rats. Under the same conditions, the synthesis 
of nucleic acid purines from radioactive glycine and the incorporation of inorganic 
orthophosphate-P*? into the uterine RNA and phospholipids of puromycin-treated 
rats continued at the control level or was stimulated somewhat. Thus, while 
protein synthesis was interrupted, a number of other basic synthetic processes 
continued at near the norma! level. However, in all cases, the early and character- 
istic accelerating action of estradiol on phospho- 
lipid and RNA synthesis was prevented by the 
simultaneous action of puromycin in vivo. Ac- 
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Bie a re estrogen action is protein synthesis. The data re- 
begin? _ ported in this paper are in accord with the selective 

Fic. 5.—The effect of puromycin 
on the imbibition of water in 
uteri of estrogen-treated rats. Same point to the usefulness of this agent in the investi- 
experiment as in Figure 1. Each 
bar is the average wet weight of four 
uteri 4 hr after initiation of the in Che observation that the characteristic imbibi- 
vivo treatments. Results are ex- 
pressed as per cent of the average 
control value which was 65 mg per be blocked by puromycin treatment in vivo argues 
uterus, 


inhibition of protein synthesis by puromycin and 
gation of other metabolic problems. 
tion of water into estrogen-treated uteri can also 


decisively against any direct alteration of cellular 
permeability by the hormone. The results, in- 
stead, suggest that even this aspect of the estrogen response is dependent on a 
puromycin-sensitive reaction which appears to involve protein synthesis. The 
conclusion that estrogens do not act by altering permeability had been arrived at 
previously by Halkerston ef al.!> through direct studies on uterine permeability. 

In these studies, it would be of interest if some estrogen action could be demon- 
strated which was not blocked by puromycin. So far, evidence of such a process 
has not been forthcoming. Similarly, all attempts to demonstrate an estrogen- 
sensitive uterine transhydrogenase!® and hence its participation in the estrogenic 
response in the rat uterus have been negative. 


In a previous discussion of the mechanism of action of estrogens, it has been 
pointed out that an understanding of the problem awaits both the identification of 
the primary molecular interaction of the hormone with the cell receptor and the 
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elucidation of the system whereby the effects of this interaction are amplified.’ 
The data presented here suggest that the process of protein synthesis is the cellular 


amplifier and that the estrogen response is expressed through the production of rate- 


limiting enzymes. Accordingly, attention is now directed to a study of the manner 
in which estrogens accelerate protein synthesis in the uterus. Whether estrogens 
act primarily to activate or inhibit some existing enzyme which in turn contributes 
to the improvement of the intracellular environment for protein synthesis in general 
or whether it induces the de novo synthesis of such an enzyme remains to be ascer- 
tained. 

Summary.— The administration of puromycin to ovariectomized rats blocked 
uterine protein synthesis in vivo without inhibition of phospholipid or ribonucleic 
acid synthesis. Under the same conditions, puromycin treatment also prohibited 
the early acceleration of phospholipid and ribonucleic acid synthesis by estradiol 
and prevented the hormone-induced imbibition of water in the uterus. The essen- 


tiality of protein synthesis for early estrogen action is discussed. 
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A KINETIN-LIKE COMPOUND IN MAIZE* 
By Carbos O. MILLER 
DEPARTMENT OF BOTANY, INDIANA UNIVERSITY 
Communicated by Ralph E. Cleland, December 27, 1960 


Since the discovery and isolation of kinetin! and its characterization as 6-fur- 
furylaminopurine,? many analogues active in cell division have been synthesized 
and numerous effects of kinetin have been described. The question as to whether 
kinetin or a closely related compound occurs naturally in plants, however, has 
remained open. Last spring, in the hope of obtaining an answer to this question, 
we began working with materials from Zea mays kernels still in the milk stage. As 
many others have demonstrated,* the milky endosperm of maize contains sub- 
stances which promote plant cell growth. Our progress in identifying one such 
substance is summarized herein. 

Assay.—We observed that water or alcohol extracts of the immature kernels 
substituted very effectively for kinetin in the growth of soybean (Glycine max) callus 
tissue which had been derived from cotyledons of sterile plants. The variety Acme 
was used, although several other varieties worked quite well. The stock tissues 
were cultured on a medium containing (mg/l): KH.PO,, 300; KNOs, 1,000; 
NHyNOs3, 1,000; Ca(NO3)2-4H2O, 500; MgSO,, 35; KCl, 65; KI, 0.8; MnSQ,, 
14; ZnSO,, 1.5; H3BOs3, 1.6; glycine, 2; thiamin-HCi, 0.1; nicotinic acid, 0.5; 
pyridoxine: HCl, 0.1; indole-3-acetic acid, 5; kinetin, 0.5; sucrose, 30,000; Bacto- 
agar, 10,000. The medium was adjusted to pH 5.8. To test a preparation for 
kinetin-like activity, it was added to this medium and the kinetin was omitted. 
Ordinarily each of four 125-ml Erlenmeyer flasks containing 50 ml of hardened 
medium was planted with three pieces of soybean tissue which had been growing 
on the medium with kinetin. The tissues were grown at 27°C and constantly ex- 
posed to fluorescent lighting at about 40 ft.-c. The pieces were weighed after 3 
or 4 weeks. 

Purification of Factor —To prepare extracts, ears of maize with grains in the 
milky stage were scraped with a simple creaming device. This removed the endo- 
sperm tissue along with a large percentage of embryos from the cob. The watery 
mixture is hereafter referred to as creamed maize. In an early experiment, a 65 


per cent ethanol extract of creamed maize was prepared at room temperature. The 


alcohol was then removed under vacuum at 50°C. The resulting product, called 
alcohol extract in Table 1, was quite effective in replacing kinetin. This extract 
showed a division of activity when run through a Dowex 50 W-X8 (H+, 50-100 
mesh) column having a bed volume of 600 ml. Some active material came through 
the column and was part of the column effluent mentioned in Table 1. Other active 
material was recovered by eluting with 6 N NH,OH after the resin had been washed 
with a large volume of distilled water and the ammonia was removed under vacuum 
at 60°C. For each gram of fresh creamed maize, the dry weight in this eluate 
amounted to only 1.9 mg. This preparation was tested for promotion of growth 
of soybean callus and it was very active (Table 1). Microscopic examination re- 
vealed that cells of the treated tissues were as small as or smaller than those of 
control tissues, thus indicating that the factor was promoting cell division. Al- 


170 
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though an unequivocal interpretation of the data presented is not possible, it seems 
that the column treatment did not produce active material since the crude alcohol 
extract showed more activity than either the effluent or eluate when each was 
tested at a level equivalent to 10 gm of creamed maize per 100 ml media. Other 


purifications, in which the temperature was kept below 35°C and in which no NH,OH 


or ethanol was used, all yielded about the same levels of activity as obtained with 
the above procedure. lIurthermore, the removal of ethanol before the resin 
treatment was found to have no effect. The likelihood that the active material 
was produced by our extraction and purification procedures therefore seems small. 

The active material retained by the column behaved much as kinetin; that is, 
it was held by the resin under slightly acidic conditions, was easily eluted by base, 
but was eluted only slowly by acid at rather high concentrations. Since we were 
looking specifically for a kinetin-like material, we attempted to learn the nature of 
the factor held by the resin. 

This past summer, we prepared 84 kg of creamed maize (variety Golden Bantam) 
and extracted it with 65 per cent ethanol. This extract, still containing the ethanol, 


TABLE 1 


PROMOTION OF SOYBEAN CaLLus TIssuE GROWTH BY EXTRACTS FROM MAIZE 
KERNELS 


Fresh weig 
per piece 

Additive Concentration mg 
14 
Alcohol extract 107 913 
Column effluent 10t O31 
Eluate 1t 748 
2 5T ,494 
10F Ole 
20t , 690 
Kinetin 0.5 mg/l 984 


* Initial fresh weight was 24 mg. 28 days of growth 

Tt Expressed as the equivalent number of grams of creamed maize per 100 ml of medium 
was run through the column already described in batches equivalent to 6 kg of 
maize. The column was washed thoroughly with double-distilled water, eluted 
with 1500 ml of 6. N NH,OH and the ammonia was removed under vacuum at 60°C. 
The acidity of the preparation was next adjusted to pH 1.5 with HCl. The acidified 
solution was then added in batches equivalent to 10 kg creamed maize to a smaller 
Dowex 50 (H~*) column having a bed volume of 25 ml. The column was washed 
with 100 ml water and then eluted (0.8 ml/min) with 600 ml 1.5 NV HCl followed 
by 500 ml 4 V HCl. Fractions of 7 ml were collected. Shortly after the start of 
the 4 N HCI (tube 86), a peak of material absorbing ultraviolet light in the 260 
270 mu region was detected; this material consisted of several compounds. The sub- 
stance active in the soybean growth test was collected in fractions shortly after 
this peak of absorbing materials. In most runs, tubes 93 to 130 contained the 
factor, and their contents were pooled. 

These combined fractions were then diluted with an equal volume of water and 
again run through a small Dowex 50 (H+) column. After the column had been 
washed with water, the factor was eluted with 6 N NH,OH, and the ammonia was 
removed under vacuum. The remaining material was then concentrated and 
applied in stripes onto sheets of Whatman #1 chromatography paper. 
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did the more gentle treatment. The derivative exhibited an Re of 0.66 in the 3 
t-butanol:1 H.O:1 NH,OH solvent system and 0.50 in the 4 n-butanol:1 
H.0:1 glacial acetic acid system. Bringing a 6 N NH,OH solution of the factor 
to a boil did not affect the quenching compound. Neither was any change achieved 
by heating a 0.1 N NaOH solution of the factor at 100°C for 5.5 hr. 

Both the quenching compound and the derivative obtained by the 1 NV HCl treat- 
ment apparently formed complexes with mercury since both yielded detectable 
spots when chromatograms were sprayed with mercuric acetate followed by di- 


phenylearbazone.® The factor also complexed with silver and was precipitated 


at pH 5.5 but apparently not at pH 1. Neither compound produced any color 
when sprayed with a ninhydrin solution or with the modified Dische reagent’ which 
is effective on kinetin. 

Conclusions.—Our results demonstrate that young maize kernels contain more 
than one substance which can cause the same type of growth effect as kinetin. 
The chemical relationship of these factors is presently unknown. However, if the 
factor we have purified is a purine somewhat similar to kinetin, it might occur in 
several nucleic acid derivatives and thus account for the various active substances. 


TABLE 3 
ABSORPTION CHARACTERISTICS OF THE QUENCHING MATERIAL ASSOCIATED WITH ACTIVITY 
95% Ethanol 0.1 N HCI 1 N NaOH 

Absorption max (my 269 273 ; , 2744 
Absorption min (my 233 : 24. 
O.D. max* 0.71 0 
O.D. min* 0.19 es 
O.D. 260 m : 

a 1.61 
O.D. 250 my 
O.D. 280 mu 
O.D. 260 mu 


0.80 


* Fora solution containing 21 micrograms of purest preparatio 
t Aslight shoulder in the 280-285 my region was evident 
The factor which we have greatly purified is like kinetin jn that (a) both are 
held on Dowex 50 under the same conditions and are removed from the resin by 
essentially similar treatments, (b) both strongly absorb ultraviolet light, showing 
peaks in the 270 my region, and exhibit similar shifts when the pH is changed, (c) 
both stimulate cell division in the soybean test at low concentrations, and (d@) both 
apparently form complexes with mercuric and silver ions, as do many purine 
derivatives. The two compounds apparently are different in that (a) the factor 
is changed by a fairly mild acid treatment, whereas kinetin is not affected, (b) 
the factor produces no color with the Dische reagent, whereas kinetin gives a pink 
color thus indicating that the factor does not contain deoxyribose or a sugar 
derivative such as the furfuryl group found in kinetin, (c) the two show different 
Ry values and can be separated by chromatography, and (d@) the factor does not 
dissolve in ether as does kinetin. Obviously, the two compounds are different 
chemically. They nevertheless may be closely related since our data support the 
idea that the factor is a substituted purine. Such a relationship can be verified 


only by further analysis of the maize factor. 


I wish to thank Patricia Hall and John Miller who have ably and cheerfully assisted me in this 


work. 
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STABILIZATION OF NH IN HYDROCARBON MATRICES 
AND ITS RELATION TO COMETARY PHENOMENA 


By Jay A. GLASEL 
SCHOOL OF SCIENCE AND ENGINEERING, UNIVERSITY OF CALIFORNIA AT SAN DIEGO 


Communicated by Harold C. Urey, December 16, 1960 


The work described in this paper is the result of some initial studies of the low- 
temperature properties of trapped free radicals resulting from the breakup of some 
of the simple hydrogen compounds of nitrogen, oxygen, and carbon. The purpose 


in beginning these studies was to explore the possible significance of free radicals 
in astronomical matter, such as comets and interstellar grains. Their importance 
to these problems was first pointed out by Urey! and Donn and Urey.” 


Experimental.—The low-temperature optical dewar used in these experiments was a hybrid 
of several others.* It featured an absorption cell which could be rotated inside the dewar without 
loss of vacuum in order to (1) face the sample beam of a double-beam recording spectrophotometer, 
(2) face a gas discharge source of free radicals, (3) face a photolysing beam of ultraviolet light. 

The outer, room-temperature, optical cell contained windows of calcium fluoride. This mate- 
rial, while somewhat brittle, is easily ground and polished by using the chelating action of EDTA 
on the calcium ion. Sapphire was used as the optical material for the inner absorption cell. 
This cell consisted of a single window, 1 mm thick, upon which the molecules to be studied were 
condensed. Sapphire was chosen because of its good thermal conductivity at low temperature 
and, as will be discussed below, because of its high index of refraction. The window was set in a 
small copper block which was attached to a boiling refrigerant bath. Further details of the cell 
are not pertinent to this discussion, other than the fact that a copper-constantan thermocouple 
attached directly to the window was used to measure temperatures when liquid hydrogen was used 
as the coolant, and a carbon resistor thermometer,' set in the block, was used for the same purpose 
at liquid helium temperatures. Thermal contact between the window and the metal block was 
achieved by a paste of aluminum powder and stopcock grease. 

The products reported here were produced by a 175-input-watt RF transmitter coupled tightly 
to a discharge tube of 10-mm Pyrex tubing through which the parent gases traveled. This tube 
was inserted into the outer cell of the apparatus so that the inner sapphire window could be ro- 
tated to face it. The end of the discharge tube facing the sapphire window was formed into a 1-mm 
hole through which the discharged gases passed and traveled 3 cm to the window. Gases not 
condensed were immediately pumped away by a liquid nitrogen trap and mercury diffusion pump 
system. The residual pressure in the cell while these two competing processes were going on was 
measured by anion gauge. The flow rate of the parent gases into the discharge tube was regulated 
by a needle valve. 

The spectrophotometric apparatus consisted of a Beckman DK-1 recording prism spectro- 
photometer. The instrument has a range of 1850 A to 3.5 microns. The resolution in the UV 
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region of present interest Is about 3 A The spectrophotometer was calibrated in the visible and 
UV with a mercury are lamp 

The CH, used was Phillips research grade which was analyzed by « mass spectrometer and found 
to be 99.6297 CHy, 0.35% No, and 0.08°7 CsH¢s. This was used directly to make up CHy-N» and 
CH,-No-A mixtures with no further addition of No. The effect of the C,H, on the experiments 
was assumed to be negligible. The argon was Linde research grade and was used without further 
purification. The NH; was Matheson anhydrous grade, which was purified in the case of the 
20°K experiments by condensing it in a flask cooled with liquid nitrogen and pumping on it for 
several hours and, in the case of the 6°K experiments, by collecting it from a sodium metal solu- 
tion. Mass spectrometer analyses were unable to detect NH; for technical reasons, but samples 
of NH;-A run through the mass spectrometer revealed no traces of other gases except No in small 
amounts. No oxygen could be found in any of the mass spectrometer analyses. H.O could be 
detected by its near-infrared spectrum in very small amounts (less than 0.10 ) but was not found 
in any of the experiments discussed here. Nitrogen and hydrogen gases used in one of the mixtures 
were Linde research grade used without further purification. All bands reported were not present 
in blank runs with the undischarged gases 

The fact that sapphire has a high index of refraction makes it possible to utilize the reflection 
interference patterns resulting from the films of condensed material in order to calculate their 
thickness. The patterns are observed in absorption with the spectrophotometer as it plots ab- 
sorption versus wavelength for a given film in the visible re gion of the spectra. In the thicknesses 
used, usually about 20 peaks are measurable. If Av is the difference in wave-number between 


successive absorption maxima, then, 
Av 


In order to calculate the thickness, d, of the film, a somewhat arbitrary estimate must be made of 
the refractive index, n, of the film. In this paper, the films were all assumed to have an index of 
refraction 1.4. The oscillator strength of NH has been found experimentally to be 0.008,° and 
this may then be combined with the film thickness to give a value for the concentration of the 
radicals to +20°;. The error arises mainly from the refractive index assumptions. Typical 
values are shown in Table 3. 

Warmup samples of the materials frozen on the window were obtained by sealing off the vacuum 
pumps and allowing the dewar to warm up to room temperature. The gases liberated were 
allowed to enter an evacuated sample tube. The sample was then analyzed with a commercial 
mass spectrometer. Thus, an analysis represents only materials volatile at room temperature and 
includes gases that were given off by the walls of the metal dewar at room temperature. Materials 
that were absorbed by the walls were not represented in the analyses. Blanks of pure argon run 
through the equivalent process of sample-taking showed that the only gas given off by the walls 
was nitrogen. This accounts for the increase in nitrogen from the original sample to the warmup 
products shown in Table 2. The quantitative figures given in this should only be taken as repre- 


senting approximately the proportions of materials found in the sample 


TABLE 1 
PEAKS IN SPECTRA OF DISCHARGED MIXTURES 


Mixture , Peak » Assignment 

NH:-A 29,637 em NH, (0,0 
10,683 4 

29 949 NH, (0.0) 
39,761 ? 

10,000 


30,231 


CH,-N> 


NH, 0.0) 


CH,-N2-A 


y 


Le RS Oe Oe ene) 


N.-H -A 
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TABLE 2 


{EPRESENTATIVE MASS SPECTROMETER ANALYSES OF INITIAL SAMPLES AND 
Warmup Propucts 
Warm products 


Contents of sample in order o 
Mixture before discharge concentration 


CH,-N> CH, 99 6°; Methane 12.5 
N2 0.3% Acetylene l 

thane 0.08% Kthylene 
Nitrogen 
Ethane 
1,3 Butadiene 
Hydrogen 
Propylene 
Propane 
Butenes 
Methane 
Acetylene OL: 
Nitrogen 0.! 
Kthylene 0 
Argon 4 
Not analyzed 


—ho ke Ore 
bot ves 


A 99 9°; 
(NH; not detected) ratio 
of NH; to argon was 1:30 


Experiments at 20°K: The initial experiments in this series were done with the sapphire 
window cooled to 20°K with liquid hydrogen. The mixtures passed through the discharge tube 
were NH;-A, CH.-No-A, and CHy-N.. The ratios of the various molecules before discharge are 
shown in Table 2. Figures 1 and 2 and Tables 1 and 2 give the results of these experiments. The 
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Fic. 1.—Spectra of NH resulting from RF discharge through (1) CHy-N2-A and (2) 
CHy-N». Results from discharge in NH;-A are similar but not shown. Left side of 
figure shows result, (1’), of a diffusion experiment on film (1). Spectrophotometer on 
I/Iy scale with electronic magnification. Each unit of the vertical scale is 1% T. The 
scale runs from 90% T on the bottom to 100°, T on the top. Spectra arbitrarily displayed 
for illustrative purposes. 
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Fig. 2.—Spectra of unknown species resulting from RF discharge through mixtures given 
in Figure 1. Same films as Figure 1. Spectrophotometer on I/I, scale with no magnifica- 
tion. Bottom is 0°, T and top is 100°, si Spectra arbitrarily displaced for illustrative 
purposes 


figures show that the peak observed by previous workers in NH;-A’ and assigned to the (0,0) band 
of NH occurs in all the discharged mixtures. There are two new peaks which differ slightly in 
wavelength maximum in different matrices (Fig. 2). The films were all deposited over a time of 
about 3 hours, at a pressure in the dewar cell of 10~> mm. 

It should be noted that the warmup products of the CH,-N» mixture indicate that hydrogen gas 
was trapped in the condensing film. This was also observed in an experiment using a mixture of 
hydrogen and nitrogen 

Diffusion experiments, which are not all illustrated, were made on the films of the discharged 
mixtures. These involved allowing a film to warm to 40°K for a few seconds and then recooling 
it to 20°K. The spectra taken before and after the warming were then compared. In all cases, 
the peaks decreased noticeably in total absorption. This included the NH peaks and the peaks 
around 2450 A 

Experiment at 6°K: Using the same technique as above, mixtures of NH;-A and H.-N» were 
discharged and condensed on the window which was now cooled with liquid helium. The results 
are shown in Figure 3. Several features stand out from the previous experiments. In none of 
the experiments at 6°K were the peaks at 2450 A evident. As shown in Table 3, the apparent 
concentration of the radical NH in the N.-H» mixture was higher by a factor of ten than any of the 
other mixtures. The NH band was much broader than in the 20°K experiments and was par 
ticularly broad in the mixture of N»-H The instrumental slit width was much narrower than 
the width of these bands. No mass spectrometry was attempted with these particular mixtures 


TABLE 3 
APPROXIMATE EXPERIMENTAL VALUES OF NH CoNCENTRATION IN VARIOUS 
FILMS 
Mixture discharged NH Mole 
and condensation temperature concentration in condensed film 
CH,y-N2(20°K 0.011 
NH.-A(6°K 0.015 
N.-H2(6° RK) 0.18 
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Fic. 3.—Spectra of NH at 6°K resulting from discharges through the mixtures indicated. 
Spectrophotometer on I/Ip scale with no magnification. Bottom is 0°) T and top is 100°) T 
Note therefore that these bands are much stronger than those of Figure 1. Spectra arbi- 
trarily displaced for illustrative purposes. 

Discussion of Observations.—The fact that a chemically inert matrix of a rare gas 
is not needed to stabilize NH is illustrated by the CHy-N. experiments. Moreover, 
the presence of the larger organic unsaturates in the warmup products indicates 
that the radical can coexist with these in the film. In the case of the major one of 
these, acetylene, its spectrum was obtained in the near infrared at the same time 
the spectrum in the UV region showed the NH. This indicates that the organics 
are not formed in the warmup process but are formed in the discharge. The re- 
sults also show that large amounts of nitrogen in the parent gases are not necessary 
to the production of observable amounts of NH. No energy releases were ob- 
served in the films either upon standing at 20°K or in warmups. 

On the basis of these experiments, a definite assignment cannot be given to the 
peaks at 2450 A. A mixture of Hy-N» discharged and then condensed at 20°K 
indicated that the peaks at 2450 A appeared but that the NH peak did not. The 
peaks due to the NH, radical reported by Robinson and McCarty*® were not ob- 
served in these experiments, so that the peaks are not associated with this radical. 
It should be noted that the splitting of 240 em~! between the two peaks in the 
CH,-Ne experiments is of the same order of magnitude as may be expected for a 
hydrogen bending vibration. The cleanliness of the discharge tube had no effect 
on these peaks. The possibility that they were due to the y band of NO was in- 
vestigated by depositing a mixture of NO in argon (matrix/gas = 150) on the cold 
window and observing the spectra from 6° K to 30°K. A typical spectrum is shown 
in Figure 4. Since the spectra of the NO are not similar to the unknown bands, 
and since the splitting of the latter does not correspond to the expected 2IIs,, — *II,,, 
splitting (121 em~), it is improbable that they are the NO band. 
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The results of the concentration determinations show that the lowering of the 
condensation temperature alone does not bring about greater concentrations of 
stabilized NH. The important factor in these experiments seems to be a larger 
concentration of available N and H in the discharge mixture. Thus, the No-H» 
mixture probably represents the optimum for production of NH while the matrix 
makes little difference. It must be concluded that NH is by far the most stable free 
radical yet studied. Carbon monosulfide, which may not be a radical, is the only 
species thus far which exceeds the stability of NH. From these studies, it is clear 
that only small concentrations of NH may be stabilized in condensed films at low 
temperatures. 

Implications concerning Cometary Phenomena.—In relation to the original prob- 
lem, therefore, it is concluded that small concentrations of NH may be built up in 
frozen solids of thermodynamically unstable molecules in an environment similar 
to that existing at the surfaces of comets. The presence of large proportions of 
acetylene, ethylene, ethane, propane, and propylene in the condensed solids sug- 
gests that these compounds are responsible for the C2 and C; bands observed in the 
spectra of cometary comae. Since the gas density in a comet’s coma is so low as 
to preclude chemical reactions, the C. and C; observed must have their source in 
evaporation products from the vicinity of the nucleus. Considering the tempera- 
ture at the nucleus (<150°K) when the spectra are observed, the most likely source 
is the two- and three-carbon hydrocarbons which are evaporated from the surface 
and stripped of their hydrogen by photolysis. If this interpretation is true, it 
would be expected that the spectra of comet comae should be very rich in the bands 
of Hy and H+ since hydrogen must be lost from the surfaces at 150°K following 
the reactions which produce acetylene, ete. Since these bands are in the far UV 


they have not yet been observed. 
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The first part of the Donn-Urey hypothesis concerning the nature of cometary 
outbursts has been demonstrated. That is, it is possible to condense small quanti- 
ties of very reactive species along with molecules with very high free energy at low 
temperatures. The addition of H.O to the reaction mixture gives some amount of 
condensable H2O,. There are two possible subsequent courses for the peroxide in 
the solid. It is possible that it reacts slowly and continuously with the hydro- 
carbons yielding compounds with low free energy. However, it is also possible 
that at 150°K there is a certain amount of segregation of microcrystals going on 
which would keep the crystals of unsaturated hydrocarbons and those of the H,O. 
separated. Such a mixture would be able to store a very large amount of energy.® 
It now remains to investigate the behavior of these mixtures when they are bom- 
barded with beams of positive and negative particles of astronomically possible 
energy and density. If the hypothesis is correct, bombardment under the appro- 
priate conditions should lead to violent reactions of the sort observed at the heads 
of comets. 
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CORRELATED CHANGES IN ENVIRONMENT AND LETHAL 
FREQUENCY IN A NATURAL POPULATION OF DROSOPHILA 
MELANOGASTER 
By H. T. Banp!: 2 anp P. T. Ives 
DEPARTMENT OF BIOLOGY, AMHERST COLLEGE 


Communicated by I. Michael Lerner, December 5, 1960 


The frequency of recessive lethal chromosomes in American populations of Dro- 


sophila melanogaster has been shown to increase as one moves southward.*: 4 How- 
ever, within populations this frequency is not as high as it could be on the basis of 
assumed mutation rate versus measured amount of allelism.* Furthermore, if 


mutators sometimes cause even higher than assumed mutation rates,> the observed 
frequency of lethals within populations should be even greater. On the contrary, it 
has been found that in the 1947-1952 period, the frequency of recessive lethals 
dropped sharply in both northern and southern areas.4- To account for the sub- 
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maximum frequency of recessive lethals, the hypothesis was presented that adverse 
environmental conditions in nature might initiate selection against lethal heterozy- 
gotes, and one instance of such selection by temperature was cited.’ 

The present paper will report a significant decrease in lethal frequency in the 
South Amherst area within one month and will present evidence that the adverse 
environmental condition in part responsible for this and previous declines is daily 
temperature range in the natural environment during the week prior to the collec- 
tion of flies. 

Materials and Methods.—In 1958, male flies were collected from the South 
Amherst population in September and again in October.’ Second and third chromo- 
somes were extracted by the Cy/Bl and T M 1/D techniques, respectively, but for 
present purposes, only data from the second chromosomes will be presented. Cul- 
tures were raised at 25° + 0.5°C, and at leagt 100 flies were counted for each line in 
the test generation. Lines giving less than 17 per cent non-Cy flies were scored as 
lethals in the present study; only those giving less than 2 per cent non-Cy flies were 
used in the allelism tests.’ 

Temperature readings used in the present study were begun by Ives at the col- 
lection site in 1945. After 1948, they were continued by him at this department’s 
weather station, four miles from the collection site. The small difference in tem- 


perature readings at the two locations proved to be of no consequence in our subse- 
quent calculations and did not affect daily temperature range. Hence, the tem- 
peratures shown in Table 2 refer to the actual temperatures observed. 


Although it is not known how long adults usually survive in nature, in our labora- 
tory the life span of 200 flies from an inbred line averaged 70 days and a few lived 
beyond 100 days.* Hence, it is possible only to seek a correlation between environ- 
ment and lethal frequency during the time the collected flies were already adults or 
in the late pupal stage. The time period chosen was the week preceding the first 
day of each coilection. The average high, low, and mean temperatures and the 
difference between daily maximum and minimum temperatures (temperature range) 
were calculated for each of the appropriate weeks. The averages thus represent a 
general statement of these environmental components during the particular time 
intervals. 

To evaluate the relationship between lethal frequency and the various tempera- 
ture components, the Spearman rank-correlation test. was used.’ A 2.5 per cent 
level of significance was set and one-tailed tests of significance employed. For 
Chi-square tests on lethal frequency, the usual 5 per cent level of significance was 
used. 

Results.—-Table | contains the data on second chromosomes for the fall of 1958. 
Lethal frequency in September is comparable to past findings, but the allelism fre- 
quency is somewhat lower than the level maintained in the 1947-1952 period.‘ 
In October, though, the lethal frequency has dropped to the lowest recorded to date 
and the amount of allelism nearly so. This 10 per cent reduction in lethal fre- 
quency within less than a month’s time is significant beyond the 5 per cent level by 
Chi-square. When both second and third chromosomes are considered, the decline 
in October frequency is even more significant. Of the 525 chromosomes analyzed 
in September, 199 (38°) were classified as lethal; in October, only 71 (25.8%) of the 
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275 chromosomes were in this category. Chi-square comparison of the two samples 
indicate the decline in lethal frequency to be highly significant beyond the 0.06% 
level. 


TABLE 1 
FREQUENCY OF LETHALS AND ALLELES IN Two 1958 CoLLECTIONS IN SouTH AMHERST, 
MASSACHUSETTS 
Month Tests Lethals % Crosstests Alleles 
Sept. 27: 98 390.9 1386 8 
Oct. E 39 25.8 373 1 
Sept. X Oct. - 961 4 


TABLE 2 


LETHAL FREQUENCY AND Datty TEMPERATURE AVERAGES FOR THE WEEK PRECEDING THE 
COLLEGTION PERIOD 

Chromo- Lethal Temperature averages in °F 

somes frequency 
Year Mo/day tested (%) High Mean Range 
1945 10/8 190 45. ; 56 47 19 
1946 9/29 79 49; : 74 66 18 
1947 9/13 185 38. ie 84 74 20 
1947 10/10 76 39.5 70 58 22 
1948 9/4 148 oo. ‘ 81 69 25 
1948 10/8 162 of. . 53 25 
1949 9/27 198 30.¢ ay 69 34 
1950 9/20 163 < ( 4: 57 27 
1951 9/24; 9/30 182 37 .¢ 5% 65 26 
1952 9/17 190 37 .¢ 56 67 25 
1953 /8 79 7 Hy 66 28 
1953 9/2 158 32. 5 72 60 
1958 9/16 273 35. ¢ E 74 62 
1958 £ 151 2 3$ 52 
1959 7/10; 7/23 118 29. 5 87 , 1 


The short interval between collections, the significant difference in lethal fre- 


quency, and the divergent allelism frequencies suggest that the two samples actually 


represent different populations even though the collection site remained constant. 
Crosstests between lethals in the two collections, however, indicate that both repre- 
sent the same population. Furthermore, one lethal which was present four times 
in the September sample was present once in the October sample also. 

The sudden sharp decline in lethal frequency in what might well have been a 
later sample of the same group of flies suggested a selective force capable of rather 
drastic effects. But as lethal frequency in this population has been declining since 
1947, we hoped to uncover some factor which would prove applicable to the entire 
period, not just the one interlude between samples. A period of cool weather inter- 
vening between the two collections hinted that temperature might be such a factor, 
particularly low temperature. However, the actual day-to-day weather records 
suggested temperature range. This is not unlikely if wide fluctuations between 
maximum and minimum temperatures represent adverse environmental conditions. 

Calculations performed are based on the data given in Table 2. This table con- 
tains the second chromosome data from South Amherst collections for which there 
are corresponding temperature records. Data on lethal frequency from 1945-1952 
have been published previously,‘ except that here each collection in 1947 and 1948 is 
treated separately. The two 1953 collections were analyzed by Ives subsequent to 
that report, while data on the 1958 collections has also been given in Table 1. 
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Analysis of the latest collection, July, 1959, is still in progress; yet it is complete 
enough to be included in the statistical analysis. All but three collections contain 
over 100 chromosomes. Only one, 1951, extended beyond the usual collecting 
period of one to three days. The 1959 collection is actually two—a small early 
July and larger later July—analyzed separately, then pooled. 

tank correlation analysis confirmed the suspected relationship of lethal frequency 
and temperature range. The correlation is —65.7% (Spearman rank statistic 
K = —2.47; P = 0.007). For low, high, and mean temperatures, no significant 
correlations were found. For low, the correlation is 1.1% (K = 0.05; P = 0.487); 
for high, the correlation is —23.9% (K = —0.89; P = 0.187); for mean, the corre- 
lation is —18.6% (K = —0.69; P = 0.245). 
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Fic. 1.—The relationship of lethal frequency (—) to average temperature range (—-—) and mean 
temperature (——) isshown. The temperatures apply to the week preceding collection of the flies. 


Figure 1 presents graphically the relationship of lethal frequency to temperature 
range and to mean temperature. It is easy to see here the significantly negative, 
almost inverse relationship between lethal frequency and temperature range, 
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Somewhat harder to visualize is the negative relationship between mean tempera- 
ture and lethal frequency; erratic ups and downs between 1947 and 1950 are re- 
sponsible for this. It is also interesting to note that during this period, 1947-1950, 
the lethal allelism frequency more than doubled, and it remained consistently higher 
than pre-1947 values‘ until October, 1958. 

Discussion and Conclusions.—Climatie changes of a measurable amount have 
occurred in the northeastern United States during the past two decades.” In 
Amherst, meteorological observations since 1889 at the Massachusetts Agricultural 
Experiment Station show that the annual mean temperature, on the rise since 
about 1905, jumped to a level about 1°C above the average of the previous 60 years 
in a seven-year period beginning with the fall months of 1948. Also, in the 1939 
1957 period there was a moderately different precipitation pattern from that of the 
previous 50 years, lower in winter and summer, higher in spring and late fall. This 
was sufficient to bring about. the introduction of irrigation in favorable agricultural 
locations. In the past three years, both mean temperature and precipitation 
pattern have been more nearly normal. 

Although subtle by human standards, these shifts in climate, especially the 
changes in temperature conditions, may have been instrumental in effecting a strik- 
ing genetic reorganization in the South Amherst D. melanogaster population. The 
evidence from the study of lethal chromosomes suggests that the former relatively 
stable environment has shifted to instability with a consequent reduction in lethal 
frequency, population size, and increase in inbreeding. The upheaval in population 
structure and reorganization at a new genetic level has apparently enabled the 
population to survive, adapt, and flourish once again. From 1950 on, lethal fre- 


quency again follows the behavior of mean temperature, as is expected on the basis 
4 


of the geographical gradient of increasing lethal frequency both in melanogaster*: 
and pseudoobscura.''!: '- (Temperature is thought to be related to this.)'* Hence, 
the major reorganization appears to have occurred between 1949 and 1950. Ives‘ 
reached similar conclusions on the basis of the lethal-frequency—lethal-allelism data. 

Previous reports of such drastic reorganization in the gene pool of a population 
have been confined to inversion-characterized species which invade marginal 
habitats and manage to survive in them.'4~'* In such habitats, inversions are lost, 
so that the population may utilize better the process of recombination and the in- 
finite array of gene arrangements resulting from it. In most populations of D. 
melanogaster, though, recombination is rarely impeded by inversions." Further- 
more, in the present case it is the climate at a given location which has changed, not 
the location. Yet, perhaps it is only a marginal environment which will minimize 
the serious consequences of chance homozygosity while maximizing the genetic 
potential to adapt to a new or altered environment. 

With the successful readjustment to other temperature factors having been 
maintained since 1950, continued marginality of the environment may be ques- 
tioned. Still, the allelism frequencies have yet to return consistently to pre-1947 
levels, indicating either continuation of smaller population size, tolerance for chance 
homozygosis, or both. 

A recent paper has reported a drop in lethal frequency for D. pseudoobscura in 
Texas and California,’* And the decline in lethal frequency in Florida D. melano- 
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gaster populations has already been mentioned.‘ It is not known if climatic 


changes, similar to those in South Amherst, may also have occurred in these areas. 

From the South Amherst population, though, we have obtained both further 
evidence for the dynamic relationship between genetic structure and environment 
and a demonstration that lethal frequency may be under selective environmental 
control just as has been found for inversions. '4~!® 1%~?! 

Summary.—One of the causative agents responsible for reduction in lethal 
frequency in the South Amherst D. melanogaster population has been found to be 
temperature range in the natural environment during the week preceding sample 
collection. It is postulated that the increased temperature range, which has ini- 
tiated a reorganization of the genetic structure of the population, represents a 
marginal situation for D. melanogaster. 
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THE FUNDAMENTAL FREE MODE OF THE EARTH’S INNER CORE 
By Louts B. SLICHTER 
INSTITUTE OF GEOPHYSICS, UNIVERSITY OF CALIFORNIA AT LOS ANGELES 
Communicated December 20, 1960 


The large Chilean earthquake of May 22, 1960, furnished the first fully convincing 
observational evidence concerning the free modes of oscillation of the earth. This 
earthquake excited large response at long periods on both a LaCoste-Romberg 
earth-tide gravimeter! and on Benioff strain seismometers.? In the spectral analy- 
ses by Ness! of the gravimeter observations a sequence of over forty free modes 
has been identified each at a period in accord with the predictions of theory.* 

However, the first spectral peak (with period about 86 minutes) appearing at the 
low frequency end of the spectrum (Fig. 1) was not theoretically predicted. This 
peak is defined by 4 or 5 points and appears both in the “filtered”? spectrum (in 
which the relatively enormous amplitudes due to the semidiurnal and diurnal 
gravity tides were suppressed) and in the corresponding unfiltered spectrum. The 
peak also occurred in the spectral analyses of several of the sub-intervals of the 
1.6-day gravimeter record;' and, significantly, it is absent in the analysis of the 
quiet period (see Fig. 1) a month later. Its “energy’’ (amplitude-squared) is 

POWER SPECTRA ANALYSIS-UGLA EARTH TIDE GRAVIMETER 
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Fig. 1.—Spectral peak (shown by arrow) with period 86 minutes. (Portion of Fig. 4, reference 1, 
reproduced by courtesy of the authors. The numbers above the peaks indicate the period, in 
minutes. ) 
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about twice that in the adjoining spectral regions; the associated increase in am- 
plitude is about 0.64 microgals (1 gal = 1 cm/sec’). It is known that elastic 
oscillations of the inner core*® with a period about 100 minutes may occur, but in 
this mode “. .. the displacements, as well as the perturbations in the gravitational 
potential, are confined to the core .. .”’ 


186 
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Thus an explanation for an 86-minute mode capable of observation by a gravim- 
eter seems to be needed. It is suggested that a gravitational perturbation of the 


approximate observed magnitude and period could reasonably be produced by 


oscillations of the inner core as a rigid body in the surrounding fluid core. Because 
no first-order strain would be associated with such rigid body displacements, this 
mode would escape detection by a strain meter. In the present note some results 
concerning the dynamics of the mode will be summarized, and the significance of 
the mode as a means of providing new information about the inner and outer core 
will be indicated. A more extensive discussion of the dynamics of the core-fluid- 
mantle system* is being offered elsewhere. 

The Three Normal Frequencies.—The mantle (and also a gravimeter on its surface) 
must receive the counter-acceleration required to preserve constant the linear 
momentum of the earth-core system. ‘This associated acceleration of the mantle 
produces a gravimeter reading comparable with the direct gravity effect due to 
the core oscillations, and is additive in phase. For the mode under discussion the 
inner core—fluid core—mantle system behaves as a simple linear vibrator. If the 
circular frequency of this system for displacements parallel to the earth’s rotation 
axis iS wo, then it is easy to show that two eigenvalues, w, and we, exist for oscilla- 
tions executed in the plane of the equator, where a = wy) + Q, w = w — 2. Here 
denotes the angular velocity of the earth’s rotation. Accordingly, if the 86- 
minute period represents the w mode, one would expect to find associated periods 
of 91 minutes and 81 minutes. At present the resolution of these three periods has 
not been accomplished; but with the aid of more extensive obseravtions it may be 
feasible to identify the three associated modes. 

Energy and Momentum in Core Mode.—The amplitude of the core’s displacement 
required to produce the 0.64-microgal perturbation associated with the 86-minute 
period is about 48 cm. The energy associated with a core vibration of this ampli- 
tude, at period 86 minutes, is about 10°% ergs. The energy released in the Chilean 
earthquakes of last May in the form of seismic waves is still uncertain, but pre- 
liminary estimates (C. P. Richter, private communication) indicate that this total 
energy may be of order 10*> ergs. Thus the required energy in the core vibration 
is possibly 1 per cent of the energy released in the Chilean quakes. One may also 
examine the adequacy of the momentum transferred to the inner solid core from 
incident seismic waves. Partly because of the excellent match in acoustical im- 
pedance at the boundary between outer core and mantle, coupled with a severe 
mismatch at the inner core, the momentum transferred from incident dompressional 
seismic waves appears to be quantitatively adequate for exciting the postulated 
amplitude of oscillation of the inner core. 

The Gravitational Restoring Force-—The gravitational force restoring an inner 
core of density po, whose center is displaced a distance 6 from the geometric center 
of mantle and outer core, to its central position of equilibrium within a sphere of 
density ps, is precisely the attraction between two spheres with centers separated 
by 6, one having the size of the inner core, but of differential density pp — p;, and the 
other of density p;, and radius 6.¢ Since the mass of the latter sphere is proportional 
to 6°, and the attraction varies as 6~?, this force, then, is proportional to the separa- 
tion distance 6. If the outer core be regarded as a perfect incompressible fluid, 
it is found that the natural period of such a system is nearly independent of the 
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density of the outer core, in the broad range 9 < p; < 14, but is a rapidly decreasing 
function of the density of the inner core. To achieve as low a period as 86 minutes 
in such a system, a density po = 27 to 29 is needed, which is much higher than the 
limits 12 to 22 gm/em* usually assigned to the inner core. Granting even a core 
period as long as 100 minutes, an inner core density of 20 gm/cm* or more would 
be needed if only gravitational restoring forces are present. 

Need of Additional Restoring Forces.—To avoid the need of postulating so high a 
density in the inner core, one may assume the presence of additional restoring 
forces of other types. For example, a relatively small value of order » = 10'° 
dyne/cm? for the elastic modulus of rigidity in the outer core is competent to pro- 
vide the needed additional restoring force, while retaining more usual concepts of 
the inner core’s density. Values » < 10'! dyne/cm? for the shear modulus are re- 
quired to satisfy some recent observations concerning the reflection of shear waves 
from the mantle-outer core boundary.® 

More specifically, the equation of motion for the mode w) (whose displacements 
are in the direction of the rotation axis) is 


mo(1 + epipo!)(6 + 2Qaé) + ) 16Bayu + 4/3rymo pi (1 — pipo 1) }6 = 1) 


where mp and dp are the mass and radius, respectively, of the inner core, @ is the 
damping coefficient, y is Newton’s constant, and ¢ is a coefficient having value nearly 
'/> which takes account of the increased dynamic mass of the inner core. Its value 
depends slightly on the size of the fluid core and upon the variation of density in its 
interior. The value of the coefficient 8 is also nearly '/2. On the basis of this 
equation, if pp = 15, p: = 10, the value of u required to produce an 86-minute period 
is 4 X 10". In this case the elastic forces contribute 53 per cent of the restoring 
forces. The maximum resultant shear stress at the surface of the inner core is very 
small—4.8 X 104 dyne/em? or about 50 grams per square centimeter. It is sug- 
gested that the core fluid is capable of sustaining shear stresses of this small mag- 
nitude elastically, throughout the necessary time interval of the order of an hour. 
In comparison, the duration of the stresses associated with the slow motions in the 
fluid core needed to account for the earth’s magnetic field is extremely long. Namely, 
the time intervals used in describing secular variations of the earth’s internal mag- 
netic field are hundreds of years rather than hours. The velocity of the fluid as- 
sociated with secular magnetic changes is of order several hundredths of a em/see ;6 


9 


the maximum local velocity for the core mode under discussion is about 6 & 10 
cm/sec, i.e. up by a factorof2or3. Thus the difference in the velocity values in the 
two cases is not great and may not be invoked to explain the postulated Cual elastic 
and fluid-like behavior. The postulated elastic behavior of the fluid at periods of 
1'/, hours and its fluid-like behavior during prolonged stress applica‘) i~ is, how- 
ever, consistent with familiar concepts relative to the visco-elastic properties of some 
common materials. 

Absence of Seismic Shear Waves in Core.—Granting the suggested small value 
of the elastic shear modulus, then the nonpropagation of shear waves through the 
core may be ascribed to two causes. In the first place, the low value w» = 10!° 
with p; = 10, implies a velocity for shear waves of only 0.63 km/see. The as- 
sociated impedance mismatch at the outer core boundary implies that only 20 


per cent of the incident shear-wave energy would be transmitted into the core and 
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out again, even under the most favorable circumstance of normal incidence. sec- 


ondly, at. this slow velocity, waves traversing a significant distance in the core 


would be greatly delayed, and would arrive (with small amplitude) in the late part 
of the record, disturbed by surface waves and the complex of late arrivals. Under 
these circumstances unexpected shear phases might well fail to be recognized. 
Energy Dissipation in the Core. The rate of energy dissipation in a vibrating 
system is commonly expressed by the value of its Q factor. This factor is simply 
related to the damping coefficient @ in eq. (1), ie. Q = w/2mr. From spectral analy- 
ses showing the earth’s natural modes the Q factor is best determined from ratios 
of the observed amplitudes of the pertinent peak in successive time intervals of 
observation.! Using spectral amplitudes of the core mode for three intervals of 
observation (each of 56.7 hours duration) beginning, respectively, about 6 hours, 
33 hours, and 60 hours after the main shock, this method yielded the value, Q 
128.t This value which characterizes these pure shear strains confined to the 


) values obtained for high- 


“fluid” outer core is as large as 35 to 50 per cent of the Q 
order modes whose strains largely reside in the solid mantle. For example, the 
high-order mode / 15, whose period is 8.37 minutes or less than |)» that of the 
fundamental core mode, has a Q of about 250-350. ! 

Significance. The present observations concern the period, amplitude, and 
energy loss of a mode whose deformations are confined to the “fluid”? outer core. 
They suggest relationships indicated in eq. (1) among the values of the density of 
the inner and outer core, and the value of the elastic shear modulus of the outer 
core; the observations also reveal the small value of the damping coefficient a, and 
therefore the small degree of an-elasticity of this material when subject to small 
shear strains. The study of this inner-core oscillation provides new observational 
evidence about the earth’s central region, where adequate methods of observation 
are scarce indeed. The identification of the material of the core is still uncertain. 
It has recently been suggested that the outer core is a mixture of iron and silicon.’ 
Perhaps the major significance of studies of the core-mode will be the development 
of information concerning physical properties of the outer core which will assist in 
identifying the material of this important feature of the earth’s evolution. One 
may here observe mechanical properties of the earth’s core under a combination of 
conditions of high pressure (several million atmospheres) and high temperature 
(about 4,000°C) which is not attainable in the laboratory. 


The author desires to express appreciation to the Office of Naval Research for 
many years of support of gravitational studies at the University of California Insti- 
tute of Geophysics, and to his colleagues, Norman, Ness and J.C. Harrison. The 
former has significantly improved techniques of spectral analysis by the introdue- 
tion of the effective “Chebyshev” type of “window,” and with this method has 
made spectral analyses of the gravimeter data shown in part in Figure 1. J.C. 
Harrison developed the automatic digital recording methods used and has generally 
supervised the experimental aspects of the project. The instrument itself, which 
fundamentally consists merely of a weight suspended by a spring, yet which is ca- 
pable of sensing changes of one part in 10 billion in the local value of gravity, is the 


result of nearly thirty years of experience in the design and construction of sensi- 
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tive gravimeters by the LaCoste-Rombert Company.' Their cooperation in sup- 


plying an instrument of the amazing quality needed has obviously been basie in 
our studies of the faint manifestations of the earth’s free oscillations. 


«The assumed system is briefly the following: (1) An inner rigid core, of radius 1,250 kin 
2) a surrounding compliant spherical shell of outer radius 3,473 km; (3) the surrounding rigid 
spherical shell of the mantle, of outer radius 6,371 km. The rotation of the earth is obviously 
pertinent in the dynamics of this svstem. 
+ This result follows from well-known gravitational properties of spherical mass distributions and 
may be readily verified ab initio for the geometrical relationships that exist for the displaced core 
The Q value depends upon the noise level at normal times. This level which was read as 
300 (microgals)?/eph for the quiet period a month later was adjusted slightly to .298 to secure 
precise agreement between the two independent Q values derivable from the three intervals of 
observations 
Ness, N. FL, J.C. Harrison, and L. B. Slichter, /. Geophys. Research, 65, 2 (1961 
Jenioff, H., F. Press, and 8S. Smith, /. Geophys. Research, 65, 2 (1961 
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THE CATEGORY OF TOPOLOGICAL OBJECTS* 
By Yearon H. Ciirron anp J. WOLFGANG SMITH 
DEPARTMENT OF MATHEMATICS, MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
Communicated by Paul A. Smith, December 14, 1960 


|. Introduction. -In this note, we shall introduce a certain category @ which 
extends the category C of topological spaces. This extension was originally 
motivated by considerations of the following type, pertaining to foliated manifolds 
Given such a manifold, one finds that every point has a neighborhood in which 
the foliation is a fibering, and this circumstance permits one to construct local 
base spaces which are again manifolds. If one attempts to piece together these 
local base spaces to obtain a global base space, one will in general be forced to 
make additional identifications which destroy the manifold structure of the space 
The concept of topological object defined in section 2 has resulted from an attempt 
to preserve the local manifold structure and yet to obtain something which plays 
the role of a base space. Concomitant with the new concept of object, one requires 
a new concept of map, and these together define the category @ referred to above. 
KXvery foliation determines an object in @, and this object will be isomorphic to 
the ordinary base space whenever the foliation is a fiber space. In other cases, 
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the object cannot be regarded as a point set. Besides foliations, many other 
mathematical structures (e.g. Lie groups and homogeneous spaces) give rise to 
objects in @, and these are generally not topological spaces. 

We shall show how every subcategory Cy of C can be canonically extended to a 
larger subcategory @) of @. Generally speaking, one finds that a mathematical 
theory defined on a category Cy admits a rather natural extension to the category 
Cy. Concepts pertaining to homotopy, for example, are meaningful on @. Again 


many concepts belonging to the realm of differential geometry can be defined on 


(”-manifold objects (i.e. objects in @o, Co being the category of C’-manifolds and 
C’-maps). In the present note, we shall illustrate the problem of extending 
classical theories from a subcategory Cy to the subcategory @y by considering co- 
homology theories arising from a sheaf-valued functor. One finds in this case 
not only that the extension can be carried through in a natural way but, moreover, 
that it leads to some additional cohomology invariants, these being trivial for 
objects in Cy. It is consequently possible, by means of a cohomolegy theory, to 
measure the extent to which an object in @p deviates from a topological space. 

Our results concerning cohomology suggest that the extension of classical 
theories to a subcategory of @ will lead in general to some new phenomena. An 
application of the cohomology theory to Lorentz geometry will be indicated in the 
last section of this note. 

2. The Category @.—Let X and X’ be topological spaces. By a local map 
f:X — NX’, we shall understand a continuous function whose domain D(f) is an 
open subset of Y and which takes values in X’. Given two local maps f, g:X —~ X’ 
and a point « < D(f) n D(g), we shall say that f = g near x if there exists an open 
neighborhood U of x such that f,U = g U. We now define a preobject to be a pair 

X,/), where J is a set of local maps u:X — X satisfying the following extension 
axiom: Given ueeJ and xe Diu) n D(v), there exist a, 72/ such that wou = gov 
near 2. A topological object X (henceforth referred to as an object, for the sake of 
brevity) is a preobject (Y,/) satisfying the following conditions: (7) The identity 
map ty:X — NX belongs to /, and (ii) if @:X — X is a local map such that for 
every « ¢ D(a) there exist u’v,o’ ¢ 7 such that wu’ o a = v’ ovneara, thentiz/. We 
observe that condition (ii) implies closure of 7 under composition, restriction and 
union of local maps, as well as closure under a kind of factorization. It may 
consequently be referred to as the closure axiom for objects. An object (Y/) 
is said to be generated by a preobject (X,/) if J ¢ J, and if, for every object (XY J’) 
with J ¢ J’, 7 is a subset of J’. One can now prove 

LeMMA 1. A preobject generates a unique object. 

Our next task is to specify what is meant by a map F:% — X’ when X& = (YX,/) 
and x’ = (X’/’) are objects. We first define a premap G:% — 9%’ to be a set of 
local maps f:X — X’ such that 

U wD(f) = X, 
and satisfying the following extension axiom: Given uel; fg2eG,andx2D(u) a 

D(g) such that u(x) ¢ D(f), there exist w’,v’ J’ such that u’ of o u v’ og near x. 
A map F:% — &’ is now defined to be a premap satisfying the following closure 
axtom: If f:X — X’ is a local map such that for every x ¢ D(f), there exist uJ, 
fe F, and uw’, v’ 21’ such that v’ of = uw’ ofouwnears,thenfeF. The efficacy of 


this closure axiom for maps is attested by 
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LemMa 2. Jf F,G:% — X’ are maps and F c G, then F = G. 

A map F is said to be generated by a premap Gif Gc F. One finds 

LEMMA 3. <A premap generates a unique map. 

It remains to define the composition of maps. For this purpose, one proves 

LemMa 4. Let F:9% — X' and F’:%' + 9” be maps. Let H denote the set of all 
maps f’ of, f¢ F, and f’ ¢ F’, which are well defined. Then H:X% — %” is a premap. 

By virtue of Lemmas 3 and 4 we may define the composite F’ o F:% — X’ to be 
the map generated by H. One observes now that if 9% = (X,/) is an object, then 
1:% — & is a map, and moreover, J functions as an identity under composition. 
We shall call it the identity map of 2. 

If X is a topological space, then (X, }zx}) is clearly a preobject. Let 9x denote 
the object generated by (X, }ix}). In this manner, every topological space may 
be regarded as an object. Similarly, if f:X — X’ is a (global) map, then jf}: 
Xy > Xx’ isa premap. In this manner, every map f:X — NX’ may be regarded 
as a map F';:%, > Xx’. We can now state 

THEOREM 1. The collection @ of all objects X, ’, ... and all maps F:% — Xx’ 
constitutes a category.' Moreover, the category C of topological spaces (regarded as a 


subclass of @ under the identification described above) constitutes a subcategory of @. 
If xc = (X,J), 2’ = (X’ J’) are objects and F:9 > x’, F’:9’ — X are maps 

such that F’ o F = J and F 0 F’ = I’, then F is called an isomorphism. It is im- 

portant to note that an isomorphism F in general does not determine a one-to-one 


correspondence between the points of X and X’. On the other hand, we shall be 
concerned only with those propertiesof objects whichare invariant under isomorphism. 
It follows that a general topological object cannot be regarded as a point set. 

Let Cy be a subcategory of C. Given an object X = (X,/), let J» denote the set 
of all we J such that we Cy. Similarly, given a map F':% > &’, let Fy denote the 
set of all f ¢F such that f¢ (Cy. Let @o denote the class of all objects © = (X,/) 
such that Y ¢ Cy and (X,/o) generates %, together with the class of all maps F' such 
that fy isa premap. One easily verifies that @) is a subcategory of @, containing 
Cy. Let @o* denote the class }9} of all objects in @y which are isomorphic in Cy 
to objects in Cy, together with the class }F} of all maps F':9 — 9%’ in @p such that 
x, x’ ¢ }xX}. One finds again that @y* is a subcategory of Cy which extends (%. 
The objects of @* (i.e. all objects isomorphic with topological spaces) will be 
referred to as quasi-spaces. 

3. The Functor @—We will now define a covariant functor! @ from @ to C, 
which plays an important role in the theory of objects. Let 9% = (X,/) be an 
object, and let Q denote the set of all pairs (x,y) ¢ X &K X such that u(x) = v(y) 
for some (uv) ¢l XK J. It is easily verified that Q is an equivalence relation on X. 
Let X denote the set of equivalence classes in X, endowed with the quotient topol- 
ogy, and let p:X — 9X denote the natural projection. Then }p}:% > & is a 
premap. We shall refer to X as the basic space of the object 9%, and to the map 
Px generated by } pj as the basic map of 9. 

THEOREM 2. Given the map F:% — %’, there exists a unique map F':9% —> x’ 


A A A 


making the diagram 
P 
x —— x’ 
Py} p 4 Px 
oo 





Vou. 47, 1961 MATHEMATICS: CLIFTON AND SMITH 193 


commutative. The function 0:@ — C defined by setting (2) = X and P(F) = F 


constitutes a functor. 

Coro._iary. An object is a quasi-space if and only if its basic map ts an isomorphism. 

1. Cohomology.—Let Cy denote a subcategory of C, and let S denote a functor 
which to every space X in Co assigns a differential sheaf? S(X) on X and to every 
local map f:X — X’ in Cy assigns a homomorphism S(f):S(X’) ~ S(X) compatible 
with f. We shall construct a corresponding functor $, which to every object 
in @o assigns a differential sheaf $(9C) on P(X) and to every map F: XL —> XL’ in Cy 
assigns a homomorphism S(/’):8(%2’) — $8(%) compatible with @(F). Let © = 
(XJ) be an object in @y, and let U be an open subset of @(9). A section 
a:p-'(U) + S(X) will be called an /-section if S(u) 0 ao u = aon D(u) forall we Io 
with Diu) ¢ p-'(U). The totality of these /-sections constitutes a differential 
module, which we denote by T(S, X, U). If V ¢ U is open, the inclusion map 
i:p-'(V) — p~'(U) induces a homomorphism 7*:T(S, 9%, U) > ius, x, V). This 
defines a presheaf on (9), and we will let $(9%) denote the corresponding sheaf. 

Now let F:X — &’ be a map in @o, where X = (X,/). To define a map S(F): 
$(9¢’) — S8(2), we consider an open set U’ ¢ (&(2’), its inverse image UL’ under 
W(F), and a section o’ ¢ T(S, x’, WU’). The possibility of pulling «’ back by F to 
and /-section over p~!(U’) results from the following lemma, which in turn depends 
on the extension axiom for maps. 

LemMa 5. If wely and fae Fy, then S(fou)oa’o (fou) = S(g) oc’ ogon 
Dif ou) an Dig) a p“(U). 

The lemma permits us, in the first place, to define a section «:p~'(U) — S(X) 
by the formula ¢ = S(f) o a’ of, where f is understood to range over all maps in 
Fy such that D(f) ¢ p~'(U). The lemma also shows that ¢¢T(S, X%, U). This 
gives rise to a natural homomorphism S(F’) :8(2%’) — 8(%X) compatible with @(F). 

THEOREM 3. The function § defined above is a functor, extending S. 

Let I denote the functor which to every sheaf assigns its module of global sections 
and to every sheaf homomorphism assigns the induced map of sections. Let H 
denote the functor which to every differential sheaf (differential module) assigns 
its cohomology sheaf (cohomology module) and to every sheaf homomorphism 
(module homomorphism) assigns the corresponding homomorphism on the co- 
homology level. One now obtains the cohomology functors H o $8, T o H o Sand 
H oT o 8, defined on the category @o, the last one being an extension of the classical 
cohomology functor HoT oS arising from S. 

We recall that in standard cases, H o S is trivial (‘‘local triviality”’ condition). 
On the other hand, this does not imply that H o Sis trivial on @o, as will presently 
be shown by means of a simple example. Consequently H o $ and To H o $8 may 
be used to measure the extent to which an object in @p is not of class @o*. An 
application of this principle is given in section 6. 

It is easy to prove that there exists a natural homomorphism y¥:H oT 0 S(%) > 
To H o 8(X), X being an object in Cp. 

Lemma 6. Let F:% — X’ be a map in Cy. The resulting diagram 

H of o8(X’) - > Holo S(X) 
v' ty 
To o SX’) —— ToH 0 S(X) 


is commutative. 
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It follows that Ker y and Coker y determine two additional cohomology functors 
on @o, the former being an extension of the classical reduced cohomology functor 
arising from S. 

To illustrate the preceding discussion, let X denote a 2-dimensional cylinder, 


parametrized by a linear coordinate x and a circular coordinate 6. Let J denote the 
set of all local maps f:X — X such that f(x, 6) = (2 + a, @), where a is a real 
number depending on f. Clearly X% = (X,/) constitutes an object, whose basic 
space is a circle. We take S to be the functor which to every differentiable mani- 


fold M associates the sheaf of differential forms on M. A simple calculation 
indicates X has the H o To $ cohomology of a torus and the T o H o $ cohomology 
of a circle. 

5. Foliations——Let M be an n-dimensional manifold foliated* by leaves of 
dimension r (O<r<n). Let f = (fi, fo) be a homeomorphism from an open subset 
D(f) of M onto R" XK R"~" (where R* denotes k-dimensional real number space). 
f will be called a distinguished coordinate system of M if, for every y ¢ R"~', fo“ '(y) 
isa subset of a leaf. Let & denote the set of distinguished coordinate systems of M. 
For every f ¢ ®, let J, be the set of local maps u: D(f) > D(f) such that fo(u(x)) = 
fo(x) for alla e Diu). Let J = Upely Then (M,/) is a preobject. The object 
ym = (M,J) generated by (M,/) is called the canonical object of the foliation. 

It is easy to see that the basic space of IW can be canonically identified with the 
set of leaves in M, endowed with the quotient topology. If the foliation is a fiber 
space over a base space X, IW will be a quasi-space isomorphic with X. More 
generally, one can prove 

THEOREM 4. A foliation is regular’ if and only if its canonical object is a quasi- 
space. 

We would like to point out in this connection that B. L. Reinhart has apparently 
considered cohomology modules of type H oT o § (3M) associated with a foliation 
of M.4 

6. Lorentz Manifolds—Let V be a Lorentz manifold, i.e. a Riemannian mani- 
fold of dimension n > 1, whose metric tensor has signature (+ — ... —). Let 
B denote the sub-bundle of the tangent bundle of V, consisting of all unit time-like 
vectors (i.e. vectors v such that v-v = +1). If on B one identifies every vector 
v with —v, one obtains a new bundle M on V. Every time-like geodesic on V 
admits a canonical lifting to 17, and this defines a foliation on WM. Let 9M denote 
the canonical object of the foliation, and let S denote again the functor which to 
every differentiable manifold assigns its sheaf of differential forms. The results 
of Section 4 lead thus to a cohomology theory of Lorentz manifolds, which is 
analogous to the homotopy theory previously considered by one of the authors.° 

It can be shown that the nonexistence of closed time-like paths on V implies 
regularity of the foliation on M/. This leads by Theorem 4 to the following result : 

THEOREM 5. Nontriviality of the cohomology sheaf H o § (MN) implies existence 
of closed time-like paths on V. 

In conclusion we shall state a second result relating the new cohomology in- 
variants to geometric properties of V. 

THeoreM 6. Jf V is flat and complete and does not admit closed time-like paths, 
then H o T 0 S(M) ts naturally isomorphic to the de Rham cohomology of V. 
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DEGREE OF APPROXIMATION TO A CONTINUOUS FUNCTION ON A 
NON-ANALYTIC CURVE* 
By W. E. Seweiit 
DUKE UNIVERSITY 
Communicated by J. L. Walsh, December 19, 1960 


1. Jntroduction.—The purpose of this paper is to supply the proof of some re- 
sults which were announced in these PROCEEDINGS! more than twenty years ago. 
At that time the proof was to “appear elsewhere.’’ This is a belated appearance 
and the proof is a later (1960 rather than 1937) model, because the manuscript 
containing the original details has “disappeared” in the interim. The present proof 
does serve to establish a result which is slightly more general. 


We begin with two definitions. 


DerriniTion A. Let C by a Jordan curve in the z-plane and let C denote the close d 
limited set bounded by C. Let the function f(z) be analytic in C and continuous 
in C. Then we say that f(z) belongs to the class L(k,a) on C, where k is zero or a 
positive integer, and 0 <a S 1, if f(z), the k-th derivative (f(z) = f(z)) of f(z), 


exists on C in the one-dimensional sense and satisfies a Lipschitz condition of order a 
on C. 

It should be noted here that the derivatives are taken with respect to z on C 
and not with respect to are-length on C. Since C is not analytie there is a difference 
between the two for the higher derivatives. Also the Lipschitz inequality refers, 
of course, to the chord and not to the are on C. 

Now let C consist of one or more mutually exterior Jordan curves (1, C2 
in the z2(= x + zy)-plane, let C denote the closed limited point set bounded by C, 
and let K denote the complement of C with respect to the extended plane. Then? 
K possesses a Green’s function G(z,y) with pole at infinity, and the function w = 
o(z) =e + 4 
sarily uniformly onto the exterior of the unit circle y in the w-plane so that the points 


, where H/ is conjugate to G in K, maps K conformally but not neces- 


at infinity in the two planes correspond to each other; interior points of K corre- 
spond to exterior points of y, and exterior points of y correspond to interior points 
of K. Under the transformation w = $(z) we can, by choosing a particular branch, 


set up a correspondence between separate arcs y;,j = 1,2, ...,, of y and the curves 
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C,of C. If we preassign points P; of C;, take each are y, closed at one end and open 
at the other end, and set the closed end of y; in correspondence to P; and the 
closed end of y; + 1 in correspondence to P; + 1, we obtain a one-to-one correspond- 
ence between C and y. Now let z = y¥(w) be the inverse of the function w = 
¢(z) which defines this correspondence. Then the function ¥(w) is continuous at 
every interior point of each y,, whereas the end-points are points of discontinuity of 
the function. Our second definition concerns the behavior of the function Y(w) on 
the unit circle. 
DEFINITION B. We say that the set C is of type W, provided 
octas ly’ (w) EE a 
and 
W"(w) is continuous 


at each point of continuity of Y(w) on |\w) = 1. 


We will prove the following: 
TuHeoreM. Let the point set S consist of a finite number of limited closed regions 
‘\, of Type W. 
Let the functions f;(z),j = 1,2, ...,\, belong to the class L(k,a) on C;, respectively. 
Let f(z) =f;(z) where zisin C;,j = 1,2, ...,r. Then for eachn,n = 2, 3, .., there 


bounded by the set of mutually exterior Jordan curves C:C), Co, ..., € 


exists a polynomial P,,(z) of degree nin z such that 


log n\*** 
T(z) a P,(z — M ,ZOnen, 
n 


where M is independent of n and z. 
Our method is a natural extension of Walsh’s* proof that a function which is 


analytic in the interior of a Jordan curve C and continuous in C can be uniformly 
approximated in C by polynomials in z. For a function f(z) belonging to L(k,a) on 
S, where the C, are analytic Jordan curves, J. H. Curtiss! proved that there exist 


polynomials such that | f(z) — P,(2)| < M/n***, zon S. Mergelyan’ has proved 


that if f(z) belongs to L(k,a) on a curve C essentially of Type W and ¢« > 0, then 
P_,(z) exists such that | f(z) — P,(z)| S$ M(6)/n***7* zon €. 
2. The case k = 0, = 1.—We begin with a single Jordan curve C of Type W 


and suppose without loss of generality, that C contains the point z = 0 in its in- 
terior. Now let C,,p > 1, denote the image in the z-plane of jw! = pin the w-plane 
under the transformation z = y¥(w). For our purpose we will take p near 1 so 
that p — 1 is near zero. 

Certain properties of C and C, are needed in the proof of the theorem. 

(i) There is a constant D, independent of p, such that every point of C, is contained 
in a circle of radius D(p — 1) about some point of C, and every point of C is contained 
ina circle of radius D(p — 1) about some point of C,. This follows from the fact that 
by (1) the function p(w) satisfies on |w| = 1 a Lipschitz condition of order unity, and 
by virtue of tts continuity tt also satisfies’ the same conditionin 1 S \w| <p> 1. 

(ii) © and C, have continuously turning tangents. This follows from the con- 
tinuity of y'(w) on |\w = 1, which is a consequence of the existence of ''(w) there. 
And C, 

(iii) Jf As denotes the (shorter) arc of C or C, between 2, and 22, both points lying 


is analytic. 


p 


either on C or on Oe then 
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As 
=< mM, 
ai = Be 


where m is a constant independent of 2, and 22. 
(2) by considering the quotient dz/\ dz 
(iv) 


This follows (S, €5-66) from (1) and 
, and noting that (1) holds for jw| 2 1 


If we write the equations for C and C, in parametric form as follows: 


C: 2(0) = ye"), C,: z,(0) = ¥(pe®), w = pe” 
then directly from (1) we have 


~ 


(0)|} 2a>O0,\z',(0)) 2a> 0. 


Also, since y(w) satisfies a uniform Lipschitz condition in w) 2 1, we have 


2(@) — 2,(0)| & Dip — 1). 
And, since from (2) 


we know that y'(w) satisfies a uniform Lipschitz condition in 
w| 2 1, we have 


2/(@) — z,'(0) < D(p — 1), 
using the fact that 


dz ; dw : Zs . : ; j 

~ = w’(w) y'(w) pie”, p 2 1 and constant, and that y'(w) is uniformly 
dé dé 
bounded in modulus for |\w| 2 


l 


(v) If r(s) and r,(s) denote the tangent angles (with the positive real axis) of the 
tangents to C and C,,, expressed as functions of the are length, then 


net) = 


r(s)' S Ht, |7,(s + 8) T,(s)| S Ht, 
where H is a constant independent of s and t. This follows from the fact that dr/ds 
and dr,/ds are uniformly bounded in modulus by virtue of (1) and (2). 

Now let 2; = x,(z) map the interior of C, in the z-plane on the interior of C in the 
z,-plane so that x,(0) = 0, x,’(0) > 0. Then by virtue of i, ii, iii, iv, and v, we know 
from a theorem of Warschawski® that 


z— x,(z)} S Ni(p — 1), zin C, (4) 
where N; isa constant independent of p and z 


If f(z) belongs to L(0,a) on C, then (8, 


17) it satisfies a Lipschitz condition uni- 
formly in the closed region, hence by (4) we have 


d(z) — f(x,(2))| S LANi*(p — 1)%, z in €, 


(5) 
where 1; is the Lipschitz constant. 
We know that f(x,(z)) = F(z) is analytic interior to C,, and 


1 continuous and 
uniformly bounded in modulus in C,, consequently we can use the Hermite Inter- 


polation Formula (see, e.g., 8, 43 
I on(z)F (t)dt . 
. eine, 
2mt J cp w,(t)(t — 2) 


, 3, 4) 


where 
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2, = 1,2, ..., + 1, lying interior to C,, and P,(z) is the unique polynomial of 
degree n which interpolates to F(z) in the points 2;. 
Here we choose as the points 2; of interpolation the transforms of all the (7 + 1)-th 


roots of unity under the transformation z = ¥(w). Then we have from (6) 

1 
TyTel, zin c (8) 

Tv 


y) 


- 


F(z) — P(z)| 


where 


may |. (2) 
I, ={zonC |— : (9) 
ton C, wn(t) 


max 
I, = ( F(t) : (10) 
ton C, 


ldt 
b= (fit jzon) (11) 


The problem is to find upper bounds for J;, I2, J;, and hence for |F,(z) — P,(z)| 
for zon C. 

By virtue of (1) and (2) we know’ (see, e.g., 8, 48 ff.) that 

] ul w —1) , 
1=> We = » |W 
"Ww" — 1 

where 7, depends on ¥(w), but not on p, n, z, and ¢. 

Also by the construction of F(z) we know that 


I, < My, (13) 


where /, denotes the maximum modulus of the original function f(z) on the curve 
ic: 
Furthermore, applying (1) and (2) and moving to the w-plane, we have 


f dt oe f y af dw < A, f c-. -" 
Cp |t — 2 wi=p |\W(w) — ¥(w)| 0 |pe —e”*| 


where w’ = pe”, w = e’®, o fixed, and A; is independent of pfor pnear1. Letting 


¢ = 0, without loss of generality, we have 


{- pde | dé 
6 = £ 
ow = 0 V1+ p? — 2p cos 6 


n i" , dé in F dé \ 
= <p < = . 
Jo Vj + p? — 2pcos 6 AVY + = Spent 
The second integral in brackets is not greater than (32/8)/{sin (#/8)} since 
pe” — 1|> 2sin (r/8) for r/4 < 6 < xr. To bound the first integral in brackets we 
use the fact that 


4 ’ : ; 6? 64 66 
1 + p? — 2p cos 8 i+ ~eetin~ ot nm ee 


Z 6! 
2 p64 


oe ” ee 
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26? 
= lp I) + (1 1- ) 


, llp , 
(o— 1)? + 75 6,0 <0 


Hence 


Vo 
i Jas — 2p cos 0 sy llp \? 
eo Ey ‘) 
12 
< M,g\log (p — 1)i, 
where M, is a constant independent of pfor pnear 1. Consequently 
I; S M,\log (p — 1)|, 


where M, is independent of p, z, and t, when suitably chosen for p near 1. 
Therefore from inequalities (5), (8), (12), (13), and (14) we obtain 


le 


) J 
fle) — Pal2)| $ Lo — 1)" + ——*— MMyllog (o — 1)|, 
a 


and if we put p = 1 + (log n)/n we have 


' log n\* M, we 
oe S * = —_____ MM, og ("8 ‘) 
n j log n ) | 
(: * ) “— 


( 


n 


n 


log n\* 
<M ( "8 ") gon. (15) 
n 


for n sufficiently large, where M is independent of z. Then by the principle of the 
maximum the theorem is established for the case k = 0, and a single curve. 

3. The case k > 0, = 1.—-We continue to restrict S to a point set bounded by a 
single Jordan curve. Let us consider first the case k = 1. Since f’(z) satisfies a 
Lipschitz condition of order a, we have by the results of paragraph 2, above, 


log n\* = 
S'(2) — pa’ (e)i 3 u,("% *) ,zon C, (16) 


log n\* log n 
L + M'M,M, 
n 


n 


where p,,’(z) denotes the polynomial of degree n approximating to f’(z). Now form 
the function 


Ta+1(Z) = f(z) — pasi(z), (17) 


4 
Dn +1(2) -f Dn’ (z)dz 
Zo 


where 
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is a polynomial of degree n + 1 in z, with 2 fixed on C and the path of integration 
on (since Cis of Type W) or within C. By (17) we know from (16) that 


log n\* = 
< M, ( = "| gon-c, (18) 


n 
hence 


log n 


rn + 1(2) — Yn + 1(x,(2))| S un ( ) M.(p — 1), z on ‘ad 


n 


Now let 7, + 1(x,(2)) = Rnp(z), zon C,, then 


log n\* = ? 
Rup’ (z)| = |rn + 1'(xp(2)) (x,’(2))! & Ms ( S zon Cp, (20) 


n 
since x,/(z) is uniformly bounded in modulus for p near 1 (actually x,’(z) approaches 
1 as p approaches | from Warschawski (loc. cit.) or simply by the properties (1) and 
(2) of C and the fact that C, approaches C). Consequently, using the Hermite 
Formula as above, to interpolate to R |,,'(z), we have 

— fmax 
|Rnp'(2) — Qn(2) 5 [zon [Rup (2)| ] M,|log (p — 1)}, 
C 


p 


log n\* flog n 
< M( S ) ( S J, zone, 
n n 


by virtue of (20) and letting p = 1+ (logn)/n. Now form the polynomial 

P, + (2) = { Q,(t) dt + Rr,(20), 
with z on C and the path of integration suitably chosen. Since R,,,(z) — Px + 1(20) 
= (), we have 


IR»(2) — Pr+ r(2)| = if {Rny’(t) — Qn(t) }dt| 


0 


- a+l 
< M; (" ") ron C, (22) 


n 


where M; is independent of n and z. Combining inequalities (19) and (22), with 
p = 1+ (log n)/n, we obtain 


a 
zon C, 


log (n + “ 
n+ 1 ; 


| f(z) — On + i(2)| $M ( 


where ®, + (2) = pn + 1(2) + Py, + 1(2) is a polynomial of degree n + 1. This 
completes the proof fork = 1 and Ca single curve. 

For k > la repetition of the steps above yields the general result. 

1. The case of several curves.—¥Yor each curve Cj, 7 = 1, 2, ..., A, of the set S, 
the C; being mutually exterior and of Type W, we have proved that 


k+a 
log m ty 
<M; ( = , Pons, (23) 
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where f;(z) belongs to L(k,a) on C), and P,;(z) is a polynomial of degree m in z. 


From (23) and the properties of f;(z) it follows that 
Fu) Si, 2 Ons P= 1, 2. fe, ? (24) 


Now let z = ¥(w) be the function (of paragraph 1) which maps |w| > 1 conform- 
ally on K, the complement of S, so that the points at © correspond to each other. 
Then there exists a number R > 1 such that for every r, 1 <r S R, the circle |w) = 
r goes into A mutually exterior Jordan curves C, containing S in their interiors. It is 
well known that inequality (24) implies 


Pj(2)| S&S Mr™, z on C,. 


Using the Hermite Formula, we have 


1 f an(z) Pnj(z)dt bn 
Pnj(z) — Pn(Z) ,zons, (20) 
7 tf 


2m w,(t) (t — 2) 


J = 1,2, ..., A, where &,(z) is the polynomial of degree n which interpolates to the 
Pmj(Z) inn + 1 points on 8S. If we choose as the points of interpolation the trans- 
forms all the (n + 1)-th roots of unity (the correspondence is one-to-one by proper 
association of the circular ares with the curves of C), we know (Walsh and Sewell, 


loc. cit., or S, 48 ff.) that, 


Wy (Z) M;, 
<— 


n? 


_< zon S, ton C,, 
w,(t) / 


where VM, is independent of n, z, and ¢, with r > 1 and fixed. 
Thus by virtue of (25) and (27), and the fact that | — 2 2 d>Oforr> 1 by (1), 


we have from (26) 


M : ym 
‘ Ss 
bg = n 


r 


P,,,(z) — @,(2) zon Ss, 


where &,(z) is a polynomial of degree n in 2 
If we let n = 2m, we have 


P.»j(Z) — Pom(z)| < M; (5) : 
r 


and since r/r? = 1/r < 1, we know that for m large 


log m\**2 
f(z) — Pom(z)| S u( = zon Ss. 


m 


Defining p,(z) = Pom(z),n = 2,4, ...,and p, + 1(z) = p,(z), n = 3,5, ., we have 
£7 i 


. log n\** : 
f(z) — pr(z)| SM ,zons, 
n 


for all nm = 2, where VV is independent of n and z. Thus the proof of the theorem is 
, I I 


complete. 





202 MATHEMATICS: H. 8S. VANDIVER Proc. N. A. §. 


* Research supported by the Office of Ordnance Research, United States Army. 

+ The author is indebted to Professors J. L. Walsh, F. J. Dressel, and J. J. Gergen for valuable 
suggestions and assistance during the preparation of this paper. 

1 Sewell, W. E., these PROCEEDINGS, 23, 491-493 (1937). 

2 See, e.g., Walsh, J. L., Interpolation and Approximation, 65 ff. (Providence, 1960), Vol. 20, 
Colloquium Publications of the American Mathematical Society. 

3 Jbid. On pp. 36 ff. the development of the problem in the real and complex domains and the 
results of Weierstrass, Hilb and Szasz, and others, are described. 

‘ Curtiss, J. H., Bull. Amer. Math. Soc., 42, 873-878 (1936). 

5 See, e.g., Sewell, W. E., Degree of Approximation by Polynomials (Princeton: Ann. of Math. 
Studies, vol. 9, 1942), p. 17. Hereafter this book will be referred to as 8, followed by the page 
number. 

6 Warschawski, 8. E., Trans. Amer. Math. Soc., 82, 300-322 (1956). 

7 Walsh, J. L., and Sewell, W. E., Duke Math. Journal, 6, 658-705 (1940). An examination of 
the proof shows that the constant M/, in (12) depends only on the function y(w). Actually (12) 
is established there for several mutually exterior Jordan curves under conditions more general 
than (1) and (2). 

8 Mergelyan, 8. N., Uspehi Mat. Nauk (n.s.), 7, 31-122 (1952). Amer. Math. Soc. Translations 


No. 101 (1954). 


SOME ASPECTS OF THE FERMAT PROBLEM (FIRST PAPER) 
By H.S. VANpDIVER* 
THE UNIVERSITY OF TEXAS 
Communicated Jannary 3, 1961 
We shall consider first! the equation 
v+y'+2' =0, (1) 


where x, y, and z are nonzero rational integers, and / is an odd prime. For certain 
purposes of comparison we shall also consider the equation 


at e6'+y7' = 0, 


where a, 8, and y are nonzero integers in certain algebraic fields. 

Owing to the fact that the Fermat problem has been attacked by many first- 
rate mathematicians, it seems to have acquired practically a unique reputation 
because of the great difficulties encountered by said mathematicians in their efforts 
along this line. However, if we take / > 7, I wish to note that so far all criteria 
which I have observed for the solution of (1), which have led to proofs that (1) 
is impossible for certain values of |, have involved divisibility properties of certain 
rational integers. That being the case, we shall indicate here why the Fermat 
problem should not be regarded as at all unique mathematically. In another 


paper? I have considered a number theoretic problem which, in my opinion, is 


as difficult, if not more so, and it also involves divisibility. It has to do with the 
so-called Wilson quotient, which is 
(p-— 1)!+1 


W(p) = ; (3) 
p 
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It is an integer, since by Wilson’s theorem the numerator is divisible by p where 
pisany prime. The problem just referred to is, 7s there a finite number or an infinite 
number of primes p such that W(p) is divisible’ by p? In the same category with this 
problem is the question as to the number of even perfect numbers being finite or 
infinite, and also, the question‘ as to there being an infinity of primes p such that 


2?-! = 1 (mod p’). (3a) 


To my mind, these examples indicate that any problem which necessarily involves 
the concept of divisibility of rational integers may turn out to be of formidable 
depth and difficulty. Jt might be true that the solution of the Fermat problem neces- 
sarily depends on such questions. 

As to the problem concerning equation (3) mentioned above, for example, the 


only method of procedure to find out anything at all about it, known to me at the 


present time, is to calculate the remainder of W(p) modulo p using a rapid digital 
computer. In this way we pass from one prime to another in our computations. 
Since there is an infinity of primes, obviously this procedure cannot be continued 
indefinitely. Now perhaps the problem is of such a character that this is the only 
possible method of approach. Apparently nothing in logic allows us to conclude 
that there must necessarily exist a finite number of syllogisms leading to a proof either 
that the number of primes p satisfying W(p) = O (mod p) is finite or is infinite’ If 
not, it may be, then, that the Fermat conjecture concerning (1) may be correct 
but never will be proved. So far, of course, in proving® that (1) is impossible for 
any | < 4002, part of the procedure was to examine each | in this range separately, 
using the known criteria for the solution of (1) for / any given prime. In view of 
the last statement, it may turn out that this is the only method which will yield 
any new facts about (1). 

Some mathematicians argue that in view of the difficulties encountered in the 
attempts so far to dispose of it, the Fermat conjecture, one way or the other, is 
not correct; that is, there exist rational integers x, y, and z satisfying (1) for at 
least one value of 1. Personally, I do not worry at all as to whether the Fermat 
statement is true or false. For the last thirty years or so of the years that I have 
worked on the problem, many of the results obtained as by-products of these 
investigations have led me to the publication of a number of articles which do not 
always, on the surface, appear to have anything to do with the Fermat conjecture. 

In this paper I shall give a set of criteria for the solution of equation (1) which 
may be new, and in a second paper to follow I hope to give two other quite different 
sets of criteria, possibly also new. 

The criteria derived in the present paper for Case II of (1) are all elementary, 
involving only rational integers in their statements. On the surface, at least, 
they appear quite complicated, but then each of them has characteristics not 
possessed by either of the others or of any of the known criteria for the solution 
of (1), as far as I can observe. However, in common with a number of other 
criteria which have appeared during the history of the problem, they have not so 
far yielded any proofs of Fermat’s statement for special exponents. Also, in 
conversation with mathematicians interested in the subject of this paper, I find 
that they often have speculated a great deal about the truth or falsity of the 
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Fermat statement. The criteria given in this paper and the one to follow may 
give them more food for thought along this line. 

The possibility that in order to make use of our criteria it might be necessary 
to solve a problem more difficult than the one we started with would not be a new 
situation. 

In what follows, we shall designate by R, the ring of rational integers. We 
shall also designate by R; the ring of integers in the field defined by ¢, where 

c= e7* v 

irst we consider the case where x, y, 2 is prime to / in (1), usually called Case I 

of (1). The writer has given’ the criteria that if (1) is satisfied in this case, then if 


5 sk 
(EE) a1 (4) 
(x) 


fork = 0,1, ..., 1 — 1 where a is a primary integer in R;, such that the principal 
ideal (a) is also prime. A primary integer in R; is defined as being of the form 


(a) is prime to xyz, 


6' + od! (4a) 
where 6 and a € R,, ((6), (oA)) = 1, and A = (1 — §). In (4) the symbol on the 
left is 3° where, if (a) = p, 

(x + yo*)° = ¢* (mod p), 


and where 


N(p) -— 1 
i * 


with N(p) being the norm of a in R,. Consequently, (4) holds if and only if 
x + yt = p,' (mod p), (6) 


where p, « R-. The norm of p in R, is a rational prime which we shall designate 
by p. We shall now prove the 

LemMa. Jf we R;, and m = ayty...0,, the a’s being prime each to each, with 
a}, a, ..., a, belonging to R,, and 


w = a! (mod a)), w = a.! (mod a), ..., w = a,’ (mod a,), 
then we may determine an integer a with a € R, such that® 
w = a' (mod m). 
We may select @’s e R- such that 
B,' = 0 (mod 1m/a;), 
, s, and since a, is prime to m/a,, 
B,' = 1 (mod aq,). 


We shall now show that a value for @ in (7a) is 


ay'By' + a2'B.' + ... + a,'B,'. 
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Using the relation (7) we obtain 
By'w = (Bya,)' (mod m), (11) 


fork = 1,2 s. Then on adding all these expressions we obtain on the left, 
£ I 


> a) “9 


using relation (9) for each k, 


w >> (Bx)! (mod m). (12) 
k=1 


Now the right-hand member is congruent to 


S l 
( > bras) (13) 
k=] 


modulo m. For, if this expression is expanded, every term in the expansion 1s 
divisible by m except those which appear in (6), since in (13) any term involving 
8, and 8;, with 7 # j, is obviously divisible by m, by (8). This proves the Lemma. 

We shall now examine the relation (6) and apply our Lemma in order to generalize 
it. Let in the statement of the Lemma the m be the m = p; given in (6), and suppose 


that we have 
Pp DiPe. . Vs, (14) 
where the p’s are prime ideals in R;. The p’s are known to be distinct in R,, and 
s divides (1 — 1). 
Irom (6) we may write 
ut yé pu! = pu(mod p;) 
where py € R;, and more generally we may write, 
x + yo = piu! (mod p,) 
-. , 8s Now applying the Lemma, where the a’s are replaced by 
, p, and the a’s by the corresponding p’s, we obtain 
x + yo = 6' (mod (p)) (16) 
where 6 ¢ R,, using (14). Whence 
a+ yo = 0 = py, 
where y ¢ R;, and we may therefore write 
Y= + af t+ 
and 


6 = bo + bye + 1—2¢ : (18a) 


where the a's and b’s ¢e Rj. 

We may now expand @! in (17) by the multinomial theorem and use this expan- 
sion in the left-hand member of (17) after applying the relation ¢' = 1. This 
gives an expression of the form 


M+ Me. + M,¢'", 19) 
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where the W's are homogeneous expressions of the /-th degree in the b’s. Further, 
(19) may be put in the form 


(My — Miu) + (Mi, — Min)e +... + (Mi — Midge, (20) 
hy the use of the relation 
ht ie 


Substituting the values of @ and y given by (18), (18a), and (20) in (17), we obtain, 


since 1, ¢, ¢, ..., ¢'-* form an integral basis in R; with respect to Ri, 
M,; — M,.1= 0 (mod p) 
And also, 

My) — M,_, (med p), 
y = M, — M,-4 (mod p). 

In connection with equation (1) we know that 

arty =v', 
o from equation (22) and equation (23) we obtain 
v= My + M, — 2M, (mod p). 


From (21) we obtain also 


hee 1 
> M,= (l — 2)M,-1 (mod p). 


Also, from (18a) and (19), we have 


lt 1-2 \1 
YM, = (34). (25b) 
k=0 i=0 

We now consider the quantities v, from (24), and bo, bi, ..., bi-1 as defined in 
(18a). They can be considered as unknown rational integers, so (21) and (25) consist 
then of a set of conditional congruences in said quantities. They are all homogeneous, 
of degree 1, and 1 — 1 in number. Further, in connection with these congruences 
the prime p must be prime to x, y, and z. It is always obviously possible to find 
such a prime by taking it sufficiently large, as an infinity of primes exist of the 
form’ (4a), when we let o range over ail integers in R-. 

We shall now proceed to obtain other criteria for Case I and also some criteria 
for Case II of the theorem, but the proofs will be indicated. In obtaining (4) 
we employed the x and y in (1). We might just as well use the x and z in (1) and 
obtain another relation of the type (4) where y is replaced by z in that formula. 
Similarly, we can take another pair from the trio x, y, z, namely (y, z), and obtain 
a third relation of this same type. We may proceed with any of said relations as 
we did in connection with (4), and we shall obtain congruences of the type of (22), 
(23), and (24), except that the values x, y on the left may be replaced by any one 
of the other pairs (x, z) and (y, z). Also, of course, in this procedure in the con- 
gruences (21), (22), (23), (25), (25a), and (25b), the M/’s are to be replaced by 
other quantities corresponding to each of the above-mentioned pairs selected from 
x, y, and z. 
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The above-mentioned criteria refer only to the first Case of (1). We may now 
indicate how to obtain a similar set of criteria for Case II of the Fermat problem. 
lor this discussion we take x, y, 2 in (1) as prime each to each with y = 0 (mod J). 
It is then known that it necessarily follows’ that y = 0 (mod /*). Hence by the 
definition of a prime integer in FR; as given in (4a) it follows that 2 + y¢“is primary. 
It is also known that x' = x (mod [*), so x + y¢* can be written in the form x2! + 7d! 
since 1 = 0 (mod \‘~!); then in relation of (4) when applied to Case II we may 


state that if (1) is satisfied in Case IT, 


xt ye 
= |, (26) 


where 7 is any element of R, which is prime to x, y, z. We may then proceed as 
we did in deriving our criteria for the first case and obtain as before a set of 1 — 1 
congruences, modulo p, which are homogeneous and of degree 1, and 1 — 2 of them 
involving 1 — 2 unknowns and the other involving / — 1 unknowns. 

We noted above that in (4) when we consider (1) in Case I, we had similar 
relations involving two other pairs, (x, y) and (y, z). Similarly, from (1), the 
relation of (26) still holds if we replace x by z. If we take instead of the quantity 
(x + yf*) in (26) the quantity (7 + ¢z)/(i — ¢) a primary integer in R;, since y 
is divisible by /, we may proceed similarly as we did in connection with (26); how- 
ever, we cannot use the relation obtainable from (1) since y = 0 (mod J), 


(26a) 


as we do not know the value of m. Hence we will have one less conditional con- 
gruence which we indicate we would find as criteria in connection with (26). 
The criteria of the type of (21) through (25a), inclusive, are certainly very 


complicated, as, for one thing, they involve multinomial coefficients. However, 


I do not see why we should, on this account, necessarily discard the consideration 
of them. For example, in Kummer’s Memoir® of 1857 he gave the criteria (named 
after him) that if (1) is satisfied in nonzero rational integers x, y, z, prime to /, then 


71—92n/.. t 
a pa | - €@y) - 
a - = () (mod l), (27) 
dy?" si 


forn = 1,2, ..., (l — 3)/2 where B, is the nth Bernoulli number, B,; = '/.5, B. = 
Also, x and y in (4) can be replaced by any one of the pairs (2, 2), (y, 2), 
*), (z, 2). As they stand, these relations are quite complicated. Aside 


from the fact that Kummer was able to show from these relations that By,,—3) 2 and 
Bis, are divisible by 1, nothing was done with (27) until forty-eight years later when 
Mirimanoff’ simplified these congruences, and this work led to a long series of 
investigations of Case I of (1), leading to the result that (1) is impossible! for all 
values of 1 S 253, 747, 899 in said case. 

In order to handle our present criteria to obtain new relations, it might be 
necessary to find out something about the possible independence of the congruences 
(21) together with (25). We can immediately note that (25) is independent of any 
of the relations of (21), as the former involves v, which does not appear in any of 
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the other congruences. Hence we can say that at least two independent congruence 
criteria correspond to every prime p as modulus, which is greater than a certain 
limit, and it is also the norm of 2 primary element of R,, provided that | > 3. 

A number of other writers have considered homogeneous simultaneous congruences, 
and there is a well-known theorem'* proved by Chevalley, which states that a system of 
homogeneous polynomials with coefficients in a finite field, the number of whose variables 
is greater than the sum of the degrees of each such polynomial, has a common solution 
with not all of the variables zero. This result is obviously independent of the actual 
values of the coefficients in the polynomials, so since the residue classes modulo p, 


with p a rational prime, form a finite field, perhaps we could do something with our 


criteria along this general line, in spite of the complicated character of our coeffi- 
cients. 

There are several other features about the above mentioned criteria which 
differentiate them from any other set that I know of. First we note that Inkeri'® 
proved that (2) was impossible for a, 8, and y for nonzero elements in the ring of 
integers R, in the field defined by » = ¢ + ¢7', which obviously includes the 
rational ring. However, we cannot, using our methods at least, make any exten- 
sion of our criteria which will make them apply also to (2), with a, 8, and y in R(y). 
For if the M’s in our relation in (20) belonged to R,, then we could not infer our 
congruences in (21) since 


att e'?#+...+1=0 


is not in general irreducible in Ry, as, for example, in the case / = 5. Hence our 
relations in (21) and (25) may constitute more stringent criteria than some of 
those which were previously used® in connection with (1) and a generalization of 
which Inkeri used to obtain the result stated at the beginning of this paragraph. 

It is planned to add to the second paper referred to above some further discus- 
sions of the criteria we have obtained here. 

* The work on this paper was supported by Basic Research Grant 8238 from the National 
Science Foundation. 

1 For a partial history of the problem, cf. Vandiver, ‘‘Fermat’s Last Theorem,’? Amer. Math. 
Monthly, 53, 555-578 (1946). 

2 Vandiver, H. 8., “Divisibility Problems in Number Theory,” Scripta Mathematica, 21, 15 
16 (1955). 

3 There are just three values of p for which (3) is divisible by p. Cf. Karl Goldberg, “A Table 
of Wilson Quotients and the Third Wilson Prime,’ Jour. London Math. Soc., 28, 252-256 (1953). 

4 There are just two values of p < 26,000 satisfying the relation (3a). 

> As otherwise expressed, the Fermat statement may be undecidable. A number of mathemati- 
cians with whom I have discussed this matter regard the situation just mentioned as quite pos- 
sible, and I have noted several discussions in papers on logic which mention the Fermat Theorem as 
probably undecidable. 

6 Lehmer, D. H., Emma Lehmer, and H. 8. Vandiver, ‘‘An Application of High-Speed Comput- 
ing to Fermat’s Last Theorem,”’ these PROCEEDINGS, 40, 25-33 (1954); Selfridge, J. L., C. A. 
Nicol, and H. 8S. Vandiver, ‘Proof of Fermat’s Last Theorem for All Prime Exponents Less Than 
1002,”’ these PROCEEDINGS, 41, 970-973 (1955). 

6 Cf. Vandiver, “Examination of Methods of Attack on the Second Case of Fermat’s Last 
Theorem,’ these PROCEEDINGS, 40, 734 (1954), the paragraph just following the statement of 
Theorem VI. In effect, this says that all the criteria we discussed in that paper break down if 
the second factor of the class number of R¢ is divisible by /. 
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THE POSSIBILITY OF TSCHEBYCHEFF QUADRATURE ON INFINITE 
INTERVALS 
By Hersert 8. WILF 
UNIVERSITY OF ILLINOIS AND ARGONNE NATIONAL LABORATORY 
Communicated by J. L. Doob, December 8, 1960 
Introduction. Let d(x) be a positive measure, of total mass one, on the interval 
(a, b) of the real axis, and let N be a given positive integer. We say that 
Tschebycheff quadrature is possible for NV if there are numbers 2, x2, ..., IN 


lying interior to (a, b) such that 


eb Vv 
} f(x)dY(x) = 1/N » F(x.) (1) 
/a ’ l 


for every polynomial f(x) of degree SN. If Tschebycheff quadrature is possible 
on a sequence }n,}{ of values of N tending to infinity, we will say that the measure 
dy(x) has property 7 on (a, 6), and, in any event, the sequence }n,;{, finite or in- 
finite, of values of N for which Tschebycheff quadrature is possible will be called the 
7-sequence for dy(x). 

Example 1: d(x) = x-(1 — x*)~'*dx has property T on (—1, 1) since, in 
this case, (1) is identical with Gauss-Jacobi quadrature. The 7'-sequence is simply 
n,=jG=1,2,...). 

Example 2: dv(x) = dx does not have property 7 on (0, 1). This was proved 
by 8. Bernstein,! 
5. 6.7. 9}. 

Example 3: For d(x) = [T(a)|~'ve-'e-*de on (0, ©), nothing is known 
except? that the 7-sequence contains V 1,2 and does not contain N for3 SNS 
10, if a = 1. 

Example 4: For d(x) mr '?e-"dx on (— ©, ©), nothing is known except 
that the 7-sequence contains NV 1, 2, 3 and does not contain N for3 << N ¥ 10. 


who later? showed that the 7-sequence in this case is } 1, 2, 3, 4, 


The results of the present paper are first a simple necessary condition for a 
measure to have property 7’, second, an application of this condition to show that 
the 7-sequence of a measure on an infinite interval is, roughly, very ‘“‘sparse,” 
third, an application of the same condition to settle the cases of Examples 3 and 4, 
showing in each case that property 7 is not present and exactly determining the 7- 
sequences involved, and finally, some remarks about Bernstein’s method. 
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It is the author’s conjecture that if a measure has property 7' then it has zero 
mass outside of some finite interval, but as will be seen, our methods are not quite 
strong enough to prove this. 

Jensen’s Inequality.—Jensen’s inequality’ asserts that if &, &, ..., & are non- 
negative, and if we define 


yen 1 


then s<timplies o, 2 o;. Taking f(z) 


b N 
My = f r’dy(x) = 1/N 2, 


v 


Suppose a 2 0. If Tschebycheff quadrature is possible for NV, then zx, 2 


(v = 1, ..., N), and therefore 
or = (Nu,)*!” (yp = 1,2,...., 8%) 
is decreasing. Even if a < 0, however, we can apply Jensen’s inequality to 
j< = 
1 =1> en" 
(y=1 
and deduce that the sequence 
t, = (Nue,) a, ss.5 WD 
decreases. 


TueoreM 1. Let {n,;}° be the T-sequence of a measure with moments u,( 
Uj5 1 Y 4 
0,1, ...). Then for each fixed j(j = 1, 2, ...) the sequence 


Tr = (Nypo,)!” (ry = 1,2, ..., [n,/2)) 


decreases. If a 2 O, then actually, 


or = (nyp,)*" (po 1,2,..5. <5:989) 
decreases. 

We now investigate the density of T-sequences. Let {n,{ 7 be the 7'-sequence 
of a measure dy(x) on (—~, ©). Let {2p,}7, {2¢; + 1}7% be, respectively, the 
subsequences of even and odd integers in {n;}. Then, for each j, the sequence 

tr = (2p jMe,)*!” a + «9 ie 
decreases. Hence, in particular, 
Tpj-1 Pa 
or (2p jhop; )v%-1 = (2p jmop;) 1/77. 
Temporarily writing 4; = pe2p;, we have 
hy S (Qp,j)*P7—-1—-Y j_g?H/ 77-1, 
It is easy to see recursively that 


Ay S j(Aa)"!"1 


where the y, satisfy 
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- 9» (Pp44/P 1) Pr Dp 
Ver1 = (2px+1) "4 1/ "kh Ve RTs 


71 | (i = 
Now the solution of (10) is 


k 
l 
log ye" = Do ( (2p,) 


=) 


k 
| 
= > ( p; + O11) 


j=2 \Pjt 
Suppose 


Then, log y,!/" S 8B, 
and (9) becomes 


On the other hand, 


Mop; i) xidy(x) 
j . X X x } 
1 +f +f f x*idy(x) 
x xX X 
\ 
2 \/ aa f , xidy(x) 


x X 
X x 
ed a Plage 
' x | 


X*°jF(X). 


By i/p Pj 
a 2 ery” 
hus, W(X) < ( — ) 
X? 


and therefore dy() has zero mass outside of the finite interval 


le < ef v1! 2P\ 


teturning to the series (12), we have first 


ot) 
log PA = 
Pj-1 Pp 


and therefore, for any sequence p; we have 


: l ] 
po ( _ ) tog Pjiis 
j=2 \Dj-1 Pj 


Suppose there is an integer m such that 


Pi S Pj 
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) log pja + 


k 
Di l l 
) oe < 0(1) + (m-—1) X = log pj 
Pj-1 j=2 \Pj-1 P; 
= O11) (k > &) 
and the series converges. The example 
Doe ea) (9 2 ca), 
for which (12) is easily seen to diverge, shows that some restriction of the type (17) 
is essential. The same argument can be applied to the odd subsequence of }7n;{ 


with the identical result. We summarize with 


Turorem 2. Let }2p;{ 7, {2q; + 1} 7 be the even and odd subsequences of the T- 
sequence {nj} of a measure d(x) which has positive mass outside of every finite 
interval. Then the assertions 
p (17) 


qj (18) 


are false for every fixed integer m. 
actually 


n; (19) 


is false for every fixed integer m. 
The Classical Measures.—With Theorem 1 we can now easily settle the classical 
cases of the Laguerre and Hermite measures. Indeed, in the Laguerre case, 


; ms (a + 1r) 
ur = |e) ff a’ t*-"e—dr = . ; 
0 l(a) 


hence, by (3), the sequence 


NV(a@ + r) |!” 
O; : Gee (ay A 
I'( a) 


must decrease. But 
On—1 j NT(a + N) ( 
oN (Tia)\(a+N — 1)‘ 
which is less than unity as soon as 


rie) > l(a = N — l) ~ AN%e—" 
N (a + N — 1)" 


(N— o-), 


If a = 1, (20) holds as soon as N 2 3. 
In the Hermite case, the analogue of (20) is 
(2m — 1)"— 


. (2m — 3) = 
2m 
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if N = 2m, which holds for m = 3, or N 2 6. For N odd, we find N 2 7 by the 
same method, and using the result of reference 3, conclude that Tschebycheff 
quadrature is impossible for N 2 4. 

THEOREM 3. The measures of Laguerre and Hermite do not have property T. Lf 
a = | in the former case, the T-sequence ts } }. os. in the latter case, the T-sequence ts 
ri, 3 Bi. 

Remarks on Bernstein’s Method.—lf d(x), (a, b) are given, let }¢,(x)} > be the 
sequence of orthogonal polynomials thereby generated. If £1, < fa. < ie oe 
are the zeros of @,(x) and Aji, Ane, . Ann, the associated Christoffel numbers for 
Gauss-Jacobi quadrature, Bernstein’s criterion,” in the present terminology, is that 


Ani o(n~') it (22) 


implies the absence of property 7. We remark that Bernstein gives the proof for 
d(x) = dx on (0, 1), but his argument is perfectly general. If this result is used in 
conjunction with the estimates of Winston® for \,1, one can show that Hermite’s 
measure does not have property 7’, and that Laguerre’s, for a > | only, does not have 
property 7’. No information is obtained about the exact nature of the 7-sequences 
because of the imprecision of known estimates for \,:. I would therefore regard 
the methods of the preceding sections and those of Bernstein as being complemen- 
tary in that the usefulness of the former is restricted to infinite intervals and that of 
the latter to finite intervals, apparently. 

1 Bernstein, S., “On the Formulas of Approximate Integration of Tschebycheff,’’ Jzvestia Akad. 
Nauk. USSR, Math. and Phys. Sciences Series, 1219-1227 (1932 

? Bernstein, S8., “Sur les formules de quadrature de Cotes et de Tschebycheff,’’ Comptes Rendus 
de L’Académie des Sciences de L’URSS, 14, 323-826 (1937 

Salzer, H., “‘Equally Weighted Quadrature Formulas over Semi-infinite and Infinite Inter- 

vals,’’ Jour. Math. and Phys., 34, 54-63 (1955). 

‘Hardy, Littlewood, and Pélva, /nequalities (Cambridge University Press, 2d ed., 1952 

® Winston, C., “On Mechanical Quadratures Formulae Involving the Classical Orthogonal 
Polynomials,’ Ann. Math., 35 (3), 658-677 (1934) 


A RADIOBIOLOGICAL ANALOGY BETWEEN MEASLES VIRUS AND 
TEMPERATE PHAGES* 
By Paut H. ACKERMANT AND FRANCIS L. BLACK 
YALE UNIVERSITY SCHOOL OF MEDICINE 
Communicated by John R. Paul, December 16, 1960 
Measles is in many ways typical of viral diseases which confer an enduring im- 
munity to the host. A question of long standing is whether etiological agents of 


these diseases remain latent in the body subsequent to infection, continually re- 
leasing antigens into the blood.! This type of latency might be expected if a 


phenomenon similar to the lysogenic relationship of temperate bacteriophage 


existed? and the viral genome were carried, together with cellular genetic informa- 
tion, from one cell to its progeny. While lysogeny has been extensively studied in 
bacteria, there has only been a little evidence to support the hypothesis of an analo- 
gous system in animal cells. 
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Garen and Zinder* found an inverse correlation in bacteriophage between re- 
sistance to ultraviolet light (U.V.) and virulence. In addition, if temperate 
phage P22 was irradiated with U.V. and then plated on normal bacterial hosts, 
virus titers were higher than if the irradiated phage were tested on cells that had 
also been irradiated. Yet, cells irradiated for the same length of time were fully 
able to support growth of unirradiated phage. When the virulent phage, T2, 
was irradiated, identical titers were obtained on irradiated and unirradiated cells. 
It was suggested that the host cells possessed a radiosensitive component which was 
capable of repairing U.V.-damaged P22. 

Similarly, Rubin and Temin‘ found that Rous Sarcoma Virus (RSV) was con- 
siderably more resistant to U.V. than the cytocidal Newcastle Disease Virus 
(NDV) but only slightly more resistant to X-rays. The yield of RSV was lower 
when tested in U.V.-irradiated cells than in normal cells whereas the infectivity of 
NDV was unaffected by irradiation of its host cells. These authors have extended 
the concept of cellular repair of damaged virus to explain their results with Rous 
Sarcoma Virus. 

Powell has reported that no cellular repair of virus was demonstrable in the 
herpes virus system.° 

Experiments presented in this paper were designed to investigate the possibility 
that measles virus inactivated by U.V. might be partially reactivated by mammalian 
host cells. If this were the case, the viruses of measles, Rous Sarcoma and tem- 
perate bacteriophage would all have a basic similarity in their host-virus relation- 


ULTRAVIOLET IRRADIATION ELECTRON IRRADIATION AT 4°C 
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(e) 30 60 120 180 6 X10! 
UV DOSE (SECONDS) IONIZATIONS /cm? 
Fic. 1.—Inactivation curves for measles virus and poliovirus irradiated with ultraviolet 
light and with high-voltage electrons. The electron irradiation data are redrawn from the 
studies of Benesh et al.® 


ship. Virulent Type 1 polio virus was chosen as a cytocidal virus to contrast with 
measles virus, as both could be titrated by means of plaque assay on tissue culture 
of monkey kidney cells. 
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Vaterials and Methods.--Yhe Edmonston strain of measles virus grown in Hep-2 cultures, a 
serial line of human laryngeal epithelioma,® and the Mahoney strain of Type | poliovirus grown 
in primary rhesus kidney cultures were used in all tests. Measles virus preparations used had 
titers between 108 and 10° plaque-forming units per ml and the poliovirus pool had an initial titer of 
107 per ml. Experimental procedures with measles virus were carried out in Erythrocebus patas 
monkey kidney cultures and with poliovirus in rhesus kidney cultures. 

To obtain cell suspensions for irradiation, 6-day-old cultures grown in 3-ounce prescription 
bottles were incubated for 15 minutes at 37°C with 10 ml. of warm 0.02 per cent versene. The 
cells were then centrifuged for 10 min. at 600 rpm and resuspended in Hanks solution with 0.5 
per cent lactalbumin hydrolysate, 2 per cent calf serum solution and antibiotics (HLS). After 
irradiation, the cells were reseeded into the same number of bottles as used originally, that is, 
about 4 X 10° cells per culture. Handled in this way, they reformed complete cell sheets in 48 
hr. Exposure to U.V. for more than 30 seconds damaged some of the cell populations to the 
point where they no longer formed cell sheets adequate for plaque counting. 

For virus titrations, the procedures of Hsiung et al. were used.7)* One tenth ml of virus ino- 
culum was added to each bottle and allowed to adsorb to the cell sheet for one hour at 37°C before 
overlaying. Measles plaques were counted after 10 days incubation and poliovirus plaques after 
3 days. 

Ultraviolet irradiations were performed in darkness with constant agitation, and inoculated 
bottles were not exposed to light for at least 24 hours afterwards to prevent the possibility of 
photoreactivation. A Westinghouse type G30TS8 Sterilamp was used at a distance of 42 cm to give 
an intensity of 1590 ergs/sec/cm? as measured by a Westinghouse SM-600 meter. Cells were 
irradiated in petri dishes filled to a depth of about 1.5 mm. with HLS, and viruses were irradiated 
similarly in Eagles medium plus 10 per cent calf serum. Glassware was coated with a water re- 
pellent silicone preparation (Siliclad Clay-Adams) to decrease cell loss by adhesion to the glass. 
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UV DOSE TO POLIO VIRUS 
Fig. 2.—Inactivation curve for U.V.-irradiated poliovirus titrated on normal 
and U.V.-irradiated cells. 


Gamma irradiation was carried out with virus suspended in the Eagles-calf serum medium at 


0°C. The source was a hollow cobalt-60 cylinder of the Brookhaven type located in and used 
through the courtesy of the Biophysics Department of Yale University. The instrument pro- 
vided 270,000 roentgens per hour 
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Results.—When measles virus and poliovirus were irradiated with U.V., it was 
noted that measles was about twice as sensitive as poliovirus (lig. 1). In contrast, 
to this, a five-fold difference had been obtained previously with high-voltage elec- 
trons, a form of irradiation which causes ionization® (lig. 1). This suggested that 
the radiosensitivity of measles virus might bear similarities to temperate phage 
and that of poliovirus to virulent phage. Further experiments to check this pos- 
sibility were undertaken. 

It was found that Type 1 poliovirus was inactivated exponentially by U.V. over 
a 6-log range. When the rhesus kidney cells were irradiated with U.V. for 15 
or 30 seconds before infection with poliovirus, their capacity to grow poliovirus 
was unaffected, and irradiated and unirradiated cells gave nearly identical plaque 
counts whether or not the virus had also been irradiated (Fig. 2). In this respect, 
poliovirus behaved like lytic bacteriophage.* !° 

When measles virus was irradiated and titrated on unirradiated cells, there was 
again an exponential decrease in titer with increasing U.V. dose. However, if 
virus was similarly irradiated and tested on irradiated cells, there was a decrease in 
titer beyond that seen with normal cells (Fig. 3). This further decrease in titer 
was demonstrable with virus irradiated for as short a time as 5 seconds, and ir- 
radiation of virus for 15 seconds gave nearly maximal effect. 

The figures for the total effect of irradiation of virus for 5, 10, and 15 seconds as 
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Fig. 3.—Inactivation curves for U.V.-irradiated measles virus titrated on normal 
and U.V.-irradiated cells. 
measured on cells irradiated for 60 seconds, the maximum feasible dose to the cells, 
give an inactivation rate about five times that noted with U.V. irradiation of 
poliovirus. This ratio is similar to that found with ionizing radiation. Thus, 
the contribution of the cells to growth of U.V.-damaged measles virus is adequate 
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to explain the difference in sensitivity to the two forms of irradiation. This re- 
lationship does not hold for virus irradiation periods greater than 15 seconds. 
Curves based on titration of virus irradiated for periods greater than 15 seconds 
were straight and parallel to the line of tests on unirradiated cells. Thus, the 
initial slope was dependent on the U.V. dose given the cells, but the slope of the 
later portion of the curve was independent of this factor. 

Increasing the irradiation dose to the mammalian cells had progressively less 
effect on the virus titer. Extrapolating to the ordinate the parallel portions of the 
slopes of the curves for each cell irradiation dose in Figure 3 gives an estimate of 
the percentage of cell units still capable of growing damaged virus. Plotting these 
extrapolated values for cell irradiations of 0, 15, 30, and 60 seconds on a log scale 
(Fig. 4) shows that the curve levels off about 2 logs below the origin; that is, the 
cells may lose up to 99 per cent of their ability to grow damaged virus if irradiated 
with ultraviolet light. 

Figure 5 shows data similar to those in Figure 3 but plotted in a different form. 
Irradiation of cells with U.V. for periods up to 60 seconds produced no significant 
decrease in their capacity to support the growth of unirradiated virus. The titer 
of irradiated virus preparations is progressively less affected by increasing irradia- 
tion of the cells. 

It may be noted that exposure of virus to U.V. for 10 to 15 seconds results in a 
one-log decrease in titer tested on normal cells. The same dose applied to the 
cells reduces the titer one more log, indicating that the sensitivities of the portion 
of the virus essential for growth on intact cells and the portion of the cell essential 
for growth of U.V.-damaged virus are of the same order of magnitude. An ultra- 
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Fig. 4.—Percent of cell capacity to grow U.V.-damaged measles 
virus as a function of U.V. dose to the cells. 
violet dose of 10 seconds to the virus caused 90 per cent inactivation; a similar dose 
to the cells inactivated 90 per cent of the surviving fraction. 

Data presented in Figures 3 and 5 were obtained during six different experi- 
ments, using as many lots of cells and two virus stocks. In one experiment, cell 
irradiation was carried out immediately after trypsinization of the kidneys and 
prior to any growth in culture. The results obtained were similar to those in 
experiments with a later irradiation time. 

In order to determine whether photoreactivation might be a factor in repair of 
U.V.-damaged virus, viruses that had been irradiated with ultraviolet light for 
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15, 30, and 45 seconds were allowed to adsorb onto cell sheets for 105 min after 
which time five ml of Earles with 0.5 per cent lactalbumin and 2 per cent calf serum 
was added to bathe the cells. The bottles were placed on an ice bath and were 
then exposed to illumination by a #2 GE photoflood lamp at a distance of 12 in. 
for 30 min. No photoreactivation was demonstrated when titrations were com- 
pared with a control series of unilluminated viruses. 

When measles virus was irradiated with gamma rays in a cobalt-60 unit, the 
survival curve was the same whether tested on irradiated or unirradiated cells 


(lig. 6). 

Discussion.—The precise nature of U.V. damage to viruses is still largely unknown, 
but it can be reversed in three different situations: (1) U.V.-inactivated particles 
may be reactivated by genetic recombination with undamaged portions of other 


virus particles.'! (2) They may be reactivated in some cases by visible light.' 
(3) They may be reactivated by the action of the host cell in certain systems where 
the virus-host relationship is particularly intimate, as for example, temperate 
phage and RSV. In the present experiments, the multiplicity of infection was 
kept very low in as much as each infectious unit formed a discrete plaque on the 
cell sheet, and hence the first mechanism could not have affected the results. The 
second mechanism was eliminated experimentally. The third mechanism must be 


considered. 
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Fig 5.—Effect of varying the U.V. dose to cells on titer of unirradiated and 
U.V.-irradiated measles virus. 

Garen and Zinder have interpreted their studies on U.V. irradiation of phage 
P22 and its host bacteria to mean that normal host cells possess a genetic component 
that is homologous with part of the virus genome and capable of repairing U.V. 
damaged virus. This mechanism is itself U.V.-sensitive and may be destroyed by 
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irradiation of the host cells. The present data suggest that a similar repair mech- 
anism may be operating in the measles system. However, it would be premature 
to suggest that the mechanism operates through a genetic homology. The effect 
of the cell on measles virus might be on any part of the virus necessary for intra- 
cellular growth. 

The system studied here possesses certain differences from either the P22 or 
the RSV systems. Filtration and irradiation studies® suggest that measles virus 
is larger than either P22 or RSV, and the mammalian host cell is larger than the 
bacterial host. This difference in size was reflectéd in a relatively high sensitivity 
to U.V. Another difference from the P22 system appears in the effect of increasing 
the dose of U.V. to the cells of the test system. Additional irradiation to cells of 
the measles system had progressively less effect on their sensitivity to U.V.-damaged 
virus, whereas irradiation of the bacterial cells had progressively more effect. 
In the measles system, the effect of cell irradiation on virus titer approached a 
maximum at about 99 per cent reduction in virus titer. This limit to the effect of 
U.V. may have been related to the interval of 48 hr between the time of irradiation 
and the time of infection. The interval was necessary to allow the cells to reform 
a sheet suitable for overlay, but it may have allowed time for selection of undamaged 
cells or recovery of damaged cells. Some cell divisions occurred, but the mitotic 
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Fic. 6.—Inactivation curve for cobalt 60-irradiated measles virus titrated on normal 
and U.V.-irradiated cells. 





index during this period was less than 2 per cent, and the total cell count did not 
increase appreciably. 

In spite of these differences, there are striking similarities between the results 
with phage P22 and measles. Both P22 and measles show relative resistance to 
U.V. as compared to other agents tested in the same cells. The capacity of the 
cells to grow undamaged P22 and measles virus is much less sensitive to U.V. than 
the virus itself. Prior irradiation of host cells affected the titer of U.V.-irradiated 
»22 and measles virus but not the titer of gamma-irradiated virus. In both sys- 
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tems, the slope of curves for increasing virus irradiation with constant cell irradia- 
tion became independent. of the dose to the cells when the level of virus irradiation 
was high. The change of the slope of these curves has been interpreted by Garen 
and Zinder® as indicative of a relatively U.V.-insensitive fraction of the total popu- 
lation of repair mechanisms. 

The ability of normal mammalian cells to propagate measles virus that has been 
damaged by ultraviolet suggests a close interaction between host cell and virus. 
This is in contrast to the poliomyelitis, Newcastle disease, and herpes systems. 
By analogy to temperate bacteriophage, the close measles virus-host cell interaction 
might be expected to permit the virus to persist in the body long after initial in- 
fection. If measles persisted in the body in a state comparable to that of prophage, 
it might give a continual release of viral antigens and hence maintain stable anti- 
body levels. This mechanism for the maintainance of immunity to measles would 
provide an explanation of the stability of measles antibodies that has recently been 
noted over long intervals both in the presence and in the absence of exposure to the 
disease. !*. 14 

Summary.—Measles virus was found to be more resistant to ultraviolet (U.V.) 
than to ionizing radiations. The capacity of monkey kidney cells to grow un- 
irradiated measles virus or measles virus irradiated with gamma rays was not affected 
appreciably by U.V. irradiation of the cells. However, the capacity of cells to 


grow U.\V.-irradiated virus was reduced by U.V. irradiation of the cells. In contrast 
U.V. irradiation of cells did not affect their capacity to grow U.V.-treated polio- 
virus. These effects of radiation on the measles virus system are analogous to 
effects on the phage P22 system and are compatible with the hypothesis of cellular 


repair of damaged virus previously suggested for the temperate phage. 


* From a dissertation submitted by P.H.A. in partial fulfillment of the requirements for the 
degree of Doctor of Medicine at Yale. 
+ Aided by a grant from the U.S. Public Health Service, Institute of Allergy and Infectious 
Diseases. 
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INVESTIGATION OF DISEASE RESISTANCE AND THE 
MATHEMATICAL PROBLEM OF PATTERNS 
By RoGer J. WILLIAMS 
CLAYTON FOUNDATION BIOCHEMICAL INSTITUTE AND THE DEPARTMENT OF CHEMISTRY, 
THE UNIVERSITY OF TEXAS, AUSTIN 
Read before the Academy, November 14, 1960 

Immunology, a well-established and highly fruitful discipline, is concerned 
primarily with the investigation of acquired or induced resistance or immunity to 
infectious diseases. We are concerned here primarily with the investigation of 
disease resistances which lie outside this major interest of immunologists. These 
have received relatively little attention. 

What, for example, is the basis for the observation that some individuals are 
immune to an infectious disease from birth to death or that some can be carriers of 
an infectious disease without ever having suffered from it? Still more important to 
our discussion is the apparent fact that each individual in a population has a 
peculiar assortment of susceptibilities and resistances to all types of noninfective 
diseases. 

The fact that people are unequally affected when there is an infectious epidemic 
might possibly be explained on the basis of some previously acquired immunity on 


the part of some of the individuals. But this inequality of response is observed 


also when noninfective diseases are concerned. When pellagra-producing diets are 
consumed in a community only certain individuals are severely affected; when 
polished rice is the basic food of all, only selected individuals exhibit acute beri- 
beri. In regions where endemic goiter is prevalent, only susceptible individuals 
exhibit the diseased condition. Only certain individuals have specific tendencies 
toward gout or arthritis or diabetes, hyperinsulinism, hypo- or hyperthyroidism, 
coronary disease, nephritis, Addison’s disease, duodenal ulcers, or any one of a 
host of other diseases. Other individuals appear to exhibit specific resistance 
toward any one (or more) of these diseases. 

Innate immunity with respect to leukemia, mammary cancer, stomach cancer, 
ete., also appears to exist in some individuals who are free from these diseases 
throughout life, living in the same environment as those who are afflicted. Many 
medical and public health authorities are loath to admit the existence of potent 
genetic factors in cancer because such admission appears to pronounce doom for a 
percentage of the population. Failure to admit does not change the facts, how- 
ever, and failure to face them realistically retards the arrival of the day when 
cancers can be prevented by critical and intelligent alteration of the environment. 

This problem has come to the fore in recent years when it has become reasonably 
evident that heavy smoking is a potent factor in the production of lung cancers, 
that is, of the squamous-cell type. According to the best evidence, however, only 
about one heavy smoker in ten actually suffers from this type of lung cancer. 
Among women, the incidence: is estimated to be somewhat lower. It is 
becoming increasingly evident that some substantial percentage of the population 
is relatively immune to this type of lung cancer; these are not liable to contract it 
even as a result of heavy smoking. Wouldn’t the tobacco companies like to know 
who these invulnerable individuals are and why? 


29] 
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The recognition of innate resistance factors does not conflict with the probability 
that heavy smoking isan important precipitating factor; neither does it confliet with 
the growing probability that virus-like agents are involved in all cancers. It seems 
most likely on the basis of available evidence that only vulnerable individuals 
suffer from lung cancer and that these are most likely to suffer when heavy smoking 
allows the cancer-producing agents (possibly carcinogens and eventually virus-like 
agents) to function. 

Another area in which I am fully convinced that similar resistances and suscepti- 
bilities exist is that of mental disease. If mental disease has a biochemical basis 
(on this point, doubts grow less and less), then it is reasonable to conclude, be- 
cause of the known relationships between biochemistry and genetics, that geneti- 
cally-induced metabolic tendencies play an important part in determining who will 
and who will not be afflicted with mental disease—as well as the type of impair- 
ment.! This does not deny the importance of environmental factors any more than 
recognizing variable susceptibilities to lung cancer denies the importance of heavy 
smoking or the role of virus-like agents in the production of this disease. Again in 
the case of mental diseases, there is a tendency to think that if we admit that 
hereditary factors predispose toward them, this is equivalent to concluding that 
mental disease is inevitable in certain individuals. This is very far from the truth. 
There is a real and substantial hope that recognizing the roots of mental disease—in 
those who are vulnerable—will prepare us to prevent its occurrence. 

Ignoring the innate susceptibilities which exist in every population—or at least 
failing to explore them adequately—has been, in my opinion, a serious block to 
progress in the control of many noninfectious diseases—nutritional, metabolic, 
degenerative—including mental disease, alcoholism, addictions of other kinds, and a 
host of disorders in the mental realm which are not regarded as frank mental disease. 
If the vulnerable individuals in many areas could be spotted in advance, the disorders 
might be prevented with comparative ease. 

In our laboratories, attempts have been made to find biochemical earmarks 
which may indicate vulnerability to alcoholism. Some progress has' been made 
since the last report? but we have encountered a new obstacle in the fact that for 
certain items with which we are concerned there appear to be inter-ethnic dif- 
ferences (highly significant statistically) which probably make some of the pros- 
pective earmarks applicable only to the ethnic group within which they have been 
observed. 

A most crucial lack in connection with the investigation of innate resistances and 
susceptibilities in general is the lack of straightforward and entirely suitable mathe- 
matical techniques for dealing with patterns effectively. 


There exists in each individual because of his distinctive anatomical, physiological, 


neurological, biochemical, and psychological inheritance a host of innate, measurable 
factors.* Individually, each of these items yields, when applied to a population, : 
distribution curve. Collectively, these items constitute a pattern which is dis- 
tinctive and different for each individual and which cannot be handled by tradi- 
tional statistical methods. 

Among the numerous measurable items, some are presumably pertinent to one’s 
resistance to a specific disease, e.g., leukemia, and others would be meaningless in 
this connection but may be highly pertinent to one’s resistance to another disease 
such as arthritis. 
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In the exploration of the factors which enter into the make-up of one’s pattern it 
is essential that we distinguish between those measurements which correspond to 
innate factors and those which are transitory and insignificant from the standpoint of 
inborn characteristics. Though little attention has been directed toward this 
problem, the difficulties it offers are not insurmountable. 

It is also necessary in the exploration of the numerous measurable items to 
select those which are pertinent to any specific disease susceptibility under investiga- 
tion. This selection enters into the problem of patterns. 

Difficulties of mathematical treatment arise not in connection with individual 
items but with patterns. For these there are no averages, norms, standard devia- 
tions, or distribution curves in the usual sense. Various mathematical techniques 
have been used to compare these patterns involving psychological data‘ but none 
fully fills the need. To test out the value of these and other techniques® requires 
that they be applied to the less nebulous data obtainable from anatomical, physio- 
logical, and biochemical measurements. 

I am not, of course, in a position to speak for mathematicians or to keep abreast 
of most recent developments, but it does appear that there is a growing mathemati- 
cal interest in patterns (often for reasons quite independent of our particular 
application) and that we may hope for significant progress, particularly with the 
growing availability of magnificent computers capable of carrying out complicated 
operations with great rapidity. 

The lack of the development of suitable mathematical techniques, readily 
available to biological and medical scientists, has been and continues to be a 
strong deterrent to the collection of necessary data. These data must be extensive 
with respect to each individual studied and are both expensive and laborious to 
collect, and the inclination to collect will not be strong as long as doubt remains as 
to how the data can effectively be treated mathematically once it is collected. 

The potentiality exists for the complete exploration of the innate susceptibilities 
to every disease and to identify in advance those individuals who are vulnerable to 
each. Recognition of vulnerability and the factors which enter into it should go a 
long way in the direction of effective control. 

While it is logical that susceptibility to each disease be studied as a separate 


problem, the techniques for such study are undeveloped, and advancement with 


respect to any disease will throw much light not only on related diseases but upon 
the whole strategy of discovering susceptibilities. 

This largely unexplored problem of disease susceptibility is such an extensive one 
that no individual investigator nor any ordinary team of investigators can do more 
than scratch the surface. It is a problem requiring large and sustained support. 
The Congress of the United States has been and continues to be generous with 
funds to support research in the field of health. It seems obvious to the writer 
that some substantial part of these funds should be channeled into the investigation 
of the disease susceptibility of individual people—about which we know so little. 
We need to remind ourselves that individual people are the only kind that exist. 
Some of the funds devoted to disease susceptibility should be used if necessary to 
stimulate mathematicians to develop whatever tools may be necessary. 

' Williams, R. J., “The Biological Approach to the Study of Personality,”’ (Presented to Berkeley 
Conference on Personality Development in Childhood, University of California, May 5, 1960). 
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EFFECT OF SURFACE TENSION ON THE KELVIN-HELMHOLTZ 
INSTABILITY OF TWO ROTATING FLUIDS* 


By Z. ALTERMANT 
UNIVERSITY OF CHICAGO 
Communicated by S. Chandrasekhar, December 29, 1960 


1. Introduction.—The Kelvin-Helmholtz instability arises when two uniform 
fluids, separated by a horizontal boundary, are in relative horizontal motion. 
Without surface tension this streaming is unstable no matter how small the velocity 
difference between the layers may be. It was shown already by Kelvin! that 
surface tension will suppress the instability if the difference in velocity is sufficiently 
small. 

The effect of uniform rotation about the vertical axis on the development of 
the Kelvin-Helmholtz instability has been treated by Chandrasekhar.?}* He 
found that rotation does not stabilize the motion for any difference in the velocities 
in the absence of surface tension. We shall show that in a rotating system surface 
tension fails to stabilize the motion for any difference of velocities and any angular 
velocity of rotation. The exact effect of surface tension is discussed. 

2. The Characteristic Equation.—Let two uniform fluids of densities p; and p: 
be separated by a horizontal boundary at z = 0.—Let the density, po, of the upper 
fluid be less than the density, p;, of the lower fluid so that in the absence of stream- 
ing, the arrangement is a stable one. We shall further suppose that the two fluids 
are streaming with the constant velocities, VU, and UU», and that the fluids are in 
uniform rotation about the vertical z-axis with an angular velocity 2. Forces 
at the interface z = 0 arise from surface tension 7’. 

As a result of a disturbance, let the components of the velocity in the perturbed 
state be U + u,v, w. Analyzing the disturbance into normal modes, we seek 
solutions whose dependence on x, y, and ¢ is given by 


exp i(kre + kyy + nt). (1) 


Chandrasekhar‘ has shown that the characteristic equation for this arrangement 
is given by 
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pe 
where , : ; ; and 

pi + pe 
and that the Kelvin-Helmholtz instability is least inhibited for perturbations in 
the direction of streaming. Accordingly, let us further consider the case k, k. 


With the notations 
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equation (6) leads to 
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This equation defines a real locus in the (£, 7) plane. The part of the locus in the 
first quadrant connects the points 
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The analysis given by Chandrasekhar* applies to equation (8), and as a result 
there will be stability for 
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By equations (5), (7), and (10) this leads to 
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and after some rearrangements 
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1/2 
+5 Vi ++ Bis) - | = L(k); (18) 
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where 


T 
and a. ——, (14) 
2pi* 
Equation (13) gives the velocity differences for stability of disturbances for any 
given fixed frequency, k. 

3. Long Wave Disturbances.—Yor k > 0 


l l e 1/2 
L(k) > aR + 3(1 = 9 (Ak + Bk)) , — | 
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< z k. (15) 
This result is independent of surface tension and can be stated as follows: For 
any given velocity difference Uy — Us and angular velocity Q, there are perturbations 
of such low frequency as to cause instability. Surface tension has no stabilizing effect 
for such long waves. 
The same result can be obtained from an examination of the properties of the 
(—, n) locus. For a fixed Q, k ~ 0 implies w ~ ©, and by equations (10) 


- 
go > Ww 


Nn w® 


forw—> ©, (16) 


The end points (9) of the (&, 7) locus will coincide, and the locus reduces to four 
points, one in each quadrant; this leads to a real value of n (see equation (3)) 
for stability, for VU; — Uz = Oonly. There is no arrangement with U; — U2. ¥ 0 
and 2 ~ 0 that is unconditionally stable. 


1. Short wave perturbations. York —~ © and Q # 0, equation (13) leads to 
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U,; — U2| < 5 V 2Bk. (18) 


In the presence of surface tension, the coefficient B is different from zero, and by 
(18) the arrangement will be stable for short wavelengths. Without surface 
tension, B = 0. In this case, for any U,; — U>. and Q there will be short waves 
which cause instability. In other words, in the absence of surface tension both 
very long and very short wave perturbations cause instability. Surface tension 
stabilizes the arrangemeat only for short wave perturbations, as is to be expected. 

5. Concluding Remarks.—Surtace tension stabilizes the arrangement. for short 
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wave disturbances. Tor large wavelengths, when the system is rotating there is 
instability both with and without surface tension. It seems surprising that a 


rotating system is not stabilized by surface tension, whereas the nonrotating ar- 
rangement is stabilized. The reason for this becomes clear if we note that the 
nonrotating arrangement is unstable for only short waves, and surface tension is 
able to stabilize these instabilities. In a rotating system there are further in- 
stabilities associated with long waves and these the surface tension cannot prevent. 
It might be noted that the effect of rotation is different for various problems. 
As shown by Chandrasekhar,’ in problems of Rayleigh-Taylor instability in two 
uniform media, rotation does not change the range of stability at all. However, 
in a medium of exponentially varying density, the Rayleigh-Taylor instability is 
inhibited by rotation for long wave disturbances. In the case of Kelvin-Helmholtz 
instability in two homogeneous media, rotation has the opposite effect—it causes 
instability of long wave disturbances. This instability is not inhibited by surface 
tension. 


‘ 


In conclusion I wish to express my indebtedness to Professor 8. Chandrasekhar 
who drew my attention to this problem, for discussions and much valuable advice. 

* The research reported in this paper has in part been supported by the Geophysics Research 
Directorate of the Air Foree Cambridge Research Center, Air Research and Development Com- 
mand, under Contract AF 19(604)-2046 with the University of Chicago. 

t On leave from the Department of Applied Mathematics, Weizmann Institute, Rehovot, Israel. 
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‘ Chandrasekhar, 8., loc. cit., eq. (153 

5 Chandrasekhar, S., loc. cit., chap. x, $95. 


ON A MANY-DIMENSIONAL BWK APPROXIMATION* 
By Don H. Secrest AND JosepH O. HirscHFELDER 
UNIVERSITY OF WISCONSIN THEORETICAL CHEMISTRY LABORATORY 
Communicated November 25, 1960 


The possibility of generalizing the BWI approximation to nonseparable many- 
dimensional problems is discussed. There are two types of difficulty encountered. 
First, the solution of the nonseparable many-dimensional classical Hamilton- 
Jacobi equation may be even more difficult than the solution of the original 
Schrédinger equation. Second, a procedure must. be developed for making the 
necessary connections between the BWK solutions in the classical and nonclassical 
regions. The usual BWK connection formula would be very difficult to generalize 
or to use in a many-dimensional problem. We examine three somewhat. simpler 
types of connection formulas which might be more suited for this purpose. The 
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energy calculated by these alternative procedures differs from the standard BWK 
value only in having a slightly different zero-point energy. 

The procedure developed by Brillouin, Wentzel, and Kramers'!~* (BWK) is 
one of the most useful and interesting means for obtaining an approximate solution 
to the Schrédinger equation for either one-dimensional problems or problems 
which are separable. In the BWK treatment, the logarithm of the wave function 
is expanded into a series in powers of Planck’s constant, 


; 
ln ¥ = z >. ayn. (1) 


Here the yo is the classical action variable, S, which is the solution of the corre- 
sponding Hamilton-Jacobi equation. The yi, ye, ... are other functions which 
determine the higher orders of the BWK approximation. For the first BWK 
approximation, the series is terininated with n = 0; for the second BWK approxi- 
mation, the series includes n 1. For most functional forms for the potential 
energy, the series of equation (1) is only asymptotically convergent.  Further- 
more, equation (1) is not usable in the vicinity of the classical turning points. 
The applicability of the BWK formulation depends upon establishing a connection 
between the approximate wave function in the classical and nonclassical regions 
of configuration space. 

lor many-dimensional nonseparable quantum mechanical problems, it should 
still be possible to approximate the logarithm of the wave-function in terms of a 
power series having the same form as equation (1). This would comprise a generali- 
zation of the BWK method. For the many-dimensional case, the yo is again the 
action S, which is a solution to the Hamilton-Jacobi equation, and the yj, ye, 
satisfy the following relations: 
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The principal difficulty in applying such a generalization is the establishment of 
suitable relations for connecting the solution in the classical region with the solution 
in the nonclassical regions of configuration space. Once the BWK is generalized 
to many dimensions, there still remains the practical problem of solving the classical 
Hamilton-Jacobi equation which forms the zeroth approximation for any such 
treatment. Unfortunately, the Hamilton-Jacobi equation is nonlinear, and for 
most nonseparable many-dimensional cases it is much more difficult to solve than 
the Schrédinger equation. Indeed, Brillouin‘ discusses the problem of many- 
dimensional Hamilton-Jacobi equations and concludes that only for special cases 
can an analytic solution be found. Thus a many-dimensional BWK treatment 
may have limited utility. 
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For most problems, the BWK gives only approximations to the correct energy 
levels. In these cases, different types of procedures for connecting the classical 
and nonclassical zones lead to somewhat different values for the energy. In this 
paper, we examine three modifications of the BWK connection formulas, each of 
which is relatively simple and might be used in a many-dimensional problem. 
These modifications, although much simpler than the original BWK, are not quite 
as accurate. 

1. The One-Dimensional BWK Method.—Virst, let us consider the usual one- 
dimensional BWK method. The one-dimensional Schrédinger equation has 
the form 

h? d(x) ; 

+ V(xr)Wix) EW (x). 

Qu dx? 
If we assume that W has the form given by equation (1) and substitute it into 
equation (5), then the y) and y; are given by equations (2) and (3). For the one- 
dimensional case, equation (2) can be integrated to give 

S(z;, x) = SZ [2u(# — V)| ‘de. (6) 

Here, x; is related to the constant of integration. In the nonclassical region this 


becomes complex, so we define 


§; = S(x;, x) Sc (2u(V — E)]'"dxr = iS,. ( 
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) 
Here, S$; is real in the nonclassical zone. Thus, in the classical zone we may write 
the first BWK approximation as 
V A, sin (S,/”4 + @;). (8) 
The constant 6; is to be determined by boundary conditions. In the nonclassical 
zone, the first BWK approximation is 
V Be3i/*®# 4+ Cea 8i/%. (9) 
The constants B; and C,; are to be determined by boundary and normalization 
conditions. Equation (3) can be integrated to give 
: dy ‘te ; 
Wy '/,2a In ot in |2u(h — V) Lc (10) 
dx 
where ¢ is an integration constant. The second BWK approximations to the 
wave function in the classical and nonclassical zones are therefore 
V, = A,[2u(B — V)]~* sin (S,/A + 8) (11) 
and [(2u(V — E)]~'/*(ByeSi/*# + Cye—5i/4), (12) 


These solutions are good approximations to the solution of equation (5) in both 
the classical and the nonclassical zones, well away from the turning points. Let 


us assume that there are two turning points, | and 2, and that the classical region, 


region Ii, is connected. The nonclassical region is separated into two disjoint 
subregions, | and IIT. In Figure 1, # and V are plotted against x in the neighbor- 
hood of the turning points. In the usual BWK treatment, the potential is ex- 
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panded in a power series in x about the turning point and only the linear term is 
carried. Thus, 


V = EF — aa — a), (13) 


where x; is the position of turning point 1 and —q, is the slope of the potential at 
that. point. For this linearized potential, 


3/o 


S, = f2 [2u(h — V)]'dx = 2/3(2Quay)'“"(x — x1) (14) 
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The Schrédinger equation (5), is solved exactly for this potential, and the only 
solution which behaves properly as « approaches minus infinity is a linear combina- 
tion of Bessel functions of orders '/; and —!/3, known as the Airy function, Ai(z) :* ® 


VW, = AAil(Qua)'?(a1 — x) /h”). (16) 
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and in the classical region, x > 2, 
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Using equations (13), (14), and (15) to eliminate 2, for « < x, 
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and for x > 21, 
) A,(S,/h) ?[2u(E — V)\~ [J 1,,081/%) + J1),(81/%) J. 


lor large x, in the classical region equation (20) becomes asymptotically 
S te) | . . 
Y, = A,[2yu(L —_ V)] “* sin GS; h-+t! 47). (21) 


Equation (21) and equation (11) are apparently identical if A; and 4, are properly 
chosen. Actually, they are not identical, since S; in equation (11) is computed 
from the real potential while the S; in equation (21) is computed from the linearized 
potential. We may use the asymptotic Airy function of equation (21) for con- 
necting the wave function in the classical zone at a point far enough from the 
turning point to use the asymptotic form of equation (16) but close enough to use 
the linearized potential, equation (13). Thus, we must set 0; = 2/4. 
In Figure 2, we have graphed W, for the potential 
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(22) 


with energy EF = «/2. The solid line represents the true solution for this potential, 
the dashed line represents the solution for the linearized potential, equation (13), 
and the broken line is the BWK approximation. 

A similar situation holds at turning point 2. Here, we require S, S(xe, 2) 
as defined by equation (6). Solutions of the form of equations (11) and (12) are 
found. Thus, we obtain for the BWK approximation in region II at turning 


point 2 


W. = As[2u(E — V)]~' sin (S./# — 1/4n). (23) 


The functions given by equations (21) and (23) must agree at every point in 
region II. Thus, A» +A, Furthermore, the difference in the arguments of 


the sine functions must be an integral multiple of x: 
S1 A+! iT — So h +- l :7 L (24) 


Using equation (6) to eliminate S; and So, we obtain 
5 &Y 


2S? (Qu(B — V))]'"de = (n — 1/2)h. 25) 


This integral is used to determine the nth quantum state in the BWK approxima- 
tion. 

In the present. paper, we discuss three alternative methods of connecting the 
solution in the classical zone to the solution in the nonclassical zone. The different 
types of connection formulas lead to different values of @; and to energy relations 
of the same form as equation (25), except that the constant '/» on the right side 
is replaced by a slightly different number. The three alternative procedures 
considered in sections 2, 3, and 4 respectively are: 

(1) The inflection-inflection 1 st-order connection: Here, the first BWK approxi- 
mations, equations (8) and (9), are connected with a straight line between their 
inflection points nearest the turning point. The parameters A, B or C, and 6 
are adjusted so that the straight line connects with a continuous first derivative. 
The slope of the straight connecting line is less than the slope of the solution 
equation (16), obtained by use of the linearized potential. To correct this defect, 
we consider the following procedures. 

(2) The Atry-inflection I st-order connection: Again the first BWK approxima- 
tions in a classical and a nonclassical zone are connected with a straight line be- 
tween their inflection points nearest to the turning point, but the straight line is 
required to have the same slope as the Airy function, equation (16), at the turning 
point. The parameters are adjusted so that the first derivative is continuous in 
the classical region but, as we show, it is not possible to connect with continuous 
first derivative in the nonclassical region. To correct this defect we consider the 
following procedure. 

(3) The Airy-inflection 2nd-order connection: This formulation is the same 
as for the Airy-inflection Ist-order treatment except that now the second BWK 
approximation is used in the classical region. In the nonclassical region, however, 
there is no inflection point, so that the straight line is connected to the second 
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Turning Turnin 
point’ point 


XO 


This figure shows an arbitrary potential function, the two turning points, and 
the three zones. 


BWK function at the point where the slopes agree. This connection has the 


slope of the Airy solution at the turning point and is smooth in both classical and 
nonclassical regions. 

It is hoped that at least one of these connecting formulas will prove more easily 
generalizable to many dimensions than the usual BWK method. 


TURNING 
POINT 


“1.0 -0.5 -0 


Fic. 2.—The usual BWK solution, the true sities to the Schrédinger equation, and the 
Airy function for the linearized pote ntial calculated for the potential V = «z?/2 for x < 0 
and V = 0 for « > 0, with FE = «/2. 
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2. The Inflection-Inflection 1st-Order Connection.—The first BWIX approxima- 
tion to the wave function W; is given by equations (8) and (9). In region I, it is 
necessary that C; 0 in equation (9) in order that the wave function behave 
properly near z = —©. In region III, we take B, = O in order that the wave 
function behave properly near x = +°. The BWK approximation is plotted in 
Figure 3 for the potential given by equation (22). 

The approximate BWK wave function behaves improperly near the turning 
point, although it is a good approximation in the classical and nonclassical regions 
away from the turning points. The simplest way of connecting the BWK function 


TURNING 
POINT 


-0.5 
4 


\/ 
x4) — 


Fic. 3.—The inflection-inflection Ist-order connection calculated for the potential V 
for x O and V 0 for r > 0, with BE «/2. 


in the classical to the BWK in the nonclassical region is with a straight line joining 
the inflection points of the functions in the two regions. The constants A,, B,, 
and 6; are adjusted so as to make the slopes at the inflection points agree with the 
slope of the line. 

The inflection point of ¥; in region I may be found by setting the second deriva- 
tive to zero and solving for x,, the position of the inflection point at a in Figure 3. 
The first. derivative is 


lv [2u(V — E)]” 
= B, ts : € 


dx h 


Si/%, (26) 


We may expand the potential V in the neighborhood of the turning point in a 
power series about the turning point and carry only the linear term in x, as in 
equation (13). Using the expansion of V in equation (26), the second derivative 
of the BWK function in the neighborhood of the turning point is 

dy | —2ya, i = x1)) 2uay(x —_ Fi) 


— Bf, — @S/h _ B, e5i/f, 
dx? 2h h? 
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At the inflection point .,, the second derivative of the BWK is zero and 
(a, — ay)? = (Spa) ~'”. (28) 


Here, we have assumed that the inflection point in the BWK function occurs 
close enough to the turning point for the linearization of the potential, as given 
by equation (13), to be valid. With this assumption, at x, we obtain 


> 7 a . 1 2 9 « 1 2/ \4/ / . 
S21, X4) = jks | —2yuay(x —.x)|’dr = —? 3(2uay) “*(a, — a)” = —1/h. (29) 


Thus, at the inflection point, 
Wi (2,4) = Bie’ f (30) 


na sks oats 
and = | = B,(pa;/h?) “e~ ”. (31) 
av Ta 


In the classical region near turning point 1, the first derivative of the BWK is 


dv, { [2uay(a _— a) )" 
= Ay i 


cos (S,/h + 6), (32) 
dx h 


and the second derivative is 


ay, { [2uay(x = x1) | 
| ee 


sin (S;/h + 6) + 
dx? h? 


9 ..\7'/2 
[2ua;/(2 — 2) ] , a 
A, - cos (S,/A + 6). (33) 
2h 
We assume that equations (13), (32), and (33) are valid at 2, the inflection point 
nearest to turning point 1. At 2x,, the second derivative of the BWK given by 
equation (33) is also zero. Thus, 
: : [S )/h + 0) 34) 
1) i a - = tan [S(2%, 2) /2 i}. ( 
(Suay) “*(x%» — %) ”’ 3S(21, Xp) , 
Setting the first derivative in region I at x, equal to the first derivative in region IT 


at Zp, 


ane [2uai(a, — x)]” | 
B, (2) e "= A, : MAy(A v1) | cos [S(x1, 2p) /# + 4]. (35) 


h? h 
Connecting the inflection points 2, and x, by a straight line whose slope is given 
by equation (31), we find the value of the BWIK< function in region IT at 2, is given by 


Wi (2p) = Bye~' ; + By(uay h?)' "(Lp = Xa)e OP me A\ sin [S(x1, Xp) h + 6}. (36) 
Dividing equation (36) by equation (35) we obtain 
[h2/(uay)]? + (a — ta) = A[Quar(r, — a1)]~'7* tan [S(a1, a)/k + 6:). (37) 


Or, making use of equation (28), 


3[/3S(a1, 2)\" | Slar, a) a 
+ a 9¢ < 
9 h h dN(21, Xp) 


Solving equation (38) for S(a, 2), we obtain 


1 


S(x1, X) = 1/3(38 7% — 1)", 
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The value of 6 is determined from equations (34) and (39), and the ratio of A; to 
B, is determined from equations (35) or (36) and (39) and the value of 4. 

Connection is most easily made at turning point 2 by defining ¥. by equations 
(8) and (9) and then requiring that ¥. = W, at every point. In a manner similar 
to that. used at turning point 1, if 2g is the position of the inflection point of YW. in 
region III, we obtain 


 . 


(ta — Xo)? = h(8yas) i (40) 


and Xe, Ya) = =. (41) 


” 


In a manner similar to that used in arriving at equation (39), we obtain in 
region IT near turning point 2 


S(22, 2) 1/,(3 “7 — 1) 7%, (42) 


where x, is the location of the inflection point in region IT nearest turning point 2. 


From equations (39) and (42) it follows that to within a modulus of 27, 


6, = —b (43) 
The ratio of B, to Az may also be determined. Now we require that at every point 


V2 A» sin CS» h + 6) = A, sin CS; h + 61) = WV). (44) 


If equation (44) is to hold, then we must have 


Ay = + Ao, (45) 


and S:/h + 0, — So/h — & = 


Using equations (6) and (438), this condition becomes 


rom equation (34) and (39), the value of @ is 
6, = —'/,(37 — 1)” + tan 
Substituting this value into equation (47), 
257, (2u(E — V)) “de (n — 0.51085646)h. (49) 


This is very close to the result of the usual BWK approximation, equation (25). 
Kquation (49) gives slightly lower energy levels than does the usual BWK formula- 
tion. , 

3. The Atry-Inflection 1st-Order Connection.—In the inflection-inflection 1Ist- 
order procedure, the inflection points x, and x, of the first order BWK function in 
the nonclassical and classical zones are joined by a straight line, and it is required 
that the slope of the connecting line agree with the slope of the BWK function at 
xv, and 2. The Airy-inflection Ist-order procedure is similar except that now the 
condition that the slope of the BWK function at x, agree with the slope of the 
connecting line is replaced by the requirement that the connecting line have the 
slope of the Airy function which is the solution to the Schrédinger equation at the 
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TURNING 
POINT 


“08 v4 
— ——— 


Fic. 4.—The Airy-inflection Ist-order connection caleulated for the potential V 
for z < Oand V = Ofor z > 0, with E = «/2. 


turning point. This leads to a discontinuity in the slope of the wave function at 
r,. However, wave functions with discontinuous first derivatives are acceptable 
as approximate functions if certain precautions are taken in their use.&7 For 
example, if the energy is determined from the expectation value of the Hamiltonian, 
a correction must be added to the kinetic energy. However, since we obtain the 
energy from the connection in the classical zone, no correction is required to com- 
pensate for the discontinuity. The Airy-inflection Ist-order solution for the 
potential of equation (22) is shown in Figure 4. 

In the neighborhood of the turning point, we may expand the potential in a 


TURNING 
POINT 


fe}! i 
“15 -1.0 


-0.5 
1/4 
xB) — 


Fic. 5.—The Airy-inflection 2nd-order connection calculated for the potential V = x2x?/2 
for r < Oand V = Oforz > 0, with EF = «/2. 


2/9 
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power series and carry only one term as in equation (13). Using this potential 
in equation (5), we obtain the Airy function equation (16) as a solution. Expand- 
ing this solution about the turning point and carrying only the first two terms we 
obtain® 


A,'[l + K (ua, h?) "(xr — X1)}, (50) 


’ I'(?/3) 2 ‘, 
where K : (51) 
r(*/s) \9 


The slope of this solution at the turning point is 


dx, 
dx 


Ay/K (ua, h*) . (52) 


We define ,; in zone II by equation (8). The slope is given by equation (32), 
and the condition on the inflection point is given by equation (34). Connecting 
the BWK function, equation (8), to the straight line, equation (50), gives 

A, sin PS(24, Ly) h + A A,’ l + K (ua; h?) 
Requiring that the slope of ¥; at the inflection point agree with the slope of the line, 


equation (52), at the turning point, we obtain 


[2uay (2p — IZ) ; : rs ei he 7 
A ; cos [S(x1, a) /A + 8 A, K (ua, h?)?* (54) 
U 


Mliminating A,, Ay’, and tan [S(2, 2,)/% + 6] between equations (34), (53), and 
(54), it follows that 


: h h : 
K : ~— : — '/,K = 0. 
d3S(21, Xn) 39(21, 2p) 


The only real solution is 


h jl n | E 
3S(21, Xp) 2 K L4 


ee soe 
(7 - =] | f ; (ob) 


If at turning point 2 we define a new function , and proceed in a similar manner, 
we obtain at inflection point c 


h h a 
Reap — = ot) 
3S (oe, X-) 39(21, Xp) 
We must impose the further condition that ¥, and W. must be numerically equal 
at any point in region II. The conditions for this to be true are given in equations 
(45) and (46). Again, it is easy to show that 6; = —@ and the energy must be 
such as to satisfy equation (47). Solving equation (34) for @ gives 


S(21, Xp) h 
,=—- + tan! 
h 3S(21, Xp) 
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Thus, for a smooth connection with a straight line, which has the slope of the 
solution to the Schrédinger equation with a linearized potential at the turning point, 
we are led to the quantization rule 


25? [(2u(B — V)]'“de = (n — 0.43816371)h. (59) 


Equation (59) gives slightly higher energy levels than the usual BWK approxi- 
mation does. 

It is interesting to note that if, in the nonclassical zone, we had attempted to 
connect at any point where the slope of the BWK function agrees with the slope 
of the Airy straight line, we arrive at an equation with no solution in the nonclassical 
region. This indicates that in the nonclassical region the slope of the first BWK 
approximation never becomes as steep as that of the Airy straight line. The 
difficulty of making a smooth connection in the nonclassical region may be shown 
in the following manner. Let x, be the hypothetical point at which the slope of 
the BWK approximation is the same as the slope of the Airy function at the turning 
point. Then, if we assume that x, is near enough to the turning point 1 that we 
may linearize the potential, the slope of ¥; is given by equation (26). This must 
be equal to the slope of the Airy function, equation (52). Thus, 


— F))” 
h 


[2u(V 
B, mM 


Y 


eX(t1, Ly)/h = A,/K(ua,/h?)'” 


Since ¥; must be equal to x; at x = x,, we have 
l id 9\'/3 } 
Bye(41. Ly) hk — Aj’[] + K (ua, h?) \t, — U1) |}. 


Elimination of B, and A,’ between equations (60) and (61) leads to 


1 /—38(x, z,)\” 1 /—38(x,, x.) : 
= - _ - —~ (62) 
K h 2 h 


a , a ae oer 
Phis is a cubic equation in (—38(x, wx,)/%). It has only one real root, and that 


lies in region II. Thus, we are unable to connect smoothly in zone I. This 
difficulty can be overcome by using the second order BWK function. 

1. The Airy-Inflection 2nd-Order Connection.—The second order BWK approxi- 
mations to the wave function in the classical and nonclassical zones are given by 
equations (11) and (12) respectively. In order that the wave function behave 
properly as x approaches + ©, it is necessary that C, of equation (12) be zero in 
region I and that B, be zero in region III. The second order BWK functions 
become infinite at the turning point. They are plotted in Figure 5 for the potential 
given by equation (22). 

In the classical region, we connect the inflection point of the second BWK 
appreximation, labeled b in Figure 5, to the Airy straight line given by equation 
(50). Since the second BWK function has no inflection point in the nonclassical 
region, we connect it to the Airy straight line at the particular point y at which 
the slope of the straight line is equal to the slope of the BWK function. 

First, let us consider the nonclassical region. In the neighborhood of the turning 
point 1, we may linearize the potential (as in equation (13)) and the slope of % 
becomes 
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dv, ' , [2uay(2y — ZT) : : 
= B, : s(2ua)) (ry oo 2. oa >! hi (63) 
dx h 


The slope of the solution to the Schrédinger equation with the linearized potential 


is given by equation (52). Equating this to the slope of ¥; at some point y and 
expressing the results in terms of S(2;, ,), we obtain 


—3S(x1, r,)\“' —h : 
ba ( = Hs ) [ .( ) 4 | | sens he 4)'K. (64) 
pa,h 3S8(11, Ly) 


fequiring that the Airy straight line have the same value as the second BWK 
function at the point y, 


K —358(214, r.) ; - 
Ay’|1—- ; = B,[—3yua,S(x, x,) © Che FTF (60) 
2 h Y 


Eliminating A,’ and B, between equations (64) and (65), we obtain 


. = K = ; (66) 
2 \38(x1, 2,) —h Z = 


This is a fifth order equation in (—3S8(x, x,)/#)” with only one positive real root. 
Thus the connecting point, y, in zone I is unique. 

In the classical region, in the neighborhood of the turning points, the first and 
second derivatives of the second BWK functions are 
dv, y, [cos (S,/h + 6) sin (S,/f + *) 


A,(3yayS)) (67) 
dx L 
i hn a 


dx? 


h 6S; 


° ’ v0 9/9 ’ 3 » ‘ . ) vs) 
and A, sin (S,/h + 6,) k (ud))"(5uap;) " = (Spay) “ i (O23) 


At the inflection point x,, the second derivative of ¥,; is zero. Setting the second 
derivative equal to zero and solving equation (68), we obtain 
S(a1, %) = 5°°R/6. (69) 


In order to connect the Airy straight line, equation (50), with the second BWK 
function at the inflection point, it is necessary that 21(.) Wi(z,). Thus, 


5 5\'" K o 
A,|3uayS(x1, 2p) * sin ( + a) A,’ E + ( ) | (70) 
6 } 2 


l‘urthermore, at x» the slopes of x; and VW; must match. Therefore, 


[BuayS(x4, Xp) . ‘ - 4 5 
t, — +6) —5 sm | — + 4 
6 6 


I 


and (71), we obtain 
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In the neighborhood of turning point 2, we arrive at similar results and obtain, 
in a manner described in Section 2, the quantization rule given by equation (47). 
Solving equation (72), 6: = 0.39346190 and equation (47) becomes 


252. (Qu(k — V)}'“dxe = (n — 0.25048563)h. (73) 


Kquation (73) gives considerably higher energy levels than does the usual BWK 


treatment. 

Thus, in addition to the usual BWK quantization rule, equation (25), we obtain 
three slightly different rules depending upon the method of establishing the con- 
nection between the classical and the nonclassical regions. Which rule gives the 
best values for the energy levels depends upon the particular problem under 
consideration. 


* This research is being reported under Contract AT(11-1)-298 of the U.S. Atomic Energy 
Commission. 
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THE EFFECTS OF REDUCED COMPETITION IN BIRDS* 
By KENNETH CROWELL 
UNIVERSITY OF PENNSYLVANIA 
Communicated by G. Evelyn Hutchinson, January 3, 1961 


If one were able to remove a majority of species from a fauna, those species re- 
maining might be expected to increase their numbers through the exploitation of 
food and habitats relinquished by the missing species. From the extent to which 
this necessitated the acquisition of new behavior, the original relationships of the 
niches could in part be inferred. The avifauna of Bermuda, impoverished by 
virtue of its isolation, but including several North American species, presents an 
approximation to such a situation. 

Censuses showed the ten common resident land birds of Bermuda to achieve 
total populations at least as great as those of all species in comparable continental 
communities. Eighty per cent of such populations were composed of three species 
found in eastern North America—the cardinal (Richmondena cardinalis), catbird 
(Dumetella carolinensis), and white-eyed vireo (Vireo griseus). The method 
of feeding, feeding height, and nesting habitat were selected as important factors 
in which competition could occur and were compared in Bermuda and North 
America for these three species. 
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If the fitness of a species is plotted along some environmental variable, a dis- 
tribution curve is obtained with the base corresponding to the range of conditions 
under which survival is possible, and the peak overlying the zone of maximum 
efficiency.':? Population density, apportionment of time, and length of life are 
convenient estimates of fitness. In this study, such distribution curves were es- 
tablished for feeding methods and habitat selection. Curves of species occupying 
consecutive zones along a limiting parameter are usually thought to be adjacent but 
not overlapping. Evidence presented here indicates that there is considerable over- 
lap in fundamental abilities. 

Feeding methods were classified as to activity or location, and the time spent at 
each method was recorded. The distributions observed are shown in Figure 1. 
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Fic. 1.—A comparison of the percentage distribution of feeding methods in three passerine 
species common to Bermuda and North America. Numbers show observed time in 10-second 
intervals. Ground-—feeding on the ground without scratching. Foliage—feeding in leafy 
parts of vegetation. Bark—feeding on woody parts of vegetation. Hover—taking of insects 
from vegetation while in flight. Hawk—taking of flving insects while in flight. 


Differences in distributions for the two locations were significant in each species 
(Chi Square: p < 0.01). In the catbird and cardinal this was due to the pre- 
dominance of ground feeding in Bermuda. Diversity H (calculated by }>, pilneDi, 
where p; is the per cent of the total found in the 7th interval)*® was greater in Ber- 
muda only in the vireo. All methods of feeding observed in Bermuda in the re- 
spective species were also observed on the mainland. 

Height zones appropriate to the distribution of foliage were chosen, and the 
relative amount of foliage was measured in each.* Percentage feeding time was 
calculated for each zone and divided by the percent of total leaf area in that zone 
since the distribution of foliage might bias feeding heights. It is seen in Figure 2 
that feeding is restricted to slightly lower levels in Bermuda and that the diversity 
of per cent feeding per zone is again less in the catbird and cardinal and only slightly 
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greater in the vireo. Food of missing species is therefore utilized by those species 
present through application of innate feeding patterns. 

Habitats occupied in Bermuda were among those selected on the mainland, al- 
though drier scrub and absence of forests have produced changes in preference. 
There has been movement by all three species into sparser habitats through nesting 
in low brush and feeding in open areas.* 

In a diverse fauna, territories of a given species are distributed in selected sites 
and do not usually occupy the entire area. In Bermuda there has been little de- 
crease in territory size, and territories usually lie adjacent to others belonging to 
the same species. Increase in population density has therefore been achieved 
through closer placing of territories, whether allowed by decreased feeding range or 
greater tolerance to local differences in habitat. 
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Fic. 2.—A comparison of the distribution of feeding heights in three passerine species common 
to Bermuda and North America. Abscissae are in per cent time per relative leaf area per foot in 
height. Mean and standard error of plots observed are shown. 


Conclusion.— Absence of competition in Bermuda has allowed the few species 
present to attain far greater densities than they do in North America. Feeding 
behavior and habitat selection lie within the range of abilities exercised in North 
America. Distribution curves of feeding methods and feeding heights in Bermuda 
are significantly different from those of the mainland in the three species studied, 
but diversity of such curves is only slightly greater in one species. Adjacent species 


along critical ecological gradients therefore occupy broadly overlapping zones, 


rather than discrete segments; and limitation of the use of innate abilities is the 
result of the particular competitive situation.® 


Gratitude is expressed to Dr. Robert H. Mac Arthur, who suggested this approach and aided in 
its progress. A fuller treatment of the subject is to be published as part of a doctoral dissertation. 


* Contribution No, 280 from the Bermuda Biological Station. Field work was supported by 
grants from the Frank M. Chapman Fund of the American Museum of Natural History and from 
the National Science Foundation administered through the Bermuda Biological Station. 
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ERRATA 


In the paper by Demetrios G. Magiros entitled “A Method for Defining Principal 
Modes of Nonlinear Systems Utilizing Infinite Determinants,” which appeared in 
volume 46, number 12, pages 1608-1611, the author wishes to make the following 
corrections: 

1. On page 1609, line 2, substitute “anchor spring” for ‘“‘coupling,”’ so that the 
phrase reads “with nonlinear anchor spring.”’ 

2. On page 1609, section 3, lines | and 2, for the phrases “A, and AK; the spring 
constants of the anchor springs, Ke K, + wx? that of the coupling,” substitute 
“K, and K; the spring constants of the coupling and the second anchor spring, 


K, = K, + ux?’ that of the first anchor spring, . . .”’ 


In the paper entitled “On the Mutagenic Effect of Alkylating Agents,”’ by Ekke- 
hard Bautz and Ernst Freese, which appeared in volume 46, number 12, pages 
1585-1594, the authors wish to have the sentence beginning on line 2 of page 1591 
read as follows: “The smallest amount of 7-ethyladenine that could have been de- 


tected is 0.04 micromoles,” instead of 0.4 micromoles. 


In the paper entitled “‘Regulation by Coliphage Lambda of the Expression of the 
Capacity to Synthesize a Sequence of Host Enzymes,”’ by Michael B. Yarmolinsky 


and Herbert Wiesmeyer, which appeared in volume 46, number 12, pages 1626-1645, 
the authors wish to make the following corrections: 

Page 1636, line 15, for reference 44 read 54; page 1640, line 3 from bottom, for 
reference 44 read 54; page 1644, line 8, for reference 11 read 19; and page 1644, 
line 9, for reference 44 read 45. 
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